CcdPhysicsController.cpp

Go to the documentation of this file.

/*
Bullet Continuous Collision Detection and Physics Library
Copyright (c) 2003-2006 Erwin Coumans  http://continuousphysics.com/Bullet/

This software is provided 'as-is', without any express or implied warranty.
In no event will the authors be held liable for any damages arising from the use of this software.
Permission is granted to anyone to use this software for any purpose, 
including commercial applications, and to alter it and redistribute it freely, 
subject to the following restrictions:

1. The origin of this software must not be misrepresented; you must not claim that you wrote the original software. If you use this software in a product, an acknowledgment in the product documentation would be appreciated but is not required.
2. Altered source versions must be plainly marked as such, and must not be misrepresented as being the original software.
3. This notice may not be removed or altered from any source distribution.
*/

#ifndef WIN32
#include <stdint.h>
#endif

#include "CcdPhysicsController.h"
#include "btBulletDynamicsCommon.h"
#include "BulletCollision/CollisionShapes/btScaledBvhTriangleMeshShape.h"

#include "BulletCollision/CollisionShapes/btTriangleIndexVertexArray.h"

#include "PHY_IMotionState.h"
#include "CcdPhysicsEnvironment.h"
#include "RAS_MeshObject.h"
#include "KX_GameObject.h"

#include "BulletSoftBody/btSoftBody.h"
#include "BulletSoftBody//btSoftBodyInternals.h"
#include "BulletSoftBody/btSoftBodyHelpers.h"
#include "LinearMath/btConvexHull.h"
#include "BulletCollision/Gimpact/btGImpactShape.h"
#include "BulletCollision/Gimpact/btGImpactShape.h"


#include "BulletSoftBody/btSoftRigidDynamicsWorld.h"

#include "DNA_mesh_types.h"
#include "DNA_meshdata_types.h"

extern "C"{
#include "BKE_cdderivedmesh.h"
}

class BP_Proxy;


//'temporarily' global variables
//float gDeactivationTime = 2.f;
//bool  gDisableDeactivation = false;
extern float gDeactivationTime;
extern bool gDisableDeactivation;


float gLinearSleepingTreshold = 0.8f;
float gAngularSleepingTreshold = 1.0f;


btVector3 startVel(0,0,0);//-10000);

CcdPhysicsController::CcdPhysicsController (const CcdConstructionInfo& ci)
:m_cci(ci)
{
        m_prototypeTransformInitialized = false;
        m_softbodyMappingDone = false;
        m_collisionDelay = 0;
        m_newClientInfo = 0;
        m_registerCount = 0;
        m_softBodyTransformInitialized = false;
        m_parentCtrl = 0;
        // copy pointers locally to allow smart release
        m_MotionState = ci.m_MotionState;
        m_collisionShape = ci.m_collisionShape;
        // apply scaling before creating rigid body
        m_collisionShape->setLocalScaling(m_cci.m_scaling);
        if (m_cci.m_mass)
                m_collisionShape->calculateLocalInertia(m_cci.m_mass, m_cci.m_localInertiaTensor);
        // shape info is shared, increment ref count
        m_shapeInfo = ci.m_shapeInfo;
        if (m_shapeInfo)
                m_shapeInfo->AddRef();
        
        m_bulletMotionState = 0;
        
        
        CreateRigidbody();
        

/*#ifdef WIN32
        if (GetRigidBody() && !GetRigidBody()->isStaticObject())
                GetRigidBody()->setLinearVelocity(startVel);
#endif*/

}

btTransform&    CcdPhysicsController::GetTransformFromMotionState(PHY_IMotionState* motionState)
{
        static btTransform trans;
        btVector3 tmp;
        motionState->getWorldPosition(tmp.m_floats[0], tmp.m_floats[1], tmp.m_floats[2]);
        trans.setOrigin(tmp);

        float ori[12];
        motionState->getWorldOrientation(ori);
        trans.getBasis().setFromOpenGLSubMatrix(ori);
        //btQuaternion orn;
        //motionState->getWorldOrientation(orn[0],orn[1],orn[2],orn[3]);
        //trans.setRotation(orn);
        return trans;

}

class   BlenderBulletMotionState : public btMotionState
{
        PHY_IMotionState*       m_blenderMotionState;

public:

        BlenderBulletMotionState(PHY_IMotionState* bms)
                :m_blenderMotionState(bms)
        {

        }

        void    getWorldTransform(btTransform& worldTrans ) const
        {
                btVector3 pos;
                float ori[12];

                m_blenderMotionState->getWorldPosition(pos.m_floats[0],pos.m_floats[1],pos.m_floats[2]);
                m_blenderMotionState->getWorldOrientation(ori);
                worldTrans.setOrigin(pos);
                worldTrans.getBasis().setFromOpenGLSubMatrix(ori);
        }

        void    setWorldTransform(const btTransform& worldTrans)
        {
                m_blenderMotionState->setWorldPosition(worldTrans.getOrigin().getX(),worldTrans.getOrigin().getY(),worldTrans.getOrigin().getZ());
                btQuaternion rotQuat = worldTrans.getRotation();
                m_blenderMotionState->setWorldOrientation(rotQuat[0],rotQuat[1],rotQuat[2],rotQuat[3]);
                m_blenderMotionState->calculateWorldTransformations();
        }

};


btRigidBody* CcdPhysicsController::GetRigidBody()
{
        return btRigidBody::upcast(m_object);
}
btCollisionObject*      CcdPhysicsController::GetCollisionObject()
{
        return m_object;
}
btSoftBody* CcdPhysicsController::GetSoftBody()
{
        return btSoftBody::upcast(m_object);
}

#include "BulletSoftBody/btSoftBodyHelpers.h"


