KisTridiagonalSystem< T > Class Template Reference

#include <kis_cubic_curve_spline.h>

Static Public Member Functions
static QVector< T >	calculate (QList< T > &a, QList< T > &b, QList< T > &c, QList< T > &f)

Detailed Description

template<typename T>
class KisTridiagonalSystem< T >

Definition at line 21 of file kis_cubic_curve_spline.h.

Member Function Documentation

calculate()

template<typename T >

static QVector< T > KisTridiagonalSystem< T >::calculate ( QList< T > & a, QList< T > & b, QList< T > & c, QList< T > & f )

inlinestatic

Returns

- vector that is storing x[]

Check for special case when order of the matrix is equal to 1

Common case

Definition at line 38 of file kis_cubic_curve_spline.h.

                                      {
        QVector<T> x;
        QVector<T> alpha;
        QVector<T> beta;
 
        int i;
        int size = b.size();
 
        Q_ASSERT(a.size() == size - 1 &&
                 c.size() == size - 1 &&
                 f.size() == size);
 
        x.resize(size);
 
        if (size == 1) {
            x[0] = f[0] / b[0];
            return x;
        }
 
        alpha.resize(size);
        beta.resize(size);
 
 
        alpha[1] = -c[0] / b[0];
        beta[1] =  f[0] / b[0];
 
        for (i = 1; i < size - 1; i++) {
            alpha[i+1] = -c[i] /
                         (a[i-1] * alpha[i] + b[i]);
 
            beta[i+1] = (f[i] - a[i-1] * beta[i])
                        /
                        (a[i-1] * alpha[i] + b[i]);
        }
 
        x.last() = (f.last() - a.last() * beta.last())
                   /
                   (b.last() + a.last() * alpha.last());
 
        for (i = size - 2; i >= 0; i--)
            x[i] = alpha[i+1] * x[i+1] + beta[i+1];
 
        return x;
    }

---

The documentation for this class was generated from the following file:

• libs/image/kis_cubic_curve_spline.h