AnimaL

Tutorial

Documentation

---

Documentation

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

---

OctreeParticleSystemEngine.cpp

Go to the documentation of this file.

//
// C++ Implementation: OctreeParticleSystemEngine
//
// Description:
//
//
// Author: Matthieu Nesme <>, (C) 2004
//
// Copyright: See COPYING file that comes with this distribution
//
//
#include "OctreeParticleSystemEngine.h"
#include "ConstrainedVertex.h"

namespace animal
{

namespace octree
{


template<typename t_V, typename t_R, typename t_M>
OctreeParticleSystemEngine<t_V,t_R,t_M>::OctreeParticleSystemEngine(Vec3d bboxMin, Vec3d bboxMax, MFVec3f *points, MFVec3f normals, unsigned int nMaxPointsPerCell)
        : Octree(bboxMin, bboxMax, points, normals, nMaxPointsPerCell)
        , ParticleSystemEngine<t_V,t_R,t_M>()
{
    x=&_particlesPositions;
    v=&_particlesVelocities;
    m=&_particlesInvMasses;

    _springfield.links=&_springs;
    _springfield.stiffnesses=&_stiffnesses;
    _springfield.restLengths=&_restLengths;
    addForce( &_springfield );

    gravity=Vec(0,0,0);
}

template<typename t_V, typename t_R, typename t_M>
OctreeParticleSystemEngine<t_V,t_R,t_M>::~OctreeParticleSystemEngine()
{}




template<typename t_V, typename t_R, typename t_M>
int OctreeParticleSystemEngine<t_V,t_R,t_M>::allLevelsPreUpdateVerticesPositions(Cell * tmpCell,int *indice)
{
    if(tmpCell->getData()._points.size() != 0 )
    {
        for(int i=0;i<8;++i)
        {
            tmpCell->vertex(i)->_mass += pow(8,tmpCell->getData()._depth+1);
        }

        //if( tmpCell->isLeaf() )
        {
            for(int i=0;i<8;++i)
            {
                int pointeur = reinterpret_cast<int>(tmpCell->vertex(i));

                if(_verticesMap.find(pointeur)==_verticesMap.end())
                {
                    _verticesMap[pointeur] = (*indice)++;
                    _particlesPositions.push_back( Vec( tmpCell->vertex(i)->getPosition()[0],tmpCell->vertex(i)->getPosition()[1],tmpCell->vertex(i)->getPosition()[2] ));
                    _particlesVelocities.push_back( Vec(0,0,0) );
                    _particlesInvMasses.push_back( tmpCell->vertex(i)->_mass );
                }
            }

            for(int i=0;i<8;++i)
            {
                tmpCell->vertex(i)->_mass = 0; //re-init mass to 0
                int pointeuri = reinterpret_cast<int>(tmpCell->vertex(i));
                for(int j=i+1;j<8;++j)
                {
                    int pointeurj = reinterpret_cast<int>(tmpCell->vertex(j));

                    _allLevelsSprings[tmpCell->getData()._depth].push_back(_verticesMap[pointeuri]);
                    _allLevelsSprings[tmpCell->getData()._depth].push_back(_verticesMap[pointeurj]);
                    _allLevelsStiffnesses[tmpCell->getData()._depth].push_back(1000/pow(2,tmpCell->getData()._depth));

                    Real longueur = 0;


                    int directions[8][3] = { -1, -1, -1,
                                             1,-1,-1,
                                             -1,1,-1,
                                             1,1,-1,
                                             -1,-1,1,
                                             1,-1,1,
                                             -1,1,1,
                                             1,1,1};



                    if( directions[i][0] * directions[j][0] < 0 )
                        longueur += (tmpCell->getData()._initialSize[0] * tmpCell->getData()._initialSize[0]);
                    if( directions[i][1] * directions[j][1] < 0 )
                        longueur += (tmpCell->getData()._initialSize[1] * tmpCell->getData()._initialSize[1]);
                    if( directions[i][2] * directions[j][2] < 0 )
                        longueur += (tmpCell->getData()._initialSize[2] * tmpCell->getData()._initialSize[2]);
                    _allLevelsRestLengths[tmpCell->getData()._depth].push_back(sqrt(longueur));
                }
            }
        }
//        else
        {
            for(int i=0;i<8;++i)
                preUpdateVerticesPositions(tmpCell->child(i),indice);
        }
    }
}








template<typename t_V, typename t_R, typename t_M>
int OctreeParticleSystemEngine<t_V,t_R,t_M>::preUpdateVerticesPositions(Cell * tmpCell,int *indice)
{
    if(tmpCell->getData()._points.size() != 0 )
    {
        for(int i=0;i<8;++i)
        {
            tmpCell->vertex(i)->_mass += pow(8,tmpCell->getData()._depth+1);
        }

        if( tmpCell->isLeaf() )
        {
            for(int i=0;i<8;++i)
            {
                int pointeur = reinterpret_cast<int>(tmpCell->vertex(i));

                if(_verticesMap.find(pointeur)==_verticesMap.end())
                {
                    _verticesMap[pointeur] = (*indice)++;
                    _particlesPositions.push_back( Vec( tmpCell->vertex(i)->getPosition()[0],tmpCell->vertex(i)->getPosition()[1],tmpCell->vertex(i)->getPosition()[2] ));
                    _particlesVelocities.push_back( Vec(0,0,0) );
                    _particlesInvMasses.push_back( tmpCell->vertex(i)->_mass );
                }
            }

            for(int i=0;i<8;++i)
            {
                tmpCell->vertex(i)->_mass = 0; //re-init mass to 0
                int pointeuri = reinterpret_cast<int>(tmpCell->vertex(i));
                for(int j=i+1;j<8;++j)
                {
                    int pointeurj = reinterpret_cast<int>(tmpCell->vertex(j));

                    _springs.push_back(_verticesMap[pointeuri]);
                    _springs.push_back(_verticesMap[pointeurj]);
                    _stiffnesses.push_back(1000/pow(2,tmpCell->getData()._depth));

                    Real longueur = 0;


                    int directions[8][3] = { -1, -1, -1,
                                             1,-1,-1,
                                             -1,1,-1,
                                             1,1,-1,
                                             -1,-1,1,
                                             1,-1,1,
                                             -1,1,1,
                                             1,1,1};



                    if( directions[i][0] * directions[j][0] < 0 )
                        longueur += (tmpCell->getData()._initialSize[0] * tmpCell->getData()._initialSize[0]);
                    if( directions[i][1] * directions[j][1] < 0 )
                        longueur += (tmpCell->getData()._initialSize[1] * tmpCell->getData()._initialSize[1]);
                    if( directions[i][2] * directions[j][2] < 0 )
                        longueur += (tmpCell->getData()._initialSize[2] * tmpCell->getData()._initialSize[2]);
                    _restLengths.push_back(sqrt(longueur));
                }
            }
        }
        else
        {
            for(int i=0;i<8;++i)
                preUpdateVerticesPositions(tmpCell->child(i),indice);
        }
    }
}


template<typename t_V, typename t_R, typename t_M>
int OctreeParticleSystemEngine<t_V,t_R,t_M>::postUpdateVerticesPositions()
{
    for( VerticesMap::iterator it( _verticesMap.begin() ) ; it!=_verticesMap.end() ; ++it)
    {
        ConstrainedVertex * cv = reinterpret_cast<ConstrainedVertex *>((*it).first);
        cv->setPosition(Vec3d( _particlesPositions[(*it).second][0], _particlesPositions[(*it).second][1], _particlesPositions[(*it).second][2]  ));
        cv->_mass=0;
    }
}

template<typename t_V, typename t_R, typename t_M>
void OctreeParticleSystemEngine<t_V,t_R,t_M>::move(double dt)
{
    //cerr<<"--------move----------\n";
    _particlesPositions.resize(0);
    _particlesVelocities.resize(0);
    _particlesInvMasses.resize(0);
    _springs.resize(0);
    _stiffnesses.resize(0);
    _restLengths.resize(0);
    _verticesMap.clear();
    
    
    
/*  int n=10;
    _allLevelsSprings.resize(n);
    for(int i=0;i<n;++i)
    {
        _allLevelsSprings[i].resize(0);
        _allLevelsStiffnesses[i].resize(0);
        _allLevelsRestLengths[i].resize(0);
    }
    
    allLevelsPreUpdateVerticesPositions( root() );
    
    for(int i=0;i<n;++i)
    {
        if( _allLevelsSprings[i].size()!=0)
        {
            _springfield.links=&_allLevelsSprings[i];
            _springfield.stiffnesses=&_allLevelsStiffnesses[i];
            _springfield.restLengths=&_allLevelsRestLengths[i];
            ParticleSystemEngine<t_V,t_R,t_M>::move(dt);
            postUpdateVerticesPositions(  );
        }
        else
            break;
    }
*/
    
    
    preUpdateVerticesPositions( root() );
    ParticleSystemEngine<t_V,t_R,t_M>::move(dt);
    postUpdateVerticesPositions(  );

    
    movedRecLeaves(root());
}








// template<typename t_V, typename t_R, typename t_M>
// void OctreeParticleSystemEngine<t_V,t_R,t_M>::subdivideRecDeformed( Cell *cell )
// {
//  //cerr<<"OctreeParticleSystemEngine<t_V,t_R,t_M>::subdivideRecDeformed( Cell *cell )\n";
//     /*for( dfs_leaf_iterator it( cell ) ; !it.empty() ; )
//     {
//         Cell *tmpCell = ++it;
//
//         if( tmpCell->isLeaf() && tmpCell->getData()._points.size() != 0 )
//         {
//          subdivideDeformed( tmpCell );
//         }
//     }*/
//
//
//  subdivideDeformed( cell );
//
//  structureChangedRecLeaves(root());
//     movedRecLeaves(root());
// }

/*template<typename t_V, typename t_R, typename t_M>
void OctreeParticleSystemEngine<t_V,t_R,t_M>::subdivideDeformed( Cell *cell )
{
    //cerr<<"OctreeParticleSystemEngine<t_V,t_R,t_M>::subdivideDeformed( Cell *cell )\n";
    Octree::subdivideDeformed( cell );
    for(int i=0;i<8;++i)
        if(cell->child(i)->getData()._points.size() != 0)
            addSprings( cell->child(i) );
}*/

/*template<typename t_V, typename t_R, typename t_M>
void OctreeParticleSystemEngine<t_V,t_R,t_M>::subdivideDeformed( Cell *cell )
{
    if(cell)
    {
        if(cell->isLeaf())
        {
            if( cell->getData()._points.size() != 0 )
            {
                Octree::subdivideDeformed( cell );
                for(int i=0;i<8;++i)
                    if(cell->child(i)->getData()._points.size() != 0)
                        addSprings( cell->child(i) );
            }
        }
        else
        {
            for(int i=0;i<8;++i)
                subdivideDeformed(cell->child(i));
        }
    }
}*/









/*template<typename t_V, typename t_R, typename t_M>
void OctreeParticleSystemEngine<t_V,t_R,t_M>::addSprings( Cell *cell )
{
    for(int i=0;i<8;i++)
    {
        _particlesPositions.push_back( Vec( cell->vertex(i)->getPosition()[0],cell->vertex(i)->getPosition()[1],cell->vertex(i)->getPosition()[2]) );
        _particlesVelocities.push_back( Vec(0,0,0) );
        _particlesInvMasses.push_back(1);
    }   
    
    for(unsigned int i=_particlesPositions.size()-8;i<_particlesPositions.size();i++)
    {
        for(unsigned int j=i+1;j<_particlesPositions.size();++j)
        {
            _springs.push_back(i);
            _springs.push_back(j);
            _stiffnesses.push_back(1000);
            _restLengths.push_back(v_norm(_particlesPositions[i]-_particlesPositions[j]));
        }
    }
}*/


/*template<typename t_V, typename t_R, typename t_M>
void OctreeParticleSystemEngine<t_V,t_R,t_M>::simplify( Cell *cell )
{
    Octree::simplify( cell );
    
 
    _particlesPositions.resize(0);
    _particlesVelocities.resize(0);
    _particlesInvMasses.resize(0);
    
    _springs.resize(0);
    _stiffnesses.resize(0);
    _restLengths.resize(0);
    
    addSprings(root());
    
}*/



template<typename t_V, typename t_R, typename t_M>
void OctreeParticleSystemEngine<t_V,t_R,t_M>::draw()
{
    ParticleSystemEngine<t_V,t_R,t_M>::draw();
}


};

};

---