bool CcdPhysicsController::CreateSoftbody()
{
        int shapeType = m_cci.m_collisionShape ? m_cci.m_collisionShape->getShapeType() : 0;

        //disable soft body until first sneak preview is ready
        if (!m_cci.m_bSoft || !m_cci.m_collisionShape ||
                ((shapeType != CONVEX_HULL_SHAPE_PROXYTYPE)&&
                (shapeType != TRIANGLE_MESH_SHAPE_PROXYTYPE) &&
                (shapeType != SCALED_TRIANGLE_MESH_SHAPE_PROXYTYPE)))
        {
                return false;
        }

        btRigidBody::btRigidBodyConstructionInfo rbci(m_cci.m_mass,m_bulletMotionState,m_collisionShape,m_cci.m_localInertiaTensor * m_cci.m_inertiaFactor);
        rbci.m_linearDamping = m_cci.m_linearDamping;
        rbci.m_angularDamping = m_cci.m_angularDamping;
        rbci.m_friction = m_cci.m_friction;
        rbci.m_restitution = m_cci.m_restitution;
        
        btVector3 p(0,0,0);// = getOrigin();
        //btSoftBody*   psb=btSoftBodyHelpers::CreateRope(worldInfo,    btVector3(-10,0,i*0.25),btVector3(10,0,i*0.25), 16,1+2);
        btSoftBody* psb  = 0;
        btSoftBodyWorldInfo& worldInfo = m_cci.m_physicsEnv->getDynamicsWorld()->getWorldInfo();

        if (m_cci.m_collisionShape->getShapeType() == CONVEX_HULL_SHAPE_PROXYTYPE)
        {
                btConvexHullShape* convexHull = (btConvexHullShape* )m_cci.m_collisionShape;
                {
                        int nvertices = convexHull->getNumPoints();
                        const btVector3* vertices = convexHull->getPoints();

                        HullDesc                hdsc(QF_TRIANGLES,nvertices,vertices);
                        HullResult              hres;
                        HullLibrary             hlib;/*??*/ 
                        hdsc.mMaxVertices=nvertices;
                        hlib.CreateConvexHull(hdsc,hres);
                        
                        psb=new btSoftBody(&worldInfo,(int)hres.mNumOutputVertices,
                                &hres.m_OutputVertices[0],0);
                        for(int i=0;i<(int)hres.mNumFaces;++i)
                        {
                                const int idx[]={       hres.m_Indices[i*3+0],
                                        hres.m_Indices[i*3+1],
                                        hres.m_Indices[i*3+2]};
                                if(idx[0]<idx[1]) psb->appendLink(      idx[0],idx[1]);
                                if(idx[1]<idx[2]) psb->appendLink(      idx[1],idx[2]);
                                if(idx[2]<idx[0]) psb->appendLink(      idx[2],idx[0]);
                                psb->appendFace(idx[0],idx[1],idx[2]);
                        }
                        hlib.ReleaseResult(hres);
                }
        } else
        {
                int numtris = 0;
                if (m_cci.m_collisionShape->getShapeType() ==SCALED_TRIANGLE_MESH_SHAPE_PROXYTYPE)
                {
                        btScaledBvhTriangleMeshShape* scaledtrimeshshape = (btScaledBvhTriangleMeshShape*) m_cci.m_collisionShape;
                        btBvhTriangleMeshShape* trimeshshape = scaledtrimeshshape->getChildShape();

                        if (trimeshshape->getMeshInterface()->getNumSubParts()==1)
                        {
                                unsigned char* vertexBase;
                                PHY_ScalarType vertexType;
                                int numverts;
                                int vertexstride;
                                unsigned char* indexbase;
                                int indexstride;
                                PHY_ScalarType indexType;
                                trimeshshape->getMeshInterface()->getLockedVertexIndexBase(&vertexBase,numverts,vertexType,vertexstride,&indexbase,indexstride,numtris,indexType);
                                
                                psb = btSoftBodyHelpers::CreateFromTriMesh(worldInfo,(const btScalar*)vertexBase,(const int*)indexbase,numtris,false);
                        }
                } else
                {
                        btTriangleMeshShape* trimeshshape = (btTriangleMeshShape*) m_cci.m_collisionShape;
                        if (trimeshshape->getMeshInterface()->getNumSubParts()==1)
                        {
                                unsigned char* vertexBase;
                                PHY_ScalarType vertexType;
                                int numverts;
                                int vertexstride;
                                unsigned char* indexbase;
                                int indexstride;
                                PHY_ScalarType indexType;
                                trimeshshape->getMeshInterface()->getLockedVertexIndexBase(&vertexBase,numverts,vertexType,vertexstride,&indexbase,indexstride,numtris,indexType);
                                
                                psb = btSoftBodyHelpers::CreateFromTriMesh(worldInfo,(const btScalar*)vertexBase,(const int*)indexbase,numtris,false);
                        }
                }
                // store face tag so that we can find our original face when doing ray casting
                btSoftBody::Face* ft;
                int i;
                for (i=0, ft=&psb->m_faces[0]; i<numtris; ++i, ++ft)
                {
                        // Hack!! use m_tag to store the face number, normally it is a pointer
                        // add 1 to make sure it is never 0
                        ft->m_tag = (void*)((uintptr_t)(i+1));
                }
        }
        if (m_cci.m_margin > 0.f)
        {
                psb->getCollisionShape()->setMargin(m_cci.m_margin);
                psb->updateBounds();
        }
        m_object = psb;
        
        //btSoftBody::Material* pm=psb->appendMaterial();
        btSoftBody::Material*   pm=psb->m_materials[0];
        pm->m_kLST                              =       m_cci.m_soft_linStiff;
        pm->m_kAST                              =       m_cci.m_soft_angStiff;
        pm->m_kVST                              =       m_cci.m_soft_volume;
        psb->m_cfg.collisions = 0;

        if (m_cci.m_soft_collisionflags & CCD_BSB_COL_CL_RS)
        {
                psb->m_cfg.collisions   +=      btSoftBody::fCollision::CL_RS;
        } else
        {
                psb->m_cfg.collisions   +=      btSoftBody::fCollision::SDF_RS;
        }
        if (m_cci.m_soft_collisionflags & CCD_BSB_COL_CL_SS)
        {
                psb->m_cfg.collisions += btSoftBody::fCollision::CL_SS;
        } else
        {
                psb->m_cfg.collisions += btSoftBody::fCollision::VF_SS;
        }


        psb->m_cfg.kSRHR_CL = m_cci.m_soft_kSRHR_CL;            /* Soft vs rigid hardness [0,1] (cluster only) */
        psb->m_cfg.kSKHR_CL = m_cci.m_soft_kSKHR_CL;            /* Soft vs kinetic hardness [0,1] (cluster only) */
        psb->m_cfg.kSSHR_CL = m_cci.m_soft_kSSHR_CL;            /* Soft vs soft hardness [0,1] (cluster only) */
        psb->m_cfg.kSR_SPLT_CL = m_cci.m_soft_kSR_SPLT_CL;      /* Soft vs rigid impulse split [0,1] (cluster only) */

        psb->m_cfg.kSK_SPLT_CL = m_cci.m_soft_kSK_SPLT_CL;      /* Soft vs rigid impulse split [0,1] (cluster only) */
        psb->m_cfg.kSS_SPLT_CL = m_cci.m_soft_kSS_SPLT_CL;      /* Soft vs rigid impulse split [0,1] (cluster only) */
        psb->m_cfg.kVCF = m_cci.m_soft_kVCF;                    /* Velocities correction factor (Baumgarte) */
        psb->m_cfg.kDP = m_cci.m_soft_kDP;                      /* Damping coefficient [0,1] */

        psb->m_cfg.kDG = m_cci.m_soft_kDG;                      /* Drag coefficient [0,+inf] */
        psb->m_cfg.kLF = m_cci.m_soft_kLF;                      /* Lift coefficient [0,+inf] */
        psb->m_cfg.kPR = m_cci.m_soft_kPR;                      /* Pressure coefficient [-inf,+inf] */
        psb->m_cfg.kVC = m_cci.m_soft_kVC;                      /* Volume conversation coefficient [0,+inf] */

        psb->m_cfg.kDF = m_cci.m_soft_kDF;                      /* Dynamic friction coefficient [0,1] */
        psb->m_cfg.kMT = m_cci.m_soft_kMT;                      /* Pose matching coefficient [0,1] */
        psb->m_cfg.kCHR = m_cci.m_soft_kCHR;                    /* Rigid contacts hardness [0,1] */
        psb->m_cfg.kKHR = m_cci.m_soft_kKHR;                    /* Kinetic contacts hardness [0,1] */

        psb->m_cfg.kSHR = m_cci.m_soft_kSHR;                    /* Soft contacts hardness [0,1] */
        psb->m_cfg.kAHR = m_cci.m_soft_kAHR;                    /* Anchors hardness [0,1] */

        if (m_cci.m_gamesoftFlag & CCD_BSB_BENDING_CONSTRAINTS)//OB_SB_GOAL)
        {
                psb->generateBendingConstraints(2,pm);
        }

        psb->m_cfg.piterations = m_cci.m_soft_piterations;
        psb->m_cfg.viterations = m_cci.m_soft_viterations;
        psb->m_cfg.diterations = m_cci.m_soft_diterations;
        psb->m_cfg.citerations = m_cci.m_soft_citerations;

        if (m_cci.m_gamesoftFlag & CCD_BSB_SHAPE_MATCHING)//OB_SB_GOAL)
        {
                psb->setPose(false,true);//
        } else
        {
                psb->setPose(true,false);
        }
        
        psb->randomizeConstraints();

        if (m_cci.m_soft_collisionflags & (CCD_BSB_COL_CL_RS+CCD_BSB_COL_CL_SS))
        {
                psb->generateClusters(m_cci.m_soft_numclusteriterations);
        }

        psb->setTotalMass(m_cci.m_mass);
        
        psb->setCollisionFlags(0);

        {
                RAS_MeshObject* rasMesh= GetShapeInfo()->GetMesh();

                if (rasMesh && !m_softbodyMappingDone)
                {
                        //printf("apply\n");
                        RAS_MeshSlot::iterator it;
                        RAS_MeshMaterial *mmat;
                        RAS_MeshSlot *slot;
                        size_t i;

                        //for each material
                        for (int m=0;m<rasMesh->NumMaterials();m++)
                        {
                                mmat = rasMesh->GetMeshMaterial(m);

                                slot = mmat->m_baseslot;
                                for(slot->begin(it); !slot->end(it); slot->next(it))
                                {
                                        int index = 0;
                                        for(i=it.startvertex; i<it.endvertex; i++,index++) 
                                        {
                                                RAS_TexVert* vertex = &it.vertex[i];
                                                //search closest index, and store it in vertex
                                                vertex->setSoftBodyIndex(0);
                                                btScalar maxDistSqr = 1e30;
                                                btSoftBody::tNodeArray&   nodes(psb->m_nodes);
                                                btVector3 xyz = btVector3(vertex->getXYZ()[0],vertex->getXYZ()[1],vertex->getXYZ()[2]);
                                                for (int n=0;n<nodes.size();n++)
                                                {
                                                        btScalar distSqr = (nodes[n].m_x - xyz).length2();
                                                        if (distSqr<maxDistSqr)
                                                        {
                                                                maxDistSqr = distSqr;
                                                                
                                                                vertex->setSoftBodyIndex(n);
                                                        }
                                                }
                                        }
                                }
                        }
                }
        }
        m_softbodyMappingDone = true;

        btTransform startTrans;
        rbci.m_motionState->getWorldTransform(startTrans);

        m_MotionState->setWorldPosition(startTrans.getOrigin().getX(),startTrans.getOrigin().getY(),startTrans.getOrigin().getZ());
        m_MotionState->setWorldOrientation(0,0,0,1);

        if (!m_prototypeTransformInitialized)
        {
                m_prototypeTransformInitialized = true;
                m_softBodyTransformInitialized = true;
                psb->transform(startTrans);
        }
        m_object->setCollisionFlags(m_object->getCollisionFlags() | m_cci.m_collisionFlags);
        if (m_cci.m_do_anisotropic)
                m_object->setAnisotropicFriction(m_cci.m_anisotropicFriction);
        return true;
}


void CcdPhysicsController::CreateRigidbody()
{

        //btTransform trans = GetTransformFromMotionState(m_MotionState);
        m_bulletMotionState = new BlenderBulletMotionState(m_MotionState);

        if (CreateSoftbody())
                // soft body created, done
                return;

        //create a rgid collision object
        btRigidBody::btRigidBodyConstructionInfo rbci(m_cci.m_mass,m_bulletMotionState,m_collisionShape,m_cci.m_localInertiaTensor * m_cci.m_inertiaFactor);
        rbci.m_linearDamping = m_cci.m_linearDamping;
        rbci.m_angularDamping = m_cci.m_angularDamping;
        rbci.m_friction = m_cci.m_friction;
        rbci.m_restitution = m_cci.m_restitution;
        m_object = new btRigidBody(rbci);
        
        //
        // init the rigidbody properly
        //
        
        //setMassProps this also sets collisionFlags
        //convert collision flags!
        //special case: a near/radar sensor controller should not be defined static or it will
        //generate loads of static-static collision messages on the console
        if (m_cci.m_bSensor)
        {
                // reset the flags that have been set so far
                GetCollisionObject()->setCollisionFlags(0);
                // sensor must never go to sleep: they need to detect continously
                GetCollisionObject()->setActivationState(DISABLE_DEACTIVATION);
        }
        GetCollisionObject()->setCollisionFlags(m_object->getCollisionFlags() | m_cci.m_collisionFlags);
        btRigidBody* body = GetRigidBody();

        if (body)
        {
                body->setGravity( m_cci.m_gravity);
                body->setDamping(m_cci.m_linearDamping, m_cci.m_angularDamping);

                if (!m_cci.m_bRigid)
                {
                        body->setAngularFactor(0.f);
                }
                body->setContactProcessingThreshold(m_cci.m_contactProcessingThreshold);

        }
        if (m_object && m_cci.m_do_anisotropic)
        {
                m_object->setAnisotropicFriction(m_cci.m_anisotropicFriction);
        }
                
}

