AnimaL

Tutorial

Documentation

---

Documentation

• Overview
• Modules
• Namespace List
• Class Hierarchy
• Alphabetical List
• Compound List
• File List
• Namespace Members
• Compound Members
• File Members
• Examples

---

octree_instanciation.cpp

Go to the documentation of this file.

#include <octree_instanciation.h>

namespace animal
{
namespace octree
{


Vec3d HVecToVec3d( const HVec v )
{
    Require( v[3] != 0 );

    Vec3d res;
    res[0] = v[0];
    res[1] = v[1];
    res[2] = v[2];
    res /= v[3];

    return res;
}

HVec Vec3dToHVec( const Vec3d v )
{
    HVec res(4);

    res[0] = v[0];
    res[1] = v[1];
    res[2] = v[2];
    res[3] = 1;

    return res;
}


FloatingPointType *getMatrixInverse( unsigned short size, FloatingPointType *data )
{
    HMat A(size,size);

    for( unsigned int j=0 ; j<size ; ++j )
    {
        for(unsigned int i=0 ; i<size ; ++i )
        {
            A[j][i] = data[i+j*size];
        }
    }

    HMat LU(size,size), LU_1(size,size);
    m_eq_ludcmp( LU, A );

    HVec aux(size);
    m_eq_luinv( LU_1, LU, aux );

    FloatingPointType *res = new FloatingPointType(size*size);

    for( unsigned int j=0 ; j<size ; ++j )
    {
        for(unsigned int i=0 ; i<size ; ++i )
        {
            res[i+j*size] = LU_1[j][i];
        }
    }

    return res;
}

HMat getMatrixInverse( HMat A )
{
    Require( A.nrows() == A.ncols() );

    unsigned short size = A.ncols();
    HMat LU(size,size), LU_1(size,size);
    m_eq_ludcmp( LU, A );

    HVec aux(size);
    m_eq_luinv( LU_1, LU, aux );

    return LU_1;
}

}
}

---