void _GeneralSwarm_Init(
void* swarm,
EscapedRoutine* escapedRoutine )
{
GeneralSwarm* self = (GeneralSwarm*)swarm;
GlobalParticle globalParticle;
self->swarmAdvector = NULL; /* If we're using a SwarmAdvector, it will 'attach' itself later on. */
self->escapedRoutine = escapedRoutine;
self->particleCoordVariable = Swarm_NewVectorVariable( self, (Name)"Position", GetOffsetOfMember( globalParticle, coord ),
Variable_DataType_Double,
self->dim,
"PositionX",
"PositionY",
"PositionZ" );
LiveComponentRegister_Add( LiveComponentRegister_GetLiveComponentRegister(), (Stg_Component*)self->particleCoordVariable->variable );
LiveComponentRegister_Add( LiveComponentRegister_GetLiveComponentRegister(), (Stg_Component*)self->particleCoordVariable );
/* init members */
self->previousIntSwarmMap = NULL;
/* lets init this guy with one spot for convenience */
self->intSwarmMapList = List_New("intSwarmMapList");
List_SetItemSize(self->intSwarmMapList, sizeof(SwarmMap*));
}
void _IntegrationPointsSwarm_Init(
void* swarm,
FeMesh* mesh )
{
IntegrationPointsSwarm* self = (IntegrationPointsSwarm*)swarm;
LocalParticle localParticle;
IntegrationPoint particle;
self->mesh = mesh;
self->weightVariable = Swarm_NewScalarVariable( self, (Name)"Weight", GetOffsetOfMember( particle , weight ),
Variable_DataType_Double );
LiveComponentRegister_Add( LiveComponentRegister_GetLiveComponentRegister(), (Stg_Component*)self->weightVariable );
LiveComponentRegister_Add( LiveComponentRegister_GetLiveComponentRegister(), (Stg_Component*)self->weightVariable->variable );
self->localCoordVariable = Swarm_NewVectorVariable( self, (Name)"LocalElCoord", GetOffsetOfMember( localParticle , xi ),
Variable_DataType_Double, self->dim, "Xi", "Eta", "Zeta" );
LiveComponentRegister_Add( LiveComponentRegister_GetLiveComponentRegister(), (Stg_Component*)self->localCoordVariable );
LiveComponentRegister_Add( LiveComponentRegister_GetLiveComponentRegister(), (Stg_Component*)self->localCoordVariable->variable );
/* _Construct calls _Swarm_Init */
/* Lock down the extension manager.
* It doesn't make sense for the IntegrationPointsSwarm to allow IntegrationPoints to be extended
* This means attempts to extend integration points are firewalled to pickup errors.
* -- Alan 20060506
*/
ExtensionManager_SetLockDown( self->particleExtensionMgr, True );
self->swarmsMappedTo = List_New("swarmsMappedTo");
List_SetItemSize(self->swarmsMappedTo, sizeof(SwarmMap*));
}
void ListSuite_Setup( ListSuiteData* data ) {
Index idx;
data->list = List_New();
List_SetItemSize( data->list, sizeof(int) );
for( idx = 0; idx < NUM_ITEMS; idx++ ) {
data->arrayData[idx] = idx;
}
}
Mesh* _Mesh_New( MESH_DEFARGS ) {
Mesh* self;
/* Allocate memory */
assert( _sizeOfSelf >= sizeof(Mesh) );
self = (Mesh*)_Stg_Component_New( STG_COMPONENT_PASSARGS );
self->topo = (MeshTopology*)IGraph_New( "" );
self->vertices = NULL;
self->verticesVariable = NULL;
self->vGlobalIdsVar = NULL;
self->verticesgid = NULL;
self->e_n = NULL;
self->enMapVar = NULL;
self->elgid = NULL;
self->eGlobalIdsVar = NULL;
self->vars = List_New();
List_SetItemSize( self->vars, sizeof(MeshVariable*) );
self->minSep = 0.0;
self->minAxialSep = NULL;
self->minLocalCrd = NULL;
self->maxLocalCrd = NULL;
self->minDomainCrd = NULL;
self->maxDomainCrd = NULL;
self->minGlobalCrd = NULL;
self->maxGlobalCrd = NULL;
self->algorithms = Mesh_Algorithms_New( "", NULL );
Mesh_Algorithms_SetMesh( self->algorithms, self );
self->nElTypes = 0;
self->elTypes = NULL;
self->elTypeMap = NULL;
self->topoDataSizes = UIntMap_New();
self->topoDataInfos = NULL;
self->topoDatas = NULL;
self->info = ExtensionManager_New_OfExistingObject( "mesh_info", self );
self->vertGridId = (unsigned)-1;
self->elGridId = (unsigned)-1;
self->periodicId = (unsigned)-1;
self->localOriginId = (unsigned)-1;
self->localRangeId = (unsigned)-1;
self->generator = NULL;
self->emReg = NULL;
self->isDeforming = False;
self->isRegular = False;
self->parentMesh = NULL;
return self;
}