static void DeleteBulletShape(btCollisionShape* shape, bool free)
{
        if (shape->getShapeType() == TRIANGLE_MESH_SHAPE_PROXYTYPE)
        {
                // shapes based on meshes use an interface that contains the vertices.
                btTriangleMeshShape* meshShape = static_cast<btTriangleMeshShape*>(shape);
                btStridingMeshInterface* meshInterface = meshShape->getMeshInterface();
                if (meshInterface)
                        delete meshInterface;
        }
        if(free) {
                delete shape;
        }
}

bool CcdPhysicsController::DeleteControllerShape( )
{
        if (m_collisionShape)
        {
                // collision shape is always unique to the controller, can delete it here
                if (m_collisionShape->isCompound())
                {
                        // bullet does not delete the child shape, must do it here
                        btCompoundShape* compoundShape = (btCompoundShape*)m_collisionShape;
                        int numChild = compoundShape->getNumChildShapes();
                        for (int i=numChild-1 ; i >= 0; i--)
                        {
                                btCollisionShape* childShape = compoundShape->getChildShape(i);
                                DeleteBulletShape(childShape, true);
                        }
                }
                DeleteBulletShape(m_collisionShape, true);

                return true;
        }

        return false;
}

bool CcdPhysicsController::ReplaceControllerShape(btCollisionShape *newShape)
{
        
        /* Note, deleting the previous collision shape must be done already */
        /* if (m_collisionShape) DeleteControllerShape(); */

        m_object->setCollisionShape(newShape);
        m_collisionShape= newShape;
        m_cci.m_collisionShape= newShape;
        
        if (GetSoftBody()) {
                // soft body must be recreated
                m_cci.m_physicsEnv->removeCcdPhysicsController(this);
                delete m_object;
                m_object = NULL;
                // force complete reinitialization
                m_softbodyMappingDone = false;
                m_prototypeTransformInitialized = false;
                m_softBodyTransformInitialized = false;
                CreateSoftbody();
                assert(m_object);
                // reinsert the new body
                m_cci.m_physicsEnv->addCcdPhysicsController(this);
        }
        
        /* Copied from CcdPhysicsEnvironment::addCcdPhysicsController() */
        
        /* without this, an object can rest on the old physics mesh
         * and not move to account for the physics mesh, even with 'nosleep' */ 
        btSoftRigidDynamicsWorld* dw= GetPhysicsEnvironment()->getDynamicsWorld();
        btCollisionObjectArray &obarr= dw->getCollisionObjectArray();
        btCollisionObject *ob;
        btBroadphaseProxy* proxy;

        for(int i= 0; i < obarr.size(); i++) {
                ob= obarr[i];
                if (ob->getCollisionShape() == newShape) {
                        proxy = ob->getBroadphaseHandle();
                        
                        if(proxy)
                                dw->getPairCache()->cleanProxyFromPairs(proxy,dw->getDispatcher());
                }
        }
        
        return true;
}

CcdPhysicsController::~CcdPhysicsController()
{
        //will be reference counted, due to sharing
        if (m_cci.m_physicsEnv)
                m_cci.m_physicsEnv->removeCcdPhysicsController(this);

        if (m_MotionState)
                delete m_MotionState;
        if (m_bulletMotionState)
                delete m_bulletMotionState;
        delete m_object;

        DeleteControllerShape();

        if (m_shapeInfo)
        {
                m_shapeInfo->Release();
        }
}


bool            CcdPhysicsController::SynchronizeMotionStates(float time)
{
        //sync non-static to motionstate, and static from motionstate (todo: add kinematic etc.)

        btSoftBody* sb = GetSoftBody();
        if (sb)
        {
                if (sb->m_pose.m_bframe) 
                {
                        btVector3 worldPos = sb->m_pose.m_com;
                        btQuaternion worldquat;
                        btMatrix3x3     trs = sb->m_pose.m_rot*sb->m_pose.m_scl;
                        trs.getRotation(worldquat);
                        m_MotionState->setWorldPosition(worldPos[0],worldPos[1],worldPos[2]);
                        m_MotionState->setWorldOrientation(worldquat[0],worldquat[1],worldquat[2],worldquat[3]);
                }
                else 
                {
                        btVector3 aabbMin,aabbMax;
                        sb->getAabb(aabbMin,aabbMax);
                        btVector3 worldPos  = (aabbMax+aabbMin)*0.5f;
                        m_MotionState->setWorldPosition(worldPos[0],worldPos[1],worldPos[2]);
                }
                m_MotionState->calculateWorldTransformations();
                return true;
        }

        btRigidBody* body = GetRigidBody();

        if (body && !body->isStaticObject())
        {
                
                if ((m_cci.m_clamp_vel_max>0.0) || (m_cci.m_clamp_vel_min>0.0))
                {
                        const btVector3& linvel = body->getLinearVelocity();
                        float len= linvel.length();
                        
                        if((m_cci.m_clamp_vel_max>0.0) && (len > m_cci.m_clamp_vel_max))
                                        body->setLinearVelocity(linvel * (m_cci.m_clamp_vel_max / len));
                        
                        else if ((m_cci.m_clamp_vel_min>0.0) && btFuzzyZero(len)==0 && (len < m_cci.m_clamp_vel_min))
                                body->setLinearVelocity(linvel * (m_cci.m_clamp_vel_min / len));
                }
                
                const btTransform& xform = body->getCenterOfMassTransform();
                const btMatrix3x3& worldOri = xform.getBasis();
                const btVector3& worldPos = xform.getOrigin();
                float ori[12];
                worldOri.getOpenGLSubMatrix(ori);
                m_MotionState->setWorldOrientation(ori);
                m_MotionState->setWorldPosition(worldPos[0],worldPos[1],worldPos[2]);
                m_MotionState->calculateWorldTransformations();

                float scale[3];
                m_MotionState->getWorldScaling(scale[0],scale[1],scale[2]);
                btVector3 scaling(scale[0],scale[1],scale[2]);
                GetCollisionShape()->setLocalScaling(scaling);
        } else
        {
                btVector3 worldPos;
                btQuaternion worldquat;

/*              m_MotionState->getWorldPosition(worldPos[0],worldPos[1],worldPos[2]);
                m_MotionState->getWorldOrientation(worldquat[0],worldquat[1],worldquat[2],worldquat[3]);
                btTransform oldTrans = m_body->getCenterOfMassTransform();
                btTransform newTrans(worldquat,worldPos);
                                
                SetCenterOfMassTransform(newTrans);
                //need to keep track of previous position for friction effects...
                
                m_MotionState->calculateWorldTransformations();
*/
                float scale[3];
                m_MotionState->getWorldScaling(scale[0],scale[1],scale[2]);
                btVector3 scaling(scale[0],scale[1],scale[2]);
                GetCollisionShape()->setLocalScaling(scaling);
        }
        return true;

}

void            CcdPhysicsController::WriteMotionStateToDynamics(bool nondynaonly)
{
        btTransform& xform = CcdPhysicsController::GetTransformFromMotionState(m_MotionState);
        SetCenterOfMassTransform(xform);
}

void            CcdPhysicsController::WriteDynamicsToMotionState()
{
}
                // controller replication
void            CcdPhysicsController::PostProcessReplica(class PHY_IMotionState* motionstate,class PHY_IPhysicsController* parentctrl)
{
        
        m_softBodyTransformInitialized=false;
        m_MotionState = motionstate;
        m_registerCount = 0;
        m_collisionShape = NULL;

        // always create a new shape to avoid scaling bug
        if (m_shapeInfo)
        {
                m_shapeInfo->AddRef();
                m_collisionShape = m_shapeInfo->CreateBulletShape(m_cci.m_margin, m_cci.m_bGimpact, !m_cci.m_bSoft);

                if (m_collisionShape)
                {
                        // new shape has no scaling, apply initial scaling
                        //m_collisionShape->setMargin(m_cci.m_margin);
                        m_collisionShape->setLocalScaling(m_cci.m_scaling);
                        
                        if (m_cci.m_mass)
                                m_collisionShape->calculateLocalInertia(m_cci.m_mass, m_cci.m_localInertiaTensor);
                }
        }

        // load some characterists that are not 
        btRigidBody* oldbody = GetRigidBody();
        m_object = 0;
        CreateRigidbody();
        btRigidBody* body = GetRigidBody();
        if (body)
        {
                if (m_cci.m_mass)
                {
                        body->setMassProps(m_cci.m_mass, m_cci.m_localInertiaTensor * m_cci.m_inertiaFactor);
                }

                if (oldbody)
                {
                        body->setLinearFactor(oldbody->getLinearFactor());
                        body->setAngularFactor(oldbody->getAngularFactor());
                        if (oldbody->getActivationState() == DISABLE_DEACTIVATION)
                                body->setActivationState(DISABLE_DEACTIVATION);
                }
        }       
        // sensor object are added when needed
        if (!m_cci.m_bSensor)
                m_cci.m_physicsEnv->addCcdPhysicsController(this);


/*      SM_Object* dynaparent=0;
        SumoPhysicsController* sumoparentctrl = (SumoPhysicsController* )parentctrl;
        
        if (sumoparentctrl)
        {
                dynaparent = sumoparentctrl->GetSumoObject();
        }
        
        SM_Object* orgsumoobject = m_sumoObj;
        
        
        m_sumoObj       =       new SM_Object(
                orgsumoobject->getShapeHandle(), 
                orgsumoobject->getMaterialProps(),                      
                orgsumoobject->getShapeProps(),
                dynaparent);
        
        m_sumoObj->setRigidBody(orgsumoobject->isRigidBody());
        
        m_sumoObj->setMargin(orgsumoobject->getMargin());
        m_sumoObj->setPosition(orgsumoobject->getPosition());
        m_sumoObj->setOrientation(orgsumoobject->getOrientation());
        //if it is a dyna, register for a callback
        m_sumoObj->registerCallback(*this);
        
        m_sumoScene->add(* (m_sumoObj));
        */



}

