LinearApproximation.js 2.8 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071
  1. import defined from './defined.js';
  2. import DeveloperError from './DeveloperError.js';
  3. /**
  4. * An {@link InterpolationAlgorithm} for performing linear interpolation.
  5. *
  6. * @exports LinearApproximation
  7. */
  8. var LinearApproximation = {
  9. type : 'Linear'
  10. };
  11. /**
  12. * Given the desired degree, returns the number of data points required for interpolation.
  13. * Since linear interpolation can only generate a first degree polynomial, this function
  14. * always returns 2.
  15. * @param {Number} degree The desired degree of interpolation.
  16. * @returns {Number} This function always returns 2.
  17. *
  18. */
  19. LinearApproximation.getRequiredDataPoints = function(degree) {
  20. return 2;
  21. };
  22. /**
  23. * Interpolates values using linear approximation.
  24. *
  25. * @param {Number} x The independent variable for which the dependent variables will be interpolated.
  26. * @param {Number[]} xTable The array of independent variables to use to interpolate. The values
  27. * in this array must be in increasing order and the same value must not occur twice in the array.
  28. * @param {Number[]} yTable The array of dependent variables to use to interpolate. For a set of three
  29. * dependent values (p,q,w) at time 1 and time 2 this should be as follows: {p1, q1, w1, p2, q2, w2}.
  30. * @param {Number} yStride The number of dependent variable values in yTable corresponding to
  31. * each independent variable value in xTable.
  32. * @param {Number[]} [result] An existing array into which to store the result.
  33. * @returns {Number[]} The array of interpolated values, or the result parameter if one was provided.
  34. */
  35. LinearApproximation.interpolateOrderZero = function(x, xTable, yTable, yStride, result) {
  36. //>>includeStart('debug', pragmas.debug);
  37. if (xTable.length !== 2) {
  38. throw new DeveloperError('The xTable provided to the linear interpolator must have exactly two elements.');
  39. } else if (yStride <= 0) {
  40. throw new DeveloperError('There must be at least 1 dependent variable for each independent variable.');
  41. }
  42. //>>includeEnd('debug');
  43. if (!defined(result)) {
  44. result = new Array(yStride);
  45. }
  46. var i;
  47. var y0;
  48. var y1;
  49. var x0 = xTable[0];
  50. var x1 = xTable[1];
  51. //>>includeStart('debug', pragmas.debug);
  52. if (x0 === x1) {
  53. throw new DeveloperError('Divide by zero error: xTable[0] and xTable[1] are equal');
  54. }
  55. //>>includeEnd('debug');
  56. for (i = 0; i < yStride; i++) {
  57. y0 = yTable[i];
  58. y1 = yTable[i + yStride];
  59. result[i] = (((y1 - y0) * x) + (x1 * y0) - (x0 * y1)) / (x1 - x0);
  60. }
  61. return result;
  62. };
  63. export default LinearApproximation;