Changeset 2736

Timestamp:

01/18/13 10:53:15 (12 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes CAS-4770

Ready for Test: Yes

Interface Changes: Yes/No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: Yes/No

Module(s): Module Names change impacts.

Description: Describe your changes here...

Bug fix on reusable_ attribute.

Location:

trunk/src

Files:

• CubicSplineInterpolator1D.h (modified) (1 diff)
• CubicSplineInterpolator1D.tcc (modified) (2 diffs)

trunk/src/CubicSplineInterpolator1D.h

@@ -33 +33 @@
   // @see Interpolator1D::setData
   void setData(T *x, U *y, unsigned int n);
+
+  // Override Interpolate1D::setX.
+  // @see Interpolate1D::setX()
+  void setX(T *x, unsigned int n);
 
   // Override Interpolator1D::setY.

trunk/src/CubicSplineInterpolator1D.tcc

@@ -38 +38 @@
 {
   Interpolator1D<T, U>::setData(x, y, n);
+  reusable_ = false;
+}
+
+template <class T, class U>
+void CubicSplineInterpolator1D<T, U>::setX(T *x, unsigned int n)
+{
+  Interpolator1D<T, U>::setX(x, n);
   reusable_ = false;
 }
@@ -94 +101 @@
 
   // Solve tridiagonal system.
-  // Here, tridiagonal matrix is decomposed to triangular matrix
-  // u stores upper triangular components while y2_ stores 
-  // right-hand side vector.
+  // Here, tridiagonal matrix is decomposed to upper triangular 
+  // matrix. u stores upper triangular components while y2_ stores 
+  // right-hand side vector. The diagonal elements are normalized 
+  // to 1.
   T a1 = this->x_[1] - this->x_[0];
   T a2, bi;