void    CcdPhysicsController::SetPhysicsEnvironment(class PHY_IPhysicsEnvironment *env)
{
        // can safely assume CCD environment
        CcdPhysicsEnvironment *physicsEnv = static_cast<CcdPhysicsEnvironment*>(env);

        if (m_cci.m_physicsEnv != physicsEnv) 
        {
                // since the environment is changing, we must also move the controler to the
                // new environement. Note that we don't handle sensor explicitely: this
                // function can be called on sensor but only when they are not registered
                if (m_cci.m_physicsEnv->removeCcdPhysicsController(this))
                {
                        physicsEnv->addCcdPhysicsController(this);
                }
                m_cci.m_physicsEnv = physicsEnv;
        }
}

void    CcdPhysicsController::SetCenterOfMassTransform(btTransform& xform)
{
        btRigidBody* body = GetRigidBody();
        if (body)
        {
                body->setCenterOfMassTransform(xform);
        } else
        {
                //either collision object or soft body?
                if (GetSoftBody())
                {

                } else
                {

                        if (m_object->isStaticOrKinematicObject())
                        {
                                m_object->setInterpolationWorldTransform(m_object->getWorldTransform());
                        } else
                        {
                                m_object->setInterpolationWorldTransform(xform);
                        }
                        if (body)
                        {
                                body->setInterpolationLinearVelocity(body->getLinearVelocity());
                                body->setInterpolationAngularVelocity(body->getAngularVelocity());
                                body->updateInertiaTensor();
                        }
                        m_object->setWorldTransform(xform);
                }
        }
}

                // kinematic methods
void            CcdPhysicsController::RelativeTranslate(float dlocX,float dlocY,float dlocZ,bool local)
{
        if (m_object)
        {
                m_object->activate(true);
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        // kinematic object should not set the transform, it disturbs the velocity interpolation
                        return;
                }

                // btRigidBody* body = GetRigidBody(); // not used anymore

                btVector3 dloc(dlocX,dlocY,dlocZ);
                btTransform xform = m_object->getWorldTransform();
        
                if (local)
                {
                        dloc = xform.getBasis()*dloc;
                }

                xform.setOrigin(xform.getOrigin() + dloc);
                SetCenterOfMassTransform(xform);
        }

}

void            CcdPhysicsController::RelativeRotate(const float rotval[9],bool local)
{
        if (m_object)
        {
                m_object->activate(true);
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        // kinematic object should not set the transform, it disturbs the velocity interpolation
                        return;
                }

                btMatrix3x3 drotmat(    rotval[0],rotval[4],rotval[8],
                                                                rotval[1],rotval[5],rotval[9],
                                                                rotval[2],rotval[6],rotval[10]);


                btMatrix3x3 currentOrn;
                GetWorldOrientation(currentOrn);

                btTransform xform = m_object->getWorldTransform();
                
                xform.setBasis(xform.getBasis()*(local ? 
                drotmat : (currentOrn.inverse() * drotmat * currentOrn)));

                SetCenterOfMassTransform(xform);
        }
}


void CcdPhysicsController::GetWorldOrientation(btMatrix3x3& mat)
{
        float ori[12];
        m_MotionState->getWorldOrientation(ori);
        mat.setFromOpenGLSubMatrix(ori);
}

void            CcdPhysicsController::getOrientation(float &quatImag0,float &quatImag1,float &quatImag2,float &quatReal)
{
        btQuaternion q = m_object->getWorldTransform().getRotation();
        quatImag0 = q[0];
        quatImag1 = q[1];
        quatImag2 = q[2];
        quatReal = q[3];
}
void            CcdPhysicsController::setOrientation(float quatImag0,float quatImag1,float quatImag2,float quatReal)
{
        if (m_object)
        {
                m_object->activate(true);
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        // kinematic object should not set the transform, it disturbs the velocity interpolation
                        return;
                }
                // not required
                //m_MotionState->setWorldOrientation(quatImag0,quatImag1,quatImag2,quatReal);
                btTransform xform  = m_object->getWorldTransform();
                xform.setRotation(btQuaternion(quatImag0,quatImag1,quatImag2,quatReal));
                SetCenterOfMassTransform(xform);
                // not required
                //m_bulletMotionState->setWorldTransform(xform);
                
                

        }

}

void CcdPhysicsController::setWorldOrientation(const btMatrix3x3& orn)
{
        if (m_object)
        {
                m_object->activate(true);
                if (m_object->isStaticObject() && !m_cci.m_bSensor)
                {
                        m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                }
                // not required
                //m_MotionState->setWorldOrientation(quatImag0,quatImag1,quatImag2,quatReal);
                btTransform xform  = m_object->getWorldTransform();
                xform.setBasis(orn);
                SetCenterOfMassTransform(xform);
                // not required
                //m_bulletMotionState->setWorldTransform(xform);
                //only once!
                if (!m_softBodyTransformInitialized && GetSoftBody())
                {
                        m_softbodyStartTrans.setBasis(orn);
                        xform.setOrigin(m_softbodyStartTrans.getOrigin());
                        GetSoftBody()->transform(xform);
                        m_softBodyTransformInitialized = true;
                }

        }

}

void            CcdPhysicsController::setPosition(float posX,float posY,float posZ)
{
        if (m_object)
        {
                m_object->activate(true);
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        // kinematic object should not set the transform, it disturbs the velocity interpolation
                        return;
                }
                // not required, this function is only used to update the physic controller
                //m_MotionState->setWorldPosition(posX,posY,posZ);
                btTransform xform  = m_object->getWorldTransform();
                xform.setOrigin(btVector3(posX,posY,posZ));
                SetCenterOfMassTransform(xform);
                if (!m_softBodyTransformInitialized)
                        m_softbodyStartTrans.setOrigin(xform.getOrigin());
                // not required
                //m_bulletMotionState->setWorldTransform(xform);
        }
}

void CcdPhysicsController::forceWorldTransform(const btMatrix3x3& mat, const btVector3& pos)
{
        if (m_object)
        {
                btTransform& xform = m_object->getWorldTransform();
                xform.setBasis(mat);
                xform.setOrigin(pos);
        }
}


void            CcdPhysicsController::resolveCombinedVelocities(float linvelX,float linvelY,float linvelZ,float angVelX,float angVelY,float angVelZ)
{
}

void            CcdPhysicsController::getPosition(PHY__Vector3& pos) const
{
        const btTransform& xform = m_object->getWorldTransform();
        pos[0] = xform.getOrigin().x();
        pos[1] = xform.getOrigin().y();
        pos[2] = xform.getOrigin().z();
}

void            CcdPhysicsController::setScaling(float scaleX,float scaleY,float scaleZ)
{
        if (!btFuzzyZero(m_cci.m_scaling.x()-scaleX) ||
                !btFuzzyZero(m_cci.m_scaling.y()-scaleY) ||
                !btFuzzyZero(m_cci.m_scaling.z()-scaleZ))
        {
                m_cci.m_scaling = btVector3(scaleX,scaleY,scaleZ);

                if (m_object && m_object->getCollisionShape())
                {
                        m_object->activate(true); // without this, sleeping objects scale wont be applied in bullet if python changes the scale - Campbell.
                        m_object->getCollisionShape()->setLocalScaling(m_cci.m_scaling);
                        
                        //printf("no inertia recalc for fixed objects with mass=0\n");
                        btRigidBody* body = GetRigidBody();
                        if (body && m_cci.m_mass)
                        {
                                body->getCollisionShape()->calculateLocalInertia(m_cci.m_mass, m_cci.m_localInertiaTensor);
                                body->setMassProps(m_cci.m_mass, m_cci.m_localInertiaTensor * m_cci.m_inertiaFactor);
                        } 
                        
                }
        }
}
                
                // physics methods
void            CcdPhysicsController::ApplyTorque(float torqueX,float torqueY,float torqueZ,bool local)
{
        btVector3 torque(torqueX,torqueY,torqueZ);
        btTransform xform = m_object->getWorldTransform();
        

        if (m_object && torque.length2() > (SIMD_EPSILON*SIMD_EPSILON))
        {
                btRigidBody* body = GetRigidBody();
                m_object->activate();
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        return;
                }
                if (local)
                {
                        torque  = xform.getBasis()*torque;
                }
                if (body)
                {
                        if      (m_cci.m_bRigid)
                        {
                                body->applyTorque(torque);
                        }
                        else
                        {
                                //workaround for incompatibility between 'DYNAMIC' game object, and angular factor
                                //a DYNAMIC object has some inconsistency: it has no angular effect due to collisions, but still has torque
                                const btVector3& angFac = body->getAngularFactor();
                                btVector3 tmpFac(1,1,1);
                                body->setAngularFactor(tmpFac);
                                body->applyTorque(torque);
                                body->setAngularFactor(angFac);
                        } 
                } 
        }
}

