Function Tridiag_Solve

src/SB/Core/x/xSpline.cpp:14–96

static F32 sBasisHermite[4][4];

void Tridiag_Solve(F32* a, F32* b, F32* c, xVec3* d, xVec3* x, S32 n)
{
    S32 j;

    F32* a_temp;
    F32* b_temp;
    F32* c_temp;
    xVec3* delta;

    S32 vec_offset;
    F32 beta;
    F32* gamma;

    F32* c_prime;
    F32* d_prime;

    c_prime = (F32*)RwMalloc(n << 2);
    d_prime = (F32*)RwMalloc(n * 0xC);

    c_prime[0] = *c / *b;

    d_prime[0] = d->x / *b;
    d_prime[1] = d->y / *b;
    d_prime[2] = d->z / *b;

    a_temp = a + 1;
    b_temp = b + 1;
    c_temp = c + 1;
    delta = d + 1;
    vec_offset = 0xC;

    if (n > 1)
    {
        for (j = 1; j < n; j += 1)
        {
            beta = b_temp[0];
            gamma = (F32*)((S32)d_prime + vec_offset - 0xc);
            vec_offset = vec_offset + 0xC;
            b_temp = b_temp + 1;
            c_temp = c_temp + 1;
            *b_temp = *b_temp - *a_temp * c_prime[0];
            c_prime[1] = *c_temp / beta;
            d_prime[3] = (delta->x - *a_temp * gamma[0]) / beta;
            d_prime[4] = (delta->y - *a_temp * gamma[1]) / beta;
            *c_temp = *a_temp;
            a_temp = a_temp + 1;
            delta = delta + 1;
            d_prime[5] = (*(&delta->z) - *c_temp * gamma[2]) / beta;
            c_prime = c_prime + 1;
            d_prime = d_prime + 3;
        }
    }

    j = n - 2;
    c_prime = d_prime + (n - 1) * 3;
    delta->x = c_prime[0];
    delta->y = c_prime[1];
    delta->z = c_prime[2];