void            CcdPhysicsController::ApplyForce(float forceX,float forceY,float forceZ,bool local)
{
        btVector3 force(forceX,forceY,forceZ);
        

        if (m_object && force.length2() > (SIMD_EPSILON*SIMD_EPSILON))
        {
                m_object->activate();
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        return;
                }
                btTransform xform = m_object->getWorldTransform();
                
                if (local)
                {       
                        force   = xform.getBasis()*force;
                }
                btRigidBody* body = GetRigidBody();
                if (body)
                        body->applyCentralForce(force);
                btSoftBody* soft = GetSoftBody();
                if (soft)
                {
                        // the force is applied on each node, must reduce it in the same extend
                        if (soft->m_nodes.size() > 0)
                                force /= soft->m_nodes.size();
                        soft->addForce(force);
                }
        }
}
void            CcdPhysicsController::SetAngularVelocity(float ang_velX,float ang_velY,float ang_velZ,bool local)
{
        btVector3 angvel(ang_velX,ang_velY,ang_velZ);
        if (m_object && angvel.length2() > (SIMD_EPSILON*SIMD_EPSILON))
        {
                m_object->activate(true);
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        return;
                }
                btTransform xform = m_object->getWorldTransform();
                if (local)
                {
                        angvel  = xform.getBasis()*angvel;
                }
                btRigidBody* body = GetRigidBody();
                if (body)
                        body->setAngularVelocity(angvel);
        }

}
void            CcdPhysicsController::SetLinearVelocity(float lin_velX,float lin_velY,float lin_velZ,bool local)
{

        btVector3 linVel(lin_velX,lin_velY,lin_velZ);
        if (m_object/* && linVel.length2() > (SIMD_EPSILON*SIMD_EPSILON)*/)
        {
                m_object->activate(true);
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        return;
                }
                
                btSoftBody* soft = GetSoftBody();
                if (soft)
                {
                        if (local)
                        {
                                linVel  = m_softbodyStartTrans.getBasis()*linVel;
                        }
                        soft->setVelocity(linVel);
                } else
                {
                        btTransform xform = m_object->getWorldTransform();
                        if (local)
                        {
                                linVel  = xform.getBasis()*linVel;
                        }
                        btRigidBody* body = GetRigidBody();
                        if (body)
                                body->setLinearVelocity(linVel);
                }
        }
}
void            CcdPhysicsController::applyImpulse(float attachX,float attachY,float attachZ, float impulseX,float impulseY,float impulseZ)
{
        btVector3 impulse(impulseX,impulseY,impulseZ);

        if (m_object && impulse.length2() > (SIMD_EPSILON*SIMD_EPSILON))
        {
                m_object->activate();
                if (m_object->isStaticObject())
                {
                        if (!m_cci.m_bSensor)
                                m_object->setCollisionFlags(m_object->getCollisionFlags() | btCollisionObject::CF_KINEMATIC_OBJECT);
                        return;
                }
                
                btVector3 pos(attachX,attachY,attachZ);
                btRigidBody* body = GetRigidBody();
                if (body)
                        body->applyImpulse(impulse,pos);
                        
        }

}
void            CcdPhysicsController::SetActive(bool active)
{
}
                // reading out information from physics
void            CcdPhysicsController::GetLinearVelocity(float& linvX,float& linvY,float& linvZ)
{
        btRigidBody* body = GetRigidBody();
        if (body)
        {
                const btVector3& linvel = body->getLinearVelocity();
                linvX = linvel.x();
                linvY = linvel.y();
                linvZ = linvel.z();
        } else
        {
                linvX = 0.f;
                linvY = 0.f;
                linvZ = 0.f;
        }

}

void            CcdPhysicsController::GetAngularVelocity(float& angVelX,float& angVelY,float& angVelZ)
{
        btRigidBody* body = GetRigidBody();
        if (body)
        {
                const btVector3& angvel= body->getAngularVelocity();
                angVelX = angvel.x();
                angVelY = angvel.y();
                angVelZ = angvel.z();
        } else
        {
                angVelX = 0.f;
                angVelY = 0.f;
                angVelZ = 0.f;
        }
}

void            CcdPhysicsController::GetVelocity(const float posX,const float posY,const float posZ,float& linvX,float& linvY,float& linvZ)
{
        btVector3 pos(posX,posY,posZ);
        btRigidBody* body = GetRigidBody();
        if (body)
        {
                btVector3 linvel = body->getVelocityInLocalPoint(pos);
                linvX = linvel.x();
                linvY = linvel.y();
                linvZ = linvel.z();
        } else
        {
                linvX = 0.f;
                linvY = 0.f;
                linvZ = 0.f;
        }
}
void            CcdPhysicsController::getReactionForce(float& forceX,float& forceY,float& forceZ)
{
}

                // dyna's that are rigidbody are free in orientation, dyna's with non-rigidbody are restricted 
void            CcdPhysicsController::setRigidBody(bool rigid)
{
        if (!rigid)
        {
                btRigidBody* body = GetRigidBody();
                if (body)
                {
                        //fake it for now
                        btVector3 inertia = body->getInvInertiaDiagLocal();
                        inertia[1] = 0.f;
                        body->setInvInertiaDiagLocal(inertia);
                        body->updateInertiaTensor();
                }
        }
}

                // clientinfo for raycasts for example
void*           CcdPhysicsController::getNewClientInfo()
{
        return m_newClientInfo;
}
void            CcdPhysicsController::setNewClientInfo(void* clientinfo)
{
        m_newClientInfo = clientinfo;
}


void    CcdPhysicsController::UpdateDeactivation(float timeStep)
{
        btRigidBody* body = GetRigidBody();
        if (body)
        {
                body->updateDeactivation( timeStep);
        }
}

bool CcdPhysicsController::wantsSleeping()
{
        btRigidBody* body = GetRigidBody();
        if (body)
        {
                return body->wantsSleeping();
        }
        //check it out
        return true;
}

PHY_IPhysicsController* CcdPhysicsController::GetReplica()
{
        // This is used only to replicate Near and Radar sensor controllers
        // The replication of object physics controller is done in KX_BulletPhysicsController::GetReplica()
        CcdConstructionInfo cinfo = m_cci;
        if (m_shapeInfo)
        {
                // This situation does not normally happen
                cinfo.m_collisionShape = m_shapeInfo->CreateBulletShape(m_cci.m_margin, m_cci.m_bGimpact, !m_cci.m_bSoft);
        } 
        else if (m_collisionShape)
        {
                switch (m_collisionShape->getShapeType())
                {
                case SPHERE_SHAPE_PROXYTYPE:
                        {
                                btSphereShape* orgShape = (btSphereShape*)m_collisionShape;
                                cinfo.m_collisionShape = new btSphereShape(*orgShape);
                                break;
                        }

                case CONE_SHAPE_PROXYTYPE:
                        {
                                btConeShape* orgShape = (btConeShape*)m_collisionShape;
                                cinfo.m_collisionShape = new btConeShape(*orgShape);
                                break;
                        }

                default:
                        {
                                return 0;
                        }
                }
        }

        cinfo.m_MotionState = new DefaultMotionState();
        cinfo.m_shapeInfo = m_shapeInfo;

        CcdPhysicsController* replica = new CcdPhysicsController(cinfo);
        return replica;
}


DefaultMotionState::DefaultMotionState()
{
        m_worldTransform.setIdentity();
        m_localScaling.setValue(1.f,1.f,1.f);
}


DefaultMotionState::~DefaultMotionState()
{

}

void    DefaultMotionState::getWorldPosition(float& posX,float& posY,float& posZ)
{
        posX = m_worldTransform.getOrigin().x();
        posY = m_worldTransform.getOrigin().y();
        posZ = m_worldTransform.getOrigin().z();
}

void    DefaultMotionState::getWorldScaling(float& scaleX,float& scaleY,float& scaleZ)
{
        scaleX = m_localScaling.getX();
        scaleY = m_localScaling.getY();
        scaleZ = m_localScaling.getZ();
}

void    DefaultMotionState::getWorldOrientation(float& quatIma0,float& quatIma1,float& quatIma2,float& quatReal)
{
        btQuaternion quat = m_worldTransform.getRotation();
        quatIma0 = quat.x();
        quatIma1 = quat.y();
        quatIma2 = quat.z();
        quatReal = quat[3];
}
                
void    DefaultMotionState::getWorldOrientation(float* ori)
{
        m_worldTransform.getBasis().getOpenGLSubMatrix(ori);
}

void    DefaultMotionState::setWorldOrientation(const float* ori)
{
        m_worldTransform.getBasis().setFromOpenGLSubMatrix(ori);
}
void    DefaultMotionState::setWorldPosition(float posX,float posY,float posZ)
{
        btVector3 pos(posX,posY,posZ);
        m_worldTransform.setOrigin( pos );
}

void    DefaultMotionState::setWorldOrientation(float quatIma0,float quatIma1,float quatIma2,float quatReal)
{
        btQuaternion orn(quatIma0,quatIma1,quatIma2,quatReal);
        m_worldTransform.setRotation( orn );
}
                
void    DefaultMotionState::calculateWorldTransformations()
{

}

// Shape constructor
std::map<RAS_MeshObject*, CcdShapeConstructionInfo*> CcdShapeConstructionInfo::m_meshShapeMap;

CcdShapeConstructionInfo* CcdShapeConstructionInfo::FindMesh(RAS_MeshObject* mesh, struct DerivedMesh* dm, bool polytope)
{
        if (polytope || dm)
                // not yet supported
                return NULL;

        std::map<RAS_MeshObject*,CcdShapeConstructionInfo*>::const_iterator mit = m_meshShapeMap.find(mesh);
        if (mit != m_meshShapeMap.end())
                return mit->second;
        return NULL;
}

bool CcdShapeConstructionInfo::SetMesh(RAS_MeshObject* meshobj, DerivedMesh* dm, bool polytope)
{
        int numpolys, numverts;

        // assume no shape information
        // no support for dynamic change of shape yet
        assert(IsUnused());
        m_shapeType = PHY_SHAPE_NONE;
        m_meshObject = NULL;
        bool free_dm = false;

        // No mesh object or mesh has no polys
        if (!meshobj || meshobj->HasColliderPolygon()==false) {
                m_vertexArray.clear();
                m_polygonIndexArray.clear();
                m_triFaceArray.clear();
                m_triFaceUVcoArray.clear();
                return false;
        }

        if (!dm) {
                free_dm = true;
                dm = CDDM_from_mesh(meshobj->GetMesh(), NULL);
        }

        MVert *mvert = dm->getVertArray(dm);
        MFace *mface = dm->getFaceArray(dm);
        numpolys = dm->getNumFaces(dm);
        numverts = dm->getNumVerts(dm);
        int* index = (int*)dm->getFaceDataArray(dm, CD_ORIGINDEX);
        MTFace *tface = (MTFace *)dm->getFaceDataArray(dm, CD_MTFACE);

        m_shapeType = (polytope) ? PHY_SHAPE_POLYTOPE : PHY_SHAPE_MESH;

        /* Convert blender geometry into bullet mesh, need these vars for mapping */
        vector<bool> vert_tag_array(numverts, false);
        unsigned int tot_bt_verts= 0;

        if (polytope)
        {
                // Tag verts we're using
                for (int p2=0; p2<numpolys; p2++)
                {
                        MFace* mf = &mface[p2];
                        RAS_Polygon* poly = meshobj->GetPolygon((index)? index[p2]: p2);

                        // only add polygons that have the collision flag set
                        if (poly->IsCollider())
                        {
                                if (vert_tag_array[mf->v1]==false) {vert_tag_array[mf->v1]= true;tot_bt_verts++;}
                                if (vert_tag_array[mf->v2]==false) {vert_tag_array[mf->v2]= true;tot_bt_verts++;}
                                if (vert_tag_array[mf->v3]==false) {vert_tag_array[mf->v3]= true;tot_bt_verts++;}
                                if (mf->v4 && vert_tag_array[mf->v4]==false) {vert_tag_array[mf->v4]= true;tot_bt_verts++;}
                        }
                }

                m_vertexArray.resize(tot_bt_verts*3);

                btScalar *bt= &m_vertexArray[0];

                for (int p2=0; p2<numpolys; p2++)
                {
                        MFace* mf = &mface[p2];
                        RAS_Polygon* poly= meshobj->GetPolygon((index)? index[p2]: p2);

                        // only add polygons that have the collisionflag set
                        if (poly->IsCollider())
                        {
                                if (vert_tag_array[mf->v1]==true)
                                {
                                        const float* vtx = mvert[mf->v1].co;
                                        vert_tag_array[mf->v1]= false;
                                        *bt++ = vtx[0];
                                        *bt++ = vtx[1];
                                        *bt++ = vtx[2];
                                }
                                if (vert_tag_array[mf->v2]==true)
                                {
                                        const float* vtx = mvert[mf->v2].co;
                                        vert_tag_array[mf->v2]= false;
                                        *bt++ = vtx[0];
                                        *bt++ = vtx[1];
                                        *bt++ = vtx[2];
                                }
                                if (vert_tag_array[mf->v3]==true)
                                {
                                        const float* vtx = mvert[mf->v3].co;
                                        vert_tag_array[mf->v3]= false;
                                        *bt++ = vtx[0];
                                        *bt++ = vtx[1];
                                        *bt++ = vtx[2];
                                }
                                if (mf->v4 && vert_tag_array[mf->v4]==true)
                                {
                                        const float* vtx = mvert[mf->v4].co;
                                        vert_tag_array[mf->v4]= false;
                                        *bt++ = vtx[0];
                                        *bt++ = vtx[1];
                                        *bt++ = vtx[2];
                                }
                        }
                }
        }
        else {
                unsigned int tot_bt_tris= 0;
                vector<int> vert_remap_array(numverts, 0);
                
                // Tag verts we're using
                for (int p2=0; p2<numpolys; p2++)
                {
                        MFace* mf = &mface[p2];
                        RAS_Polygon* poly= meshobj->GetPolygon((index)? index[p2]: p2);

                        // only add polygons that have the collision flag set
                        if (poly->IsCollider())
                        {
                                if (vert_tag_array[mf->v1]==false)
                                        {vert_tag_array[mf->v1]= true;vert_remap_array[mf->v1]= tot_bt_verts;tot_bt_verts++;}
                                if (vert_tag_array[mf->v2]==false)
                                        {vert_tag_array[mf->v2]= true;vert_remap_array[mf->v2]= tot_bt_verts;tot_bt_verts++;}
                                if (vert_tag_array[mf->v3]==false)
                                        {vert_tag_array[mf->v3]= true;vert_remap_array[mf->v3]= tot_bt_verts;tot_bt_verts++;}
                                if (mf->v4 && vert_tag_array[mf->v4]==false)
                                        {vert_tag_array[mf->v4]= true;vert_remap_array[mf->v4]= tot_bt_verts;tot_bt_verts++;}
                                tot_bt_tris += (mf->v4 ? 2:1); /* a quad or a tri */
                        }
                }

                m_vertexArray.resize(tot_bt_verts*3);
                m_polygonIndexArray.resize(tot_bt_tris);
                m_triFaceArray.resize(tot_bt_tris*3);
                btScalar *bt= &m_vertexArray[0];
                int *poly_index_pt= &m_polygonIndexArray[0];
                int *tri_pt= &m_triFaceArray[0];

                UVco *uv_pt = NULL;
                if (tface)
                {
                        m_triFaceUVcoArray.resize(tot_bt_tris*3);
                        uv_pt = &m_triFaceUVcoArray[0];
                } 
                else 
                        m_triFaceUVcoArray.clear();

                for (int p2=0; p2<numpolys; p2++)
                {
                        MFace* mf = &mface[p2];
                        MTFace* tf = (tface) ? &tface[p2] : NULL;
                        RAS_Polygon* poly= meshobj->GetPolygon((index)? index[p2]: p2);

                        // only add polygons that have the collisionflag set
                        if (poly->IsCollider())
                        {
                                MVert *v1= &mvert[mf->v1];
                                MVert *v2= &mvert[mf->v2];
                                MVert *v3= &mvert[mf->v3];

                                // the face indicies
                                tri_pt[0]= vert_remap_array[mf->v1];
                                tri_pt[1]= vert_remap_array[mf->v2];
                                tri_pt[2]= vert_remap_array[mf->v3];
                                tri_pt= tri_pt+3;
                                if (tf)
                                {
                                        uv_pt[0].uv[0] = tf->uv[0][0];
                                        uv_pt[0].uv[1] = tf->uv[0][1];
                                        uv_pt[1].uv[0] = tf->uv[1][0];
                                        uv_pt[1].uv[1] = tf->uv[1][1];
                                        uv_pt[2].uv[0] = tf->uv[2][0];
                                        uv_pt[2].uv[1] = tf->uv[2][1];
                                        uv_pt += 3;
                                }

                                // m_polygonIndexArray
                                *poly_index_pt= (index)? index[p2]: p2;
                                poly_index_pt++;

                                // the vertex location
                                if (vert_tag_array[mf->v1]==true) { /* *** v1 *** */
                                        vert_tag_array[mf->v1]= false;
                                        *bt++ = v1->co[0];
                                        *bt++ = v1->co[1];
                                        *bt++ = v1->co[2];
                                }
                                if (vert_tag_array[mf->v2]==true) { /* *** v2 *** */
                                        vert_tag_array[mf->v2]= false;
                                        *bt++ = v2->co[0];
                                        *bt++ = v2->co[1];
                                        *bt++ = v2->co[2];
                                }
                                if (vert_tag_array[mf->v3]==true) { /* *** v3 *** */
                                        vert_tag_array[mf->v3]= false;
                                        *bt++ = v3->co[0];      
                                        *bt++ = v3->co[1];
                                        *bt++ = v3->co[2];
                                }

                                if (mf->v4)
                                {
                                        MVert *v4= &mvert[mf->v4];

                                        tri_pt[0]= vert_remap_array[mf->v1];
                                        tri_pt[1]= vert_remap_array[mf->v3];
                                        tri_pt[2]= vert_remap_array[mf->v4];
                                        tri_pt= tri_pt+3;
                                        if (tf)
                                        {
                                                uv_pt[0].uv[0] = tf->uv[0][0];
                                                uv_pt[0].uv[1] = tf->uv[0][1];
                                                uv_pt[1].uv[0] = tf->uv[2][0];
                                                uv_pt[1].uv[1] = tf->uv[2][1];
                                                uv_pt[2].uv[0] = tf->uv[3][0];
                                                uv_pt[2].uv[1] = tf->uv[3][1];
                                                uv_pt += 3;
                                        }

                                        // m_polygonIndexArray
                                        *poly_index_pt= (index)? index[p2]: p2;
                                        poly_index_pt++;

                                        // the vertex location
                                        if (vert_tag_array[mf->v4]==true) { /* *** v4 *** */
                                                vert_tag_array[mf->v4]= false;
                                                *bt++ = v4->co[0];
                                                *bt++ = v4->co[1];      
                                                *bt++ = v4->co[2];
                                        }
                                }
                        }
                }


                /* If this ever gets confusing, print out an OBJ file for debugging */
#if 0
                printf("# vert count %d\n", m_vertexArray.size());
                for(i=0; i<m_vertexArray.size(); i+=1) {
                        printf("v %.6f %.6f %.6f\n", m_vertexArray[i].x(), m_vertexArray[i].y(), m_vertexArray[i].z());
                }

                printf("# face count %d\n", m_triFaceArray.size());
                for(i=0; i<m_triFaceArray.size(); i+=3) {
                        printf("f %d %d %d\n", m_triFaceArray[i]+1, m_triFaceArray[i+1]+1, m_triFaceArray[i+2]+1);
                }
#endif

        }

#if 0
        if (validpolys==false)
        {
                // should not happen
                m_shapeType = PHY_SHAPE_NONE;
                return false;
        }
#endif
        
        m_meshObject = meshobj;
        if (free_dm) {
                dm->release(dm);
                dm = NULL;
        }

        // sharing only on static mesh at present, if you change that, you must also change in FindMesh
        if (!polytope && !dm)
        {
                // triangle shape can be shared, store the mesh object in the map
                m_meshShapeMap.insert(std::pair<RAS_MeshObject*,CcdShapeConstructionInfo*>(meshobj,this));
        }
        return true;
}

#include <cstdio>

/* Updates the arrays used by CreateBulletShape(),
 * take care that recalcLocalAabb() runs after CreateBulletShape is called.
 * */
bool CcdShapeConstructionInfo::UpdateMesh(class KX_GameObject* gameobj, class RAS_MeshObject* meshobj)
{
        int numpolys;
        int numverts;

        unsigned int tot_bt_tris= 0;
        unsigned int tot_bt_verts= 0;

        int i, j;
        int v_orig;

        /* Use for looping over verts in a face as a try or 2 tris */
        const int quad_verts[7]=        {0,1,2,          0,2,3,         -1};
        const int tri_verts[4]= {0,1,2,         -1};
        const int *fv_pt;

        if(gameobj==NULL && meshobj==NULL)
                return false;
        
        if(m_shapeType != PHY_SHAPE_MESH)
                return false;

        RAS_Deformer *deformer= gameobj ? gameobj->GetDeformer():NULL;
        DerivedMesh* dm = NULL;

        if (deformer)
                dm = deformer->GetPhysicsMesh();
        
        /* get the mesh from the object if not defined */
        if(meshobj==NULL) {
                
                /* modifier mesh */
                if(dm)
                        meshobj= deformer->GetRasMesh();
                
                /* game object first mesh */
                if(meshobj==NULL) {
                        if(gameobj->GetMeshCount() > 0) {
                                meshobj= gameobj->GetMesh(0);
                        }
                }
        }
        
        if(dm && deformer->GetRasMesh() == meshobj)
        {       /*
                 * Derived Mesh Update
                 *
                 * */

                MVert *mvert = dm->getVertArray(dm);
                MFace *mface = dm->getFaceArray(dm);
                numpolys = dm->getNumFaces(dm);
                numverts = dm->getNumVerts(dm);
                int* index = (int*)dm->getFaceDataArray(dm, CD_ORIGINDEX);

                MFace *mf;
                MVert *mv;

                int flen;

                if(CustomData_has_layer(&dm->faceData, CD_MTFACE))
                {
                        MTFace *tface = (MTFace *)dm->getFaceDataArray(dm, CD_MTFACE);
                        MTFace *tf;

                        vector<bool> vert_tag_array(numverts, false);
                        vector<int> vert_remap_array(numverts, 0);

                        for(mf= mface, tf= tface, i=0; i < numpolys; mf++, tf++, i++) {
                                if(tf->mode & TF_DYNAMIC)
                                {
                                        if(mf->v4) {
                                                tot_bt_tris+= 2;
                                                flen= 4;
                                        } else {
                                                tot_bt_tris++;
                                                flen= 3;
                                        }

                                        for(j=0; j<flen; j++)
                                        {
                                                v_orig = (*(&mf->v1 + j));

                                                if(vert_tag_array[v_orig]==false)
                                                {
                                                        vert_tag_array[v_orig]= true;
                                                        vert_remap_array[v_orig]= tot_bt_verts;
                                                        tot_bt_verts++;
                                                }
                                        }
                                }
                        }

                        m_vertexArray.resize(tot_bt_verts*3);
                        btScalar *bt= &m_vertexArray[0];

                        m_triFaceArray.resize(tot_bt_tris*3);
                        int *tri_pt= &m_triFaceArray[0];

                        m_triFaceUVcoArray.resize(tot_bt_tris*3);
                        UVco *uv_pt= &m_triFaceUVcoArray[0];

                        m_polygonIndexArray.resize(tot_bt_tris);
                        int *poly_index_pt= &m_polygonIndexArray[0];

                        for(mf= mface, tf= tface, i=0; i < numpolys; mf++, tf++, i++)
                        {
                                if(tf->mode & TF_DYNAMIC)
                                {
                                        int origi = (index)? index[i]: i;

                                        if(mf->v4) {
                                                fv_pt= quad_verts;
                                                *poly_index_pt++ = origi;
                                                *poly_index_pt++ = origi;
                                                flen= 4;
                                        } else {
                                                fv_pt= tri_verts;
                                                *poly_index_pt++ = origi;
                                                flen= 3;
                                        }

                                        for(; *fv_pt > -1; fv_pt++)
                                        {
                                                v_orig = (*(&mf->v1 + (*fv_pt)));

                                                if(vert_tag_array[v_orig])
                                                {
                                                        mv= mvert + v_orig;
                                                        *bt++ = mv->co[0];
                                                        *bt++ = mv->co[1];
                                                        *bt++ = mv->co[2];

                                                        vert_tag_array[v_orig]= false;
                                                }
                                                *tri_pt++ = vert_remap_array[v_orig];
                                                uv_pt->uv[0] = tf->uv[*fv_pt][0];
                                                uv_pt->uv[1] = tf->uv[*fv_pt][1];
                                                uv_pt++;
                                        }
                                }
                        }
                }
                else {
                        /* no need for a vertex mapping. simple/fast */

                        tot_bt_verts= numverts;

                        for(mf= mface, i=0; i < numpolys; mf++, i++) {
                                tot_bt_tris += (mf->v4 ? 2:1);
                        }

                        m_vertexArray.resize(tot_bt_verts*3);
                        btScalar *bt= &m_vertexArray[0];

                        m_triFaceArray.resize(tot_bt_tris*3);
                        int *tri_pt= &m_triFaceArray[0];

                        m_polygonIndexArray.resize(tot_bt_tris);
                        int *poly_index_pt= &m_polygonIndexArray[0];

                        m_triFaceUVcoArray.clear();

                        for(mv= mvert, i=0; i < numverts; mv++, i++) {
                                *bt++ = mv->co[0]; *bt++ = mv->co[1]; *bt++ = mv->co[2];
                        }

                        for(mf= mface, i=0; i < numpolys; mf++, i++) {
                                int origi = (index)? index[i]: i;

                                if(mf->v4) {
                                        fv_pt= quad_verts;
                                        *poly_index_pt++ = origi;
                                        *poly_index_pt++ = origi;
                                }
                                else {
                                        fv_pt= tri_verts;
                                        *poly_index_pt++ = origi;
                                }

                                for(; *fv_pt > -1; fv_pt++)
                                        *tri_pt++ = (*(&mf->v1 + (*fv_pt)));
                        }
                }
        }
        else {  /*
                         * RAS Mesh Update
                         *
                         * */
                
                /* Note!, gameobj can be NULL here */

                /* transverts are only used for deformed RAS_Meshes, the RAS_TexVert data
                 * is too hard to get at, see below for details */
                float (*transverts)[3]= NULL;
                int transverts_tot= 0; /* with deformed meshes - should always be greater then the max orginal index, or we get crashes */

                if(deformer) {
                        /* map locations from the deformed array
                         *
                         * Could call deformer->Update(); but rely on redraw updating.
                         * */
                        transverts= deformer->GetTransVerts(&transverts_tot);
                }

                // Tag verts we're using
                numpolys= meshobj->NumPolygons();
                numverts= meshobj->m_sharedvertex_map.size();
                const float *xyz;


                vector<bool> vert_tag_array(numverts, false);
                vector<int> vert_remap_array(numverts, 0);

                for(int p=0; p<numpolys; p++)
                {
                        RAS_Polygon* poly= meshobj->GetPolygon(p);
                        if (poly->IsCollider())
                        {
                                for(i=0; i < poly->VertexCount(); i++)
                                {
                                        v_orig= poly->GetVertex(i)->getOrigIndex();
                                        if(vert_tag_array[v_orig]==false)
                                        {
                                                vert_tag_array[v_orig]= true;
                                                vert_remap_array[v_orig]= tot_bt_verts;
                                                tot_bt_verts++;
                                        }
                                }
                                tot_bt_tris += (poly->VertexCount()==4 ? 2:1);
                        }
                }

                m_vertexArray.resize(tot_bt_verts*3);
                btScalar *bt= &m_vertexArray[0];

                m_triFaceArray.resize(tot_bt_tris*3);
                int *tri_pt= &m_triFaceArray[0];

                /* cant be used for anything useful in this case, since we dont rely on the original mesh
                 * will just be an array like pythons range(tot_bt_tris) */
                m_polygonIndexArray.resize(tot_bt_tris);


                for(int p=0; p<numpolys; p++)
                {
                        RAS_Polygon* poly= meshobj->GetPolygon(p);

                        if (poly->IsCollider())
                        {
                                /* quad or tri loop */
                                fv_pt= (poly->VertexCount()==3 ? tri_verts:quad_verts);

                                for(; *fv_pt > -1; fv_pt++)
                                {
                                        v_orig= poly->GetVertex(*fv_pt)->getOrigIndex();

                                        if(vert_tag_array[v_orig])
                                        {
                                                if(transverts) {
                                                        /* deformed mesh, using RAS_TexVert locations would be too troublesome
                                                         * because they are use the gameob as a hash in the material slot */
                                                        *bt++ = transverts[v_orig][0];
                                                        *bt++ = transverts[v_orig][1];
                                                        *bt++ = transverts[v_orig][2];
                                                }
                                                else {
                                                        /* static mesh python may have modified */
                                                        xyz= meshobj->GetVertexLocation( v_orig );
                                                        *bt++ = xyz[0];
                                                        *bt++ = xyz[1];
                                                        *bt++ = xyz[2];
                                                }

                                                vert_tag_array[v_orig]= false;
                                        }

                                        *tri_pt++ = vert_remap_array[v_orig];
                                }
                        }

                        m_polygonIndexArray[p]= p; /* dumb counting */
                }
        }
        
#if 0
        /* needs #include <cstdio> */
        printf("# vert count %d\n", m_vertexArray.size());
        for(int i=0; i<m_vertexArray.size(); i+=3) {
                printf("v %.6f %.6f %.6f\n", m_vertexArray[i], m_vertexArray[i+1], m_vertexArray[i+2]);
        }

        printf("# face count %d\n", m_triFaceArray.size());
        for(int i=0; i<m_triFaceArray.size(); i+=3) {
                printf("f %d %d %d\n", m_triFaceArray[i]+1, m_triFaceArray[i+1]+1, m_triFaceArray[i+2]+1);
        }
#endif

        /* force recreation of the m_unscaledShape.
         * If this has multiple users we cant delete */
        if(m_unscaledShape) {
                // dont free now so it can re-allocate under the same location and not break pointers.
                // DeleteBulletShape(m_unscaledShape); 
                m_forceReInstance= true;
        }

        m_meshObject= meshobj;
        
        if (dm) {
                dm->needsFree = 1;
                dm->release(dm);
        }
        return true;
}



bool CcdShapeConstructionInfo::SetProxy(CcdShapeConstructionInfo* shapeInfo)
{
        if (shapeInfo == NULL)
                return false;
        // no support for dynamic change
        assert(IsUnused());
        m_shapeType = PHY_SHAPE_PROXY;
        m_shapeProxy = shapeInfo;
        return true;
}

btCollisionShape* CcdShapeConstructionInfo::CreateBulletShape(btScalar margin, bool useGimpact, bool useBvh)
{
        btCollisionShape* collisionShape = 0;
        btCompoundShape* compoundShape = 0;     

        if (m_shapeType == PHY_SHAPE_PROXY && m_shapeProxy != NULL)
                return m_shapeProxy->CreateBulletShape(margin, useGimpact, useBvh);

        switch (m_shapeType) 
        {
        default:
                break;

        case PHY_SHAPE_BOX:
                collisionShape = new btBoxShape(m_halfExtend);
                collisionShape->setMargin(margin);
                break;

        case PHY_SHAPE_SPHERE:
                collisionShape = new btSphereShape(m_radius);
                collisionShape->setMargin(margin);
                break;

        case PHY_SHAPE_CYLINDER:
                collisionShape = new btCylinderShapeZ(m_halfExtend);
                collisionShape->setMargin(margin);
                break;

        case PHY_SHAPE_CONE:
                collisionShape = new btConeShapeZ(m_radius, m_height);
                collisionShape->setMargin(margin);
                break;

        case PHY_SHAPE_POLYTOPE:
                collisionShape = new btConvexHullShape(&m_vertexArray[0], m_vertexArray.size()/3, 3*sizeof(btScalar));
                collisionShape->setMargin(margin);
                break;

        case PHY_SHAPE_CAPSULE:
                collisionShape = new btCapsuleShapeZ(m_radius, m_height);
                collisionShape->setMargin(margin);
                break;

        case PHY_SHAPE_MESH:
                // Let's use the latest btScaledBvhTriangleMeshShape: it allows true sharing of 
                // triangle mesh information between duplicates => drastic performance increase when 
                // duplicating complex mesh objects. 
                // BUT it causes a small performance decrease when sharing is not required: 
                // 9 multiplications/additions and one function call for each triangle that passes the mid phase filtering
                // One possible optimization is to use directly the btBvhTriangleMeshShape when the scale is 1,1,1
                // and btScaledBvhTriangleMeshShape otherwise.
                if (useGimpact)
                {                               
                                btTriangleIndexVertexArray* indexVertexArrays = new btTriangleIndexVertexArray(
                                                m_polygonIndexArray.size(),
                                                &m_triFaceArray[0],
                                                3*sizeof(int),
                                                m_vertexArray.size()/3,
                                                &m_vertexArray[0],
                                                3*sizeof(btScalar)
                                );
                                
                                btGImpactMeshShape* gimpactShape =  new btGImpactMeshShape(indexVertexArrays);
                                gimpactShape->setMargin(margin);
                                collisionShape = gimpactShape;
                                gimpactShape->updateBound();

                } else
                {
                        if (!m_unscaledShape || m_forceReInstance)
                        {
                        
                                btTriangleIndexVertexArray* indexVertexArrays = 0;

                                if (0.f != m_weldingThreshold1)
                                {
                                        btTriangleMesh* collisionMeshData = new btTriangleMesh(true,false);
                                        collisionMeshData->m_weldingThreshold = m_weldingThreshold1;
                                        bool removeDuplicateVertices=true;
                                        // m_vertexArray not in multiple of 3 anymore, use m_triFaceArray
                                        for(unsigned int i=0; i<m_triFaceArray.size(); i+=3) {
                                                btScalar *bt = &m_vertexArray[3*m_triFaceArray[i]];
                                                btVector3 v1(bt[0], bt[1], bt[2]);
                                                bt = &m_vertexArray[3*m_triFaceArray[i+1]];
                                                btVector3 v2(bt[0], bt[1], bt[2]);
                                                bt = &m_vertexArray[3*m_triFaceArray[i+2]];
                                                btVector3 v3(bt[0], bt[1], bt[2]);
                                                collisionMeshData->addTriangle(v1, v2, v3, removeDuplicateVertices);
                                        }
                                        indexVertexArrays = collisionMeshData;

                                } else
                                {
                                        indexVertexArrays = new btTriangleIndexVertexArray(
                                                        m_polygonIndexArray.size(),
                                                        &m_triFaceArray[0],
                                                        3*sizeof(int),
                                                        m_vertexArray.size()/3,
                                                        &m_vertexArray[0],
                                                        3*sizeof(btScalar));
                                }
                                
                                // this shape will be shared and not deleted until shapeInfo is deleted
                                
                                // for UpdateMesh, reuse the last memory location so instancing wont crash.
                                if(m_unscaledShape) {
                                        DeleteBulletShape(m_unscaledShape, false);
                                        m_unscaledShape->~btBvhTriangleMeshShape();
                                        m_unscaledShape = new(m_unscaledShape) btBvhTriangleMeshShape( indexVertexArrays, true, useBvh );
                                } else {
                                        m_unscaledShape = new btBvhTriangleMeshShape( indexVertexArrays, true, useBvh );
                                }
                                m_forceReInstance= false;
                        } else if (useBvh && m_unscaledShape->getOptimizedBvh() == NULL) {
                                // the existing unscaledShape was not build with Bvh, do it now
                                m_unscaledShape->buildOptimizedBvh();
                        }
                        collisionShape = new btScaledBvhTriangleMeshShape(m_unscaledShape, btVector3(1.0f,1.0f,1.0f));
                        collisionShape->setMargin(margin);
                }
                break;

        case PHY_SHAPE_COMPOUND:
                if (m_shapeArray.size() > 0)
                {
                        compoundShape = new btCompoundShape();
                        for (std::vector<CcdShapeConstructionInfo*>::iterator sit = m_shapeArray.begin();
                                 sit != m_shapeArray.end();
                                 sit++)
                        {
                                collisionShape = (*sit)->CreateBulletShape(margin, useGimpact, useBvh);
                                if (collisionShape)
                                {
                                        collisionShape->setLocalScaling((*sit)->m_childScale);
                                        compoundShape->addChildShape((*sit)->m_childTrans, collisionShape);
                                }
                        }
                        collisionShape = compoundShape;
                }
                break;
        }
        return collisionShape;
}

void CcdShapeConstructionInfo::AddShape(CcdShapeConstructionInfo* shapeInfo)
{
        m_shapeArray.push_back(shapeInfo);
        shapeInfo->AddRef();
}

CcdShapeConstructionInfo::~CcdShapeConstructionInfo()
{
        for (std::vector<CcdShapeConstructionInfo*>::iterator sit = m_shapeArray.begin();
                 sit != m_shapeArray.end();
                 sit++)
        {
                (*sit)->Release();
        }
        m_shapeArray.clear();
        if (m_unscaledShape)
        {
                DeleteBulletShape(m_unscaledShape, true);
        }
        m_vertexArray.clear();
        if (m_shapeType == PHY_SHAPE_MESH && m_meshObject != NULL) 
        {
                std::map<RAS_MeshObject*,CcdShapeConstructionInfo*>::iterator mit = m_meshShapeMap.find(m_meshObject);
                if (mit != m_meshShapeMap.end() && mit->second == this)
                {
                        m_meshShapeMap.erase(mit);
                }
        }
        if (m_shapeType == PHY_SHAPE_PROXY && m_shapeProxy != NULL)
        {
                m_shapeProxy->Release();
        }
}