mObject spcsToJSON(void) {
SP_ConstraintIter &theSPs = theDomain.getSPs();
SP_Constraint *theSP;
mObject sps;
mArray sp_val;
mValue tag, ntag, dir, ref_val;
char tag_str[15];
sps.clear();
while ((theSP = theSPs()) != 0) {
sp_val.clear();
tag = theSP->getTag();
ntag = theSP->getNodeTag();
dir = theSP->getDOF_Number();
ref_val = theSP->getValue();
sp_val.push_back(ntag);
sp_val.push_back(dir);
sp_val.push_back(ref_val);
sprintf(tag_str, "%d", tag.get_int());
sps[tag_str] = sp_val;
}
return sps;
}
void
LoadPattern::applyLoad(double pseudoTime)
{
// first determine the load factor
if (theSeries != 0 && isConstant != 0) {
loadFactor = theSeries->getFactor(pseudoTime);
loadFactor *= scaleFactor;
}
NodalLoad *nodLoad;
NodalLoadIter &theNodalIter = this->getNodalLoads();
while ((nodLoad = theNodalIter()) != 0)
nodLoad->applyLoad(loadFactor);
ElementalLoad *eleLoad;
ElementalLoadIter &theElementalIter = this->getElementalLoads();
while ((eleLoad = theElementalIter()) != 0)
eleLoad->applyLoad(loadFactor);
SP_Constraint *sp;
SP_ConstraintIter &theIter = this->getSPs();
while ((sp = theIter()) != 0)
sp->applyConstraint(loadFactor);
}
void
MultiSupportPattern::applyLoad(double time)
{
SP_Constraint *sp;
SP_ConstraintIter &theIter = this->getSPs();
while ((sp = theIter()) != 0)
sp->applyConstraint(time);
}
void
MultiSupportPattern::Print(OPS_Stream &s, int flag)
{
s << "MultiSupportPattern tag: " << this->getTag() << endln;
SP_Constraint *sp;
SP_ConstraintIter &theIter = this->getSPs();
while ((sp = theIter()) != 0)
sp->Print(s, flag);
}
SP_Constraint *
LoadPattern::removeSP_Constraint(int tag)
{
TaggedObject *obj = theSPs->removeComponent(tag);
if (obj == 0)
return 0;
SP_Constraint *result = (SP_Constraint *)obj;
result->setDomain(0);
currentGeoTag++;
return result;
}
void
UniformExcitation::setDomain(Domain *theDomain)
{
this->LoadPattern::setDomain(theDomain);
// now we go through and set all the node velocities to be vel0
// for those nodes not fixed in the dirn!
if (vel0 != 0.0) {
SP_ConstraintIter &theSPs = theDomain->getSPs();
SP_Constraint *theSP;
ID constrainedNodes(0);
int count = 0;
while ((theSP=theSPs()) != 0) {
if (theSP->getDOF_Number() == theDof) {
constrainedNodes[count] = theSP->getNodeTag();
count++;
}
}
NodeIter &theNodes = theDomain->getNodes();
Node *theNode;
Vector newVel(1);
int currentSize = 1;
while ((theNode = theNodes()) != 0) {
int tag = theNode->getTag();
if (constrainedNodes.getLocation(tag) < 0) {
int numDOF = theNode->getNumberDOF();
if (numDOF != currentSize)
newVel.resize(numDOF);
newVel = theNode->getVel();
newVel(theDof) = vel0;
theNode->setTrialVel(newVel);
theNode->commitState();
}
}
}
}
void
LoadPattern::setDomain(Domain *theDomain)
{
// if subclass does not implement .. check for 0 pointer
if (theNodalLoads != 0) {
NodalLoad *nodLoad;
NodalLoadIter &theNodalIter = this->getNodalLoads();
while ((nodLoad = theNodalIter()) != 0)
nodLoad->setDomain(theDomain);
ElementalLoad *eleLoad;
ElementalLoadIter &theElementalIter = this->getElementalLoads();
while ((eleLoad = theElementalIter()) != 0)
eleLoad->setDomain(theDomain);
SP_Constraint *theSP;
SP_ConstraintIter &theSpConstraints = this->getSPs();
while ((theSP = theSpConstraints()) != 0)
theSP->setDomain(theDomain);
}
// now we set this load patterns domain
this->DomainComponent::setDomain(theDomain);
}
int
TransformationConstraintHandler::handle(const ID *nodesLast)
{
// first check links exist to a Domain and an AnalysisModel object
Domain *theDomain = this->getDomainPtr();
AnalysisModel *theModel = this->getAnalysisModelPtr();
Integrator *theIntegrator = this->getIntegratorPtr();
if ((theDomain == 0) || (theModel == 0) || (theIntegrator == 0)) {
opserr << "WARNING TransformationConstraintHandler::handle() - ";
opserr << " setLinks() has not been called\n";
return -1;
}
// get number ofelements and nodes in the domain
// and init the theFEs and theDOFs arrays
int numMPConstraints = theDomain->getNumMPs();
// int numSPConstraints = theDomain->getNumSPs();
int numSPConstraints = 0;
SP_ConstraintIter &theSP1s = theDomain->getDomainAndLoadPatternSPs();
SP_Constraint *theSP1;
while ((theSP1 = theSP1s()) != 0)
numSPConstraints++;
numDOF = 0;
ID transformedNode(0, 64);
int i;
// create an ID of constrained node tags in MP_Constraints
ID constrainedNodesMP(0, numMPConstraints);
MP_Constraint **mps =0;
if (numMPConstraints != 0) {
mps = new MP_Constraint *[numMPConstraints];
if (mps == 0) {
opserr << "WARNING TransformationConstraintHandler::handle() - ";
opserr << "ran out of memory for MP_Constraints";
opserr << " array of size " << numMPConstraints << endln;
return -3;
}
MP_ConstraintIter &theMPs = theDomain->getMPs();
MP_Constraint *theMP;
int index = 0;
while ((theMP = theMPs()) != 0) {
int nodeConstrained = theMP->getNodeConstrained();
if (transformedNode.getLocation(nodeConstrained) < 0)
transformedNode[numDOF++] = nodeConstrained;
constrainedNodesMP[index] = nodeConstrained;
mps[index] = theMP;
index++;
}
}
// create an ID of constrained node tags in SP_Constraints
ID constrainedNodesSP(0, numSPConstraints);;
SP_Constraint **sps =0;
if (numSPConstraints != 0) {
sps = new SP_Constraint *[numSPConstraints];
if (sps == 0) {
opserr << "WARNING TransformationConstraintHandler::handle() - ";
opserr << "ran out of memory for SP_Constraints";
opserr << " array of size " << numSPConstraints << endln;
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
return -3;
}
SP_ConstraintIter &theSPs = theDomain->getDomainAndLoadPatternSPs();
SP_Constraint *theSP;
int index = 0;
while ((theSP = theSPs()) != 0) {
int constrainedNode = theSP->getNodeTag();
if (transformedNode.getLocation(constrainedNode) < 0)
transformedNode[numDOF++] = constrainedNode;
constrainedNodesSP[index] = constrainedNode;
sps[index] = theSP;
index++;
}
}
// create an array for the DOF_Groups and zero it
if ((numDOF != 0) && ((theDOFs = new DOF_Group *[numDOF]) == 0)) {
opserr << "WARNING TransformationConstraintHandler::handle() - ";
opserr << "ran out of memory for DOF_Groups";
opserr << " array of size " << numDOF << endln;
return -3;
}
for (i=0; i<numDOF; i++) theDOFs[i] = 0;
//create a DOF_Group for each Node and add it to the AnalysisModel.
// :must of course set the initial IDs
NodeIter &theNod = theDomain->getNodes();
Node *nodPtr;
int numDofGrp = 0;
int count3 = 0;
int countDOF =0;
numConstrainedNodes = 0;
numDOF = 0;
while ((nodPtr = theNod()) != 0) {
DOF_Group *dofPtr = 0;
int nodeTag = nodPtr->getTag();
int numNodalDOF = nodPtr->getNumberDOF();
int loc = -1;
int createdDOF = 0;
loc = constrainedNodesMP.getLocation(nodeTag);
if (loc >= 0) {
TransformationDOF_Group *tDofPtr =
new TransformationDOF_Group(numDofGrp++, nodPtr, mps[loc], this);
createdDOF = 1;
dofPtr = tDofPtr;
// add any SPs
if (numSPConstraints != 0) {
loc = constrainedNodesSP.getLocation(nodeTag);
if (loc >= 0) {
tDofPtr->addSP_Constraint(*(sps[loc]));
for (int i = loc+1; i<numSPConstraints; i++) {
if (constrainedNodesSP(i) == nodeTag)
tDofPtr->addSP_Constraint(*(sps[i]));
}
}
// add the DOF to the array
theDOFs[numDOF++] = dofPtr;
numConstrainedNodes++;
}
}
if (createdDOF == 0) {
loc = constrainedNodesSP.getLocation(nodeTag);
if (loc >= 0) {
TransformationDOF_Group *tDofPtr =
new TransformationDOF_Group(numDofGrp++, nodPtr, this);
int numSPs = 1;
createdDOF = 1;
dofPtr = tDofPtr;
tDofPtr->addSP_Constraint(*(sps[loc]));
// check for more SP_constraints acting on node and add them
for (int i = loc+1; i<numSPConstraints; i++) {
if (constrainedNodesSP(i) == nodeTag) {
tDofPtr->addSP_Constraint(*(sps[i]));
numSPs++;
}
}
// add the DOF to the array
theDOFs[numDOF++] = dofPtr;
numConstrainedNodes++;
countDOF+= numNodalDOF - numSPs;
}
}
// create an ordinary DOF_Group object if no dof constrained
if (createdDOF == 0) {
if ((dofPtr = new DOF_Group(numDofGrp++, nodPtr)) == 0) {
opserr << "WARNING TransformationConstraintHandler::handle() ";
opserr << "- ran out of memory";
opserr << " creating DOF_Group " << i << endln;
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
return -4;
}
countDOF+= numNodalDOF;
}
if (dofPtr == 0)
opserr << "TransformationConstraintHandler::handle() - error in logic\n";
nodPtr->setDOF_GroupPtr(dofPtr);
theModel->addDOF_Group(dofPtr);
}
// create the FE_Elements for the Elements and add to the AnalysisModel
ElementIter &theEle = theDomain->getElements();
Element *elePtr;
FE_Element *fePtr;
numFE = 0;
ID transformedEle(0, 64);
while ((elePtr = theEle()) != 0) {
int flag = 0;
if (elePtr->isSubdomain() == true) {
Subdomain *theSub = (Subdomain *)elePtr;
if (theSub->doesIndependentAnalysis() == true)
flag = 1;
}
if (flag == 0) {
const ID &nodes = elePtr->getExternalNodes();
int nodesSize = nodes.Size();
int isConstrainedNode = 0;
for (int i=0; i<nodesSize; i++) {
int nodeTag = nodes(i);
if (numMPConstraints != 0) {
int loc = constrainedNodesMP.getLocation(nodeTag);
if (loc >= 0) {
isConstrainedNode = 1;
i = nodesSize;
}
}
if (numSPConstraints != 0 && isConstrainedNode == 0) {
int loc = constrainedNodesSP.getLocation(nodeTag);
if (loc >= 0) {
isConstrainedNode = 1;
i = nodesSize;
}
}
}
if (isConstrainedNode == 1) {
transformedEle[numFE++] = elePtr->getTag();
}
}
}
// create an array for the FE_elements and zero it
if ((numFE != 0) && ((theFEs = new FE_Element *[numFE]) == 0)) {
opserr << "WARNING TransformationConstraintHandler::handle() - ";
opserr << "ran out of memory for FE_elements";
opserr << " array of size " << numFE << endln;
return -2;
}
for (i=0; i<numFE; i++) theFEs[i] = 0;
ElementIter &theEle1 = theDomain->getElements();
// int numConstraints = numMPConstraints+numSPConstraints;
int numFeEle = 0;
int numFE = 0;
while ((elePtr = theEle1()) != 0) {
int tag = elePtr->getTag();
if (elePtr->isSubdomain() == true) {
Subdomain *theSub = (Subdomain *)elePtr;
if (theSub->doesIndependentAnalysis() == false) {
if (transformedEle.getLocation(tag) < 0) {
if ((fePtr = new FE_Element(numFeEle, elePtr)) == 0) {
opserr << "WARNING TransformationConstraintHandler::handle()";
opserr << " - ran out of memory";
opserr << " creating FE_Element " << elePtr->getTag() << endln;
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
return -5;
}
} else {
if ((fePtr = new TransformationFE(numFeEle, elePtr)) == 0) {
opserr << "WARNING TransformationConstraintHandler::handle()";
opserr << " - ran out of memory";
opserr << " creating TransformationFE " << elePtr->getTag() << endln;
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
return -6;
}
theFEs[numFE++] = fePtr;
}
numFeEle++;
theModel->addFE_Element(fePtr);
theSub->setFE_ElementPtr(fePtr);
}
} else {
if (transformedEle.getLocation(tag) < 0) {
if ((fePtr = new FE_Element(numFeEle, elePtr)) == 0) {
opserr << "WARNING TransformationConstraintHandler::handle()";
opserr << " - ran out of memory";
opserr << " creating FE_Element " << elePtr->getTag() << endln;
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
return -5;
}
} else {
if ((fePtr = new TransformationFE(numFeEle, elePtr)) == 0) {
opserr << "WARNING TransformationConstraintHandler::handle()";
opserr << " - ran out of memory";
opserr << " creating TransformationFE " << elePtr->getTag() << endln;
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
return -6;
}
theFEs[numFE++] = fePtr;
}
numFeEle++;
theModel->addFE_Element(fePtr);
}
}
theModel->setNumEqn(countDOF);
// set the number of eqn in the model
// now see if we have to set any of the dof's to -3
// int numLast = 0;
if (nodesLast != 0)
for (i=0; i<nodesLast->Size(); i++) {
int nodeID = (*nodesLast)(i);
Node *nodPtr = theDomain->getNode(nodeID);
if (nodPtr != 0) {
DOF_Group *dofPtr = nodPtr->getDOF_GroupPtr();
const ID &id = dofPtr->getID();
// set all the dof values to -3
for (int j=0; j < id.Size(); j++) {
if (id(j) == -2) {
dofPtr->setID(j,-3);
count3++;
} else {
opserr << "WARNING TransformationConstraintHandler::handle() ";
opserr << " - boundary sp constraint in subdomain";
opserr << " this should not be - results suspect \n";
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
}
}
}
}
if (mps != 0) delete [] mps;
if (sps != 0) delete [] sps;
return count3;
}
int
DomainPartitioner::partition(int numParts, bool usingMain, int mainPartitionTag, int specialElementTag)
{
usingMainDomain = usingMain;
mainPartition = mainPartitionTag;
// first we ensure the partitioned domain has numpart subdomains
// with tags 1 through numparts
for (int i=1; i<=numParts; i++) {
if (i != mainPartition) {
Subdomain *subdomainPtr = myDomain->getSubdomainPtr(i);
if (subdomainPtr == 0) {
opserr << "DomainPartitioner::partition - No Subdomain: ";
opserr << i << " exists\n";
return -1;
}
}
}
// we get the ele graph from the domain and partition it
// Graph &theEleGraph = myDomain->getElementGraph();
// theElementGraph = new Graph(myDomain->getElementGraph());
theElementGraph = &(myDomain->getElementGraph());
int theError = thePartitioner.partition(*theElementGraph, numParts);
if (theError < 0) {
opserr << "DomainPartitioner::partition";
opserr << " - the graph partioner failed to partition the ";
opserr << "element graph\n";
return -10+theError;
}
/* print graph */
// opserr << "DomainPartitioner::partition - eleGraph: \n";
// theElementGraph->Print(opserr, 4);
VertexIter &theVertices1 = theElementGraph->getVertices();
Vertex *vertexPtr = 0;
bool moreThanOne = false;
vertexPtr = theVertices1();
int vertexOnePartition = 0;
if (vertexPtr != 0)
vertexOnePartition = vertexPtr->getColor();
while ((moreThanOne == false) && ((vertexPtr = theVertices1()) != 0)) {
int partition = vertexPtr->getColor();
if (partition != vertexOnePartition ) {
moreThanOne = true;
}
}
if (moreThanOne == false) {
opserr <<"DomainPartitioner::partition - too few elements for model to be partitioned\n";
return -1;
}
int specialElementColor = 1;
if (specialElementTag != 0) {
bool found = false;
VertexIter &theVerticesSpecial = theElementGraph->getVertices();
while ((found == false) && ((vertexPtr = theVerticesSpecial()) != 0)) {
int eleTag = vertexPtr->getRef();
if (eleTag == specialElementTag) {
found = true;
int vertexColor = vertexPtr->getColor();
if (vertexColor != 1)
// specialElementColor = vertexColor;
vertexPtr->setColor(1);
}
}
}
// we create empty graphs for the numParts subdomains,
// in the graphs we place the vertices for the elements on the boundaries
// we do not invoke the destructor on the individual graphs as
// this would invoke the destructor on the individual vertices
if (theBoundaryElements != 0)
delete [] theBoundaryElements;
theBoundaryElements = new Graph * [numParts];
if (theBoundaryElements == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - ran out of memory\n";
numPartitions = 0;
return -1;
}
for (int l=0; l<numParts; l++) {
theBoundaryElements[l] = new Graph(2048); // graphs can grow larger; just an estimate
if (theBoundaryElements[l] == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - ran out of memory\n";
numPartitions = 0;
return -1;
}
}
numPartitions = numParts;
// opserr << "DomainPartitioner::partition() - nodes \n";
// we now create a MapOfTaggedObjectStorage to store the NodeLocations
// and create a new NodeLocation for each node; adding it to the map object
theNodeLocations = new MapOfTaggedObjects();
if (theNodeLocations == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - ran out of memory creating MapOfTaggedObjectStorage for node locations\n";
numPartitions = 0;
return -1;
}
NodeIter &theNodes = myDomain->getNodes();
Node *nodePtr;
while ((nodePtr = theNodes()) != 0) {
NodeLocations *theNodeLocation = new NodeLocations(nodePtr->getTag());
if (theNodeLocation == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - ran out of memory creating NodeLocation for node: " << nodePtr->getTag() << endln;
numPartitions = 0;
return -1;
}
if (theNodeLocations->addComponent(theNodeLocation) == false) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - failed to add NodeLocation to Map for Node: " << nodePtr->getTag() << endln;
numPartitions = 0;
return -1;
}
}
//
// we now iterate through the vertices of the element graph
// to see if the vertex is a boundary vertex or not - if it is
// we add to the appropriate graph created above. We also set the
// value the color variable of each of the external nodes connected
// to the element to a value which will indicate that that node will
// have to be added to the subdomain.
//
VertexIter &theVertexIter = theElementGraph->getVertices();
while ((vertexPtr = theVertexIter()) != 0) {
int eleTag = vertexPtr->getRef();
int vertexColor = vertexPtr->getColor();
const ID &adjacency = vertexPtr->getAdjacency();
int size = adjacency.Size();
for (int i=0; i<size; i++) {
Vertex *otherVertex = theElementGraph->getVertexPtr(adjacency(i));
if (otherVertex->getColor() != vertexColor) {
theBoundaryElements[vertexColor-1]->addVertex(vertexPtr,false);
i = size;
}
}
Element *elePtr = myDomain->getElement(eleTag);
const ID &nodes = elePtr->getExternalNodes();
size = nodes.Size();
for (int j=0; j<size; j++) {
int nodeTag = nodes(j);
TaggedObject *theTaggedObject = theNodeLocations->getComponentPtr(nodeTag);
if (theTaggedObject == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - failed to find NodeLocation in Map for Node: " << nodePtr->getTag() << " -- A BUG!!\n";
numPartitions = 0;
return -1;
}
NodeLocations *theNodeLocation = (NodeLocations *)theTaggedObject;
theNodeLocation->addPartition(vertexColor);
}
}
// now go through the MP_Constraints and ensure the retained node is in every
// partition the constrained node is in
MP_ConstraintIter &theMPs = myDomain->getMPs();
MP_Constraint *mpPtr;
while ((mpPtr = theMPs()) != 0) {
int retained = mpPtr->getNodeRetained();
int constrained = mpPtr->getNodeConstrained();
TaggedObject *theRetainedObject = theNodeLocations->getComponentPtr(retained);
TaggedObject *theConstrainedObject = theNodeLocations->getComponentPtr(constrained);
if (theRetainedObject == 0 || theConstrainedObject == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
if (theRetainedObject == 0)
opserr << " - failed to find NodeLocation in Map for Node: " << retained << " -- A BUG!!\n";
if (theConstrainedObject == 0)
opserr << " - failed to find NodeLocation in Map for Node: " << constrained << " -- A BUG!!\n";
numPartitions = 0;
return -1;
}
NodeLocations *theRetainedLocation = (NodeLocations *)theRetainedObject;
NodeLocations *theConstrainedLocation = (NodeLocations *)theConstrainedObject;
ID &theConstrainedNodesPartitions = theConstrainedLocation->nodePartitions;
int numPartitions = theConstrainedNodesPartitions.Size();
for (int i=0; i<numPartitions; i++) {
theRetainedLocation->addPartition(theConstrainedNodesPartitions(i));
}
}
// we now add the nodes,
TaggedObjectIter &theNodeLocationIter = theNodeLocations->getComponents();
TaggedObject *theNodeObject;
while ((theNodeObject = theNodeLocationIter()) != 0) {
NodeLocations *theNodeLocation = (NodeLocations *)theNodeObject;
int nodeTag = theNodeLocation->getTag();
ID &nodePartitions = theNodeLocation->nodePartitions;
int numPartitions = theNodeLocation->numPartitions;
for (int i=0; i<numPartitions; i++) {
int partition = nodePartitions(i);
if (partition != mainPartition) {
Subdomain *theSubdomain = myDomain->getSubdomainPtr(partition);
if (numPartitions == 1) {
Node *nodePtr = myDomain->removeNode(nodeTag);
theSubdomain->addNode(nodePtr);
} else {
Node *nodePtr = myDomain->getNode(nodeTag);
theSubdomain->addExternalNode(nodePtr);
}
}
}
}
// we now move the elements
VertexIter &theVertices = theElementGraph->getVertices();
while ((vertexPtr = theVertices()) != 0) {
// move the element
int partition = vertexPtr->getColor();
if (partition != mainPartition) {
int eleTag = vertexPtr->getRef();
// opserr << "removing ele: " << eleTag << endln;
Element *elePtr = myDomain->removeElement(eleTag);
// opserr << *elePtr;
if (elePtr != 0) {
// opserr << "adding ele - start\n";
Subdomain *theSubdomain = myDomain->getSubdomainPtr(partition);
theSubdomain->addElement(elePtr);
// opserr << "adding ele - done\n";
} else {
opserr << "DomainPartitioner::partioner - element GONE! - eleTag " << eleTag << endln;
}
}
}
// now we go through the load patterns and move NodalLoad
// 1) make sure each subdomain has a copy of the partitioneddomains load patterns.
// 2) move nodal loads
// 3) move SP_Constraints
LoadPatternIter &theLoadPatterns = myDomain->getLoadPatterns();
LoadPattern *theLoadPattern;
while ((theLoadPattern = theLoadPatterns()) != 0) {
int loadPatternTag = theLoadPattern->getTag();
// check that each subdomain has a loadPattern with a similar tag and class tag
for (int i=1; i<=numParts; i++) {
if (i != mainPartition) {
Subdomain *theSubdomain = myDomain->getSubdomainPtr(i);
LoadPattern *loadPatternCopy = theSubdomain->getLoadPattern(loadPatternTag);
if (loadPatternCopy == 0) {
LoadPattern *newLoadPattern = theLoadPattern->getCopy();
if (newLoadPattern == 0) {
opserr << "DomaiPartitioner::partition - out of memory creating LoadPatterns\n";
return -1;
}
theSubdomain->addLoadPattern(newLoadPattern);
}
}
}
// now remove any nodal loads that correspond to internal nodes in a subdomain
// and add them to the appropriate loadpattern in the subdomain
NodalLoadIter &theNodalLoads = theLoadPattern->getNodalLoads();
NodalLoad *theNodalLoad;
while ((theNodalLoad = theNodalLoads()) != 0) {
int nodeTag = theNodalLoad->getNodeTag();
TaggedObject *theTaggedObject = theNodeLocations->getComponentPtr(nodeTag);
if (theTaggedObject == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - failed to find NodeLocation in Map for Node: " << nodeTag << " -- A BUG!!\n";
numPartitions = 0;
return -1;
}
NodeLocations *theNodeLocation = (NodeLocations *)theTaggedObject;
ID &nodePartitions = theNodeLocation->nodePartitions;
int numPartitions = theNodeLocation->numPartitions;
for (int i=0; i<numPartitions; i++) {
int partition = nodePartitions(i);
if (partition != mainPartition) {
if (numPartitions == 1) {
Subdomain *theSubdomain = myDomain->getSubdomainPtr(partition);
theLoadPattern->removeNodalLoad(theNodalLoad->getTag());
if ((theSubdomain->addNodalLoad(theNodalLoad, loadPatternTag)) != true)
opserr << "DomainPartitioner::partition() - failed to add Nodal Load\n";
}
}
}
}
SP_ConstraintIter &theSPs = theLoadPattern->getSPs();
SP_Constraint *spPtr;
while ((spPtr = theSPs()) != 0) {
int nodeTag = spPtr->getNodeTag();
TaggedObject *theTaggedObject = theNodeLocations->getComponentPtr(nodeTag);
if (theTaggedObject == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - failed to find NodeLocation in Map for Node: " << nodeTag << " -- A BUG!!\n";
numPartitions = 0;
return -1;
}
NodeLocations *theNodeLocation = (NodeLocations *)theTaggedObject;
ID &nodePartitions = theNodeLocation->nodePartitions;
int numPartitions = theNodeLocation->numPartitions;
for (int i=0; i<numPartitions; i++) {
int partition = nodePartitions(i);
if (partition != mainPartition) {
Subdomain *theSubdomain = myDomain->getSubdomainPtr(partition);
if (numPartitions == 1)
theLoadPattern->removeSP_Constraint(spPtr->getTag());
int res = theSubdomain->addSP_Constraint(spPtr, loadPatternTag);
if (res < 0)
opserr << "DomainPartitioner::partition() - failed to add SP Constraint\n";
}
}
}
ElementalLoadIter &theLoads = theLoadPattern->getElementalLoads();
ElementalLoad *theLoad;
while ((theLoad = theLoads()) != 0) {
int loadEleTag = theLoad->getElementTag();
SubdomainIter &theSubdomains = myDomain->getSubdomains();
Subdomain *theSub;
bool added = false;
while (((theSub = theSubdomains()) != 0) && (added == false)) {
bool res = theSub->hasElement(loadEleTag);
if (res == true) {
theLoadPattern->removeElementalLoad(theLoad->getTag());
theSub->addElementalLoad(theLoad, loadPatternTag);
if (res < 0)
opserr << "DomainPartitioner::partition() - failed to add ElementalLoad\n";
added = true;
}
}
}
}
// add the single point constraints,
SP_ConstraintIter &theDomainSP = myDomain->getSPs();
SP_Constraint *spPtr;
while ((spPtr = theDomainSP()) != 0) {
int nodeTag = spPtr->getNodeTag();
TaggedObject *theTaggedObject = theNodeLocations->getComponentPtr(nodeTag);
if (theTaggedObject == 0) {
opserr << "DomainPartitioner::partition(int numParts)";
opserr << " - failed to find NodeLocation in Map for Node: " << nodeTag << " -- A BUG!!\n";
numPartitions = 0;
return -1;
}
NodeLocations *theNodeLocation = (NodeLocations *)theTaggedObject;
ID &nodePartitions = theNodeLocation->nodePartitions;
int numPartitions = theNodeLocation->numPartitions;
for (int i=0; i<numPartitions; i++) {
int partition = nodePartitions(i);
if (partition != mainPartition) {
Subdomain *theSubdomain = myDomain->getSubdomainPtr(partition);
if (numPartitions == 1) {
myDomain->removeSP_Constraint(spPtr->getTag());
}
int res = theSubdomain->addSP_Constraint(spPtr);
if (res < 0)
opserr << "DomainPartitioner::partition() - failed to add SP Constraint\n";
}
}
}
// move MP_Constraints - add an MP_Constraint to every partition a constrained node is in
MP_ConstraintIter &moreMPs = myDomain->getMPs();
while ((mpPtr = moreMPs()) != 0) {
int constrained = mpPtr->getNodeConstrained();
TaggedObject *theConstrainedObject = theNodeLocations->getComponentPtr(constrained);
NodeLocations *theConstrainedLocation = (NodeLocations *)theConstrainedObject;
ID &theConstrainedNodesPartitions = theConstrainedLocation->nodePartitions;
int numPartitions = theConstrainedLocation->numPartitions;
for (int i=0; i<numPartitions; i++) {
int partition = theConstrainedNodesPartitions(i);
if (partition != mainPartition) {
Subdomain *theSubdomain = myDomain->getSubdomainPtr(partition);
if (numPartitions == 1)
myDomain->removeMP_Constraint(mpPtr->getTag());
int res = theSubdomain->addMP_Constraint(mpPtr);
if (res < 0)
opserr << "DomainPartitioner::partition() - failed to add MP Constraint\n";
}
}
}
// now we go through all the subdomains and tell them to update
// their analysis for the new layouts
SubdomainIter &theSubDomains = myDomain->getSubdomains();
Subdomain *theSubDomain;
while ((theSubDomain = theSubDomains()) != 0)
theSubDomain->domainChange();
// we invoke change on the PartitionedDomain
myDomain->domainChange();
myDomain->clearElementGraph();
// we are done
partitionFlag = true;
return 0;
}
GSA_Recorder::GSA_Recorder(Domain &theDom,
const char *fileName,
const char *title1,
const char *title2,
const char *title3,
const char *jobno,
const char *initials,
const char *spec,
const char *currency,
const char *length,
const char *force,
const char *temp,
double dT)
: Recorder(RECORDER_TAGS_GSA_Recorder),
theDomain(&theDom), ndm(3), ndf(6), counter(0), deltaT(dT), nextTimeStampToRecord(0.0)
{
// open file
if (theFile.setFile(fileName, OVERWRITE) < 0) {
opserr << "WARNING - GSA_Recorder::GSA_Recorder()";
opserr << " - could not open file " << fileName << endln;
exit(-1);
}
// spit out header data
if (title1 != 0)
theFile << "TITLE\t" << title1;
else
theFile << "TITLE\t" << "No Title";
if (title2 != 0)
theFile << "\t" << title2;
else
theFile << "\t" << "BLANK";
if (title3 != 0)
theFile << "\t" << title3;
else
theFile << "\t" << "BLANK";
if (jobno != 0)
theFile << "\t" << jobno;
else
theFile << "\t" << "0000";
if (initials != 0)
theFile<< "\t" << initials << endln;
else
theFile << "\t" << "ANOTHER\n";
if (spec != 0)
theFile << "SPEC\t" << spec << endln;
if (currency != 0)
theFile << "CURRENCY\t" << currency << endln;
if (length != 0)
theFile << "UNIT_DATA\tLENGTH\t" << length << endln;
if (force != 0)
theFile << "UNIT_DATA\tFORCE\t" << force << endln;
if (temp != 0)
theFile << "UNIT_DATA\tTEMP\t" << temp << endln;
// spit out nodal data
NodeIter &theNodes = theDomain->getNodes();
Node *theNode;
while ((theNode=theNodes()) != 0) {
int nodeTag = theNode->getTag();
theFile << "NODE\t" << nodeTag;
const Vector &crds = theNode->getCrds();
if (crds.Size() != ndm) {
opserr << "WARNING - GSA_Recorder::GSA_Recorder() - node: " << nodeTag ;
opserr << " has invalid number of coordinates, expecting: " << ndm << " got: " << crds.Size() << endln;
exit(-1);
}
const Vector &disp = theNode->getTrialDisp();
if (disp.Size() != ndf) {
opserr << "WARNING - GSA_Recorder::GSA_Recorder() - node: " << nodeTag ;
opserr << " has invalid number of dof, expecting: " << ndf << " got: " << disp.Size() << endln;
exit(-1);
}
for (int i=0; i<ndm; i++)
theFile << "\t" << crds(i);
theFile << endln;
}
// open file and spit out the initial data
SP_ConstraintIter &theSPs = theDomain->getSPs();
SP_Constraint *theSP;
ID theConstrainedNodes(0,6);
ID theSpMatrix(0, 6*ndf);
int numNodesWithSP = 0;
while ((theSP=theSPs()) != 0) {
int nodeTag = theSP->getNodeTag();
int location = theConstrainedNodes.getLocation(nodeTag);
if (location < 0) {
theConstrainedNodes[numNodesWithSP] = nodeTag;
for (int i=0; i<ndf; i++)
theSpMatrix[numNodesWithSP*ndf+i] = 0;
location = numNodesWithSP++;
}
int id = theSP->getDOF_Number();
theSpMatrix[location*ndf + id] = 1;
}
for (int j=0; j<numNodesWithSP; j++) {
theFile << "SPC\t" << theConstrainedNodes[j] << "\t0";
for (int i=0; i<ndf; i++)
theFile << "\t" << theSpMatrix[j*ndf+i];
theFile << endln;
}
ElementIter &theElements = theDomain->getElements();
Element *theElement;
while ((theElement=theElements()) != 0) {
theElement->Print(theFile, -1);
}
}
int
LoadPattern::recvSelf(int cTag, Channel &theChannel, FEM_ObjectBroker &theBroker)
{
// get my current database tag
// NOTE - dbTag equals 0 if not sending to a database OR has not yet been sent
int myDbTag = this->getDbTag();
// into an ID we place all info needed to determine state of LoadPattern
int numNod, numEle, numSPs;
ID lpData(11);
if (theChannel.recvID(myDbTag, cTag, lpData) < 0) {
opserr << "LoadPattern::recvSelf - channel failed to recv the initial ID\n";
return -1;
}
isConstant = lpData(7);
this->setTag(lpData(10));
if (isConstant == 0) { // we must recv the load factor in a Vector
Vector data(2);
if (theChannel.recvVector(myDbTag, cTag, data) < 0) {
opserr << "LoadPattern::recvSelf - channel failed to recv the Vector\n";
return -2;
}
loadFactor = data(0);
scaleFactor = data(1);
}
// read data about the time series
if (lpData(8) != -1) {
if (theSeries == 0) {
theSeries = theBroker.getNewTimeSeries(lpData(8));
} else if (theSeries->getClassTag() != lpData(8)) {
delete theSeries;
theSeries = theBroker.getNewTimeSeries(lpData(8));
}
if (theSeries == 0) {
opserr << "LoadPattern::recvSelf - failed to create TimeSeries\n";
return -3;
}
theSeries->setDbTag(lpData(9));
if (theSeries->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - the TimeSeries failed to recv\n";
return -3;
}
}
/*
if (theChannel.isDatastore() == 1) {
static ID theLastSendTag(1);
if (theChannel.recvID(myDbTag,0,theLastSendTag) == 0)
lastGeoSendTag = theLastSendTag(0);
}
*/
if (lastChannel != theChannel.getTag() || currentGeoTag != lpData(0) || theChannel.isDatastore() == 0) {
// clear out the all the components in the current load pattern
this->clearAll();
lastChannel = theChannel.getTag();
currentGeoTag = lpData(0);
numNod = lpData(1);
numEle = lpData(2);
numSPs = lpData(3);
dbNod = lpData(4);
dbEle = lpData(5);
dbSPs = lpData(6);
//
// now we rebuild the nodal loads
//
// first get the information from the domainData about the nodes
if (numNod != 0) {
ID nodeData(2*numNod);
// now receive the ID about the nodes, class tag and dbTags
if (theChannel.recvID(dbNod, currentGeoTag, nodeData) < 0) {
opserr << "LoadPAttern::recvSelf - channel failed to recv the NodalLoad ID\n";
return -2;
}
// now for each NodalLoad we 1) get a new node of the correct type from the ObjectBroker
// 2) ensure the node exists and set it's dbTag, 3) we invoke recvSelf on this new
// blank node and 4) add this node to the domain
int loc = 0;
for (int i=0; i<numNod; i++) {
int classTag = nodeData(loc);
int dbTag = nodeData(loc+1);
NodalLoad *theNode = theBroker.getNewNodalLoad(classTag);
if (theNode == 0) {
opserr << "LoadPattern::recv - cannot create NodalLoad with classTag " << classTag << endln;
return -2;
}
theNode->setDbTag(dbTag);
if (theNode->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - NodalLoad with dbTag " << dbTag << " failed in recvSelf\n";
return -2;
}
if (this->addNodalLoad(theNode) == false) {
opserr << "LoadPattern::recvSelf - failed adding NodalLoad tagged " << theNode->getTag() << " into LP!\n";
return -3;
}
loc+=2;
}
}
//
// now we rebuild the ElementalLoads .. same as NodalLoads above .. see comments above
//
if (numEle != 0) {
ID eleData(2*numEle);
if (theChannel.recvID(dbEle, currentGeoTag, eleData) < 0) {
opserr << "LoadPattern::recvSelf - channel failed to recv the EleLoad ID\n";
return -2;
}
int loc = 0;
for (int i=0; i<numEle; i++) {
int classTag = eleData(loc);
int dbTag = eleData(loc+1);
ElementalLoad *theEle = theBroker.getNewElementalLoad(classTag);
if (theEle == 0) {
opserr << "LoadPattern::recv - cannot create ElementalLoad with classTag " << classTag << endln;
return -2;
}
theEle->setDbTag(dbTag);
if (theEle->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - Ele with dbTag " << dbTag << " failed in recvSelf\n";
return -2;
}
if (this->addElementalLoad(theEle) == false) {
opserr << "LoadPattern::recvSelf - could not add Ele with tag " << theEle->getTag() << " into LP!\n";
return -3;
}
loc+=2;
}
}
//
// now we rebuild the SP_Constraints .. same as nodes above .. see above if can't understand!!
//
if (numSPs != 0) {
ID spData(2*numSPs);
if (theChannel.recvID(dbSPs, currentGeoTag, spData) < 0) {
opserr << "LoadPattern::recvSelf - channel failed to recv the SP_Constraints ID\n";
return -2;
}
int loc = 0;
for (int i=0; i<numSPs; i++) {
int classTag = spData(loc);
int dbTag = spData(loc+1);
SP_Constraint *theSP = theBroker.getNewSP(classTag);
if (theSP == 0) {
opserr << "LoadPattern::recv - cannot create SP_Constraint with classTag " << classTag << endln;
return -2;
}
theSP->setDbTag(dbTag);
if (theSP->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - SP_Constraint with dbTag " << dbTag << " failed in recvSelf\n";
return -2;
}
if (this->addSP_Constraint(theSP) == false) {
opserr << "LoadPattern::recvSelf - could not add SP_Constraint with tag " << theSP->getTag()
<< " into LP!\n";
return -3;
}
loc+=2;
}
}
// now set the load pattern db count
currentGeoTag = lpData(0);
lastGeoSendTag = currentGeoTag;
} else {
if (theSeries != 0)
if (theSeries->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - the TimeSeries failed to recv\n";
return -3;
}
NodalLoad *theNode;
NodalLoadIter &theNodes = this->getNodalLoads();
while ((theNode = theNodes()) != 0) {
if (theNode->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - node with tag " << theNode->getTag() << " failed in recvSelf\n";
return -7;
}
}
ElementalLoad *theEle;
ElementalLoadIter &theElements = this->getElementalLoads();
while ((theEle = theElements()) != 0) {
if (theEle->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - element with tag " << theEle->getTag() << " failed in recvSelf\n";
return -8;
}
}
SP_Constraint *theSP;
SP_ConstraintIter &theSPs = this->getSPs();
while ((theSP = theSPs()) != 0) {
if (theSP->recvSelf(cTag, theChannel, theBroker) < 0) {
opserr << "LoadPattern::recvSelf - SP_Constraint tagged " << theSP->getTag() << " failed recvSelf\n";
return -9;
}
}
}
// if we get here we are successfull
return 0;
}
int
LoadPattern::sendSelf(int cTag, Channel &theChannel)
{
// get my current database tag
// NOTE - dbTag equals 0 if not sending to a database OR has not yet been sent
int myDbTag = this->getDbTag();
// into an ID we place all info needed to determine state of LoadPattern
int numNodLd, numEleLd, numSPs;
ID lpData(11);
numNodLd = theNodalLoads->getNumComponents();
numEleLd = theElementalLoads->getNumComponents();
numSPs = theSPs->getNumComponents();
lpData(10) = this->getTag();
lpData(0) = currentGeoTag;
lpData(1) = numNodLd;
lpData(2) = numEleLd;
lpData(3) = numSPs;
if (dbNod == 0) {
dbNod = theChannel.getDbTag();
dbEle = theChannel.getDbTag();
dbSPs = theChannel.getDbTag();
}
lpData(4) = dbNod;
lpData(5) = dbEle;
lpData(6) = dbSPs;
lpData(7) = isConstant;
if (theSeries != 0) {
int dbtag = theSeries->getDbTag();
int classtag = theSeries->getClassTag();
if (dbtag == 0) {
dbtag = theChannel.getDbTag();
theSeries->setDbTag(dbtag);
}
lpData(8) = classtag;
lpData(9) = dbtag;
} else
lpData(8) = -1;
// see if we can save sending the vector containing just the load factor
// will happen in parallel if sending the loadPattern .. not in database
if (theChannel.sendID(myDbTag, cTag, lpData) < 0) {
opserr << "LoadPattern::sendSelf - channel failed to send the initial ID\n";
return -1;
}
if (isConstant == 0) {
Vector data(2);
data(0) = loadFactor;
data(1) = scaleFactor;
if (theChannel.sendVector(myDbTag, cTag, data) < 0) {
opserr << "LoadPattern::sendSelf - channel failed to send the Vector\n";
return -2;
}
}
if (theSeries != 0)
if (theSeries->sendSelf(cTag, theChannel) < 0) {
opserr << "LoadPattern::sendSelf - the TimeSeries failed to send\n";
return -3;
}
// now check if data defining the objects in the LoadPAttern needs to be sent
// NOTE THIS APPROACH MAY NEED TO CHANGE FOR VERY LARGE PROBLEMS IF CHANNEL CANNOT
// HANDLE VERY LARGE ID OBJECTS.
/*
if (theChannel.isDatastore() == 1) {
static ID theLastSendTag(1);
if (theChannel.recvID(myDbTag,0,theLastSendTag) == 0)
lastGeoSendTag = theLastSendTag(0);
else
lastGeoSendTag = -1;
}
*/
if (lastChannel != theChannel.getTag() || lastGeoSendTag != currentGeoTag || theChannel.isDatastore() == 0) {
lastChannel = theChannel.getTag();
//
// into an ID we are gonna place the class and db tags for each node so can rebuild
// this ID we then send to the channel
//
// create the ID and get the node iter
if (numNodLd != 0) {
ID nodeData(numNodLd*2);
NodalLoad *theNode;
NodalLoadIter &theNodes = this->getNodalLoads();
int loc =0;
// loop over nodes in domain adding their classTag and dbTag to the ID
while ((theNode = theNodes()) != 0) {
nodeData(loc) = theNode->getClassTag();
int dbTag = theNode->getDbTag();
// if dbTag still 0 get one from Channel;
// if this tag != 0 set the dbTag in node
if (dbTag == 0 && myDbTag != 0) {// go get a new tag and setDbTag in ele if this not 0
dbTag = theChannel.getDbTag();
if (dbTag != 0)
theNode->setDbTag(dbTag);
}
nodeData(loc+1) = dbTag;
loc+=2;
}
// now send the ID
if (theChannel.sendID(dbNod, currentGeoTag, nodeData) < 0) {
opserr << "LoadPattern::sendSelf - channel failed to send the NodalLoads ID\n";
return -4;
}
}
// we do the same for elemental loads as we did for nodal loads above .. see comments above!
if (numEleLd != 0) {
ID elementData(numEleLd*2);
ElementalLoad *theEle;
ElementalLoadIter &theElements = this->getElementalLoads();
int loc = 0;
while ((theEle = theElements()) != 0) {
elementData(loc) = theEle->getClassTag();
int dbTag = theEle->getDbTag();
if (dbTag == 0 && myDbTag != 0) {// go get a new tag and setDbTag in ele if this not 0
dbTag = theChannel.getDbTag();
if (dbTag != 0)
theEle->setDbTag(dbTag);
}
elementData(loc+1) = dbTag;
loc+=2;
}
// now send the ID
if (theChannel.sendID(dbEle, currentGeoTag, elementData) < 0) {
opserr << "Domain::send - channel failed to send the element ID\n";
return -5;
}
}
// we do the same for SP_Constraints as for NodalLoads above .. see comments above!
if (numSPs != 0) {
ID spData(numSPs*2);
SP_Constraint *theSP;
SP_ConstraintIter &theSPs = this->getSPs();
int loc = 0;
while ((theSP = theSPs()) != 0) {
spData(loc) = theSP->getClassTag();
int dbTag = theSP->getDbTag();
if (dbTag == 0 && myDbTag != 0) {// go get a new tag and setDbTag in ele if this not 0
dbTag = theChannel.getDbTag();
if (dbTag != 0)
theSP->setDbTag(dbTag);
}
spData(loc+1) = dbTag;
loc+=2;
}
if (theChannel.sendID(dbSPs, currentGeoTag, spData) < 0) {
opserr << "LoadPAttern::sendSelf - channel failed sending SP_Constraint ID\n";
return -6;
}
}
// set the lst send db tag so we don't have to do all that again
lastGeoSendTag = currentGeoTag;
if (theChannel.isDatastore() == 1) {
static ID theLastSendTag(1);
theLastSendTag(0) = lastGeoSendTag;
theChannel.sendID(myDbTag,0, theLastSendTag);
}
}
// now we invoke sendSelf on all the NodalLoads, ElementalLoads and SP_Constraints
// which have been added to the LoadCase
NodalLoad *theNode;
NodalLoadIter &theNodes = this->getNodalLoads();
while ((theNode = theNodes()) != 0) {
if (theNode->sendSelf(cTag, theChannel) < 0) {
opserr << "LoadPattern::sendSelf - node with tag " << theNode->getTag() << " failed in sendSelf\n";
return -7;
}
}
ElementalLoad *theEle;
ElementalLoadIter &theElements = this->getElementalLoads();
while ((theEle = theElements()) != 0) {
if (theEle->sendSelf(cTag, theChannel) < 0) {
opserr << "LoadPattern::sendSelf - element with tag " << theEle->getTag() << " failed in sendSelf\n";
return -8;
}
}
SP_Constraint *theSP;
SP_ConstraintIter &theSPs = this->getSPs();
while ((theSP = theSPs()) != 0) {
if (theSP->sendSelf(cTag, theChannel) < 0) {
opserr << "LoadPattern::sendSelf - SP_Constraint: " << *theSP << " failed sendSelf\n";
return -9;
}
}
// if we get here we are successfull
return 0;
}
int
RemoveRecorder::elimNode(int theNodeTag, double timeStamp)
{
// remove from domain but do not delete yet!
Node *theNode = theDomain->removeNode(theNodeTag);
// go through all load patterns and remove associated loads and constraints
LoadPatternIter &theLoadPatterns = theDomain->getLoadPatterns();
LoadPattern *thePattern;
while ((thePattern = theLoadPatterns()) != 0) {
// start with nodal laods
NodalLoadIter theLoads = thePattern->getNodalLoads();
NodalLoad *theLoad;
// ID theLoadTags(0,12);
// int cnt=0;
while ((theLoad = theLoads()) != 0) {
int NodeTag = theLoad->getNodeTag();
if (NodeTag == theNodeTag) {
// theLoadTags[cnt] = theLoad->getTag();
// cnt++;
#ifdef MMTDEBUG
opserr<<"identified load pattern "<<theLoad->getTag()<<" acting on node "<<theNode<<endln;
#endif
NodalLoad *theNodalLoad = thePattern->removeNodalLoad(theLoad->getTag());
if (theNodalLoad != 0) {
#ifdef MMTDEBUG
opserr<<"deleting nodal load pattern "<<theLoad->getTag()<<endln;
#endif
delete theNodalLoad;
}
}
}
// for (int i=0; i<cnt; i++) {
// NodalLoad *theNodalLoad = thePattern->removeNodalLoad(theLoadTags[i]);
// if (theNodalLoad != 0) {
// delete theNodalLoad;
// #ifdef MMTDEBUG
// opserr<<"deleting nodal load pattern "<<theLoadTags(i)<<endln;
// #endif
// }
// }
// follow with sp constraints
SP_ConstraintIter &theSPs = thePattern->getSPs();
SP_Constraint *theSP;
while ((theSP = theSPs()) != 0) {
int spNode = theSP->getNodeTag();
if (spNode == theNodeTag) {
// theSPTags[cnt] = theSP->getTag();
// cnt++;
#ifdef MMTDEBUG
opserr<<"identified SP_Constraint "<<theSP->getTag()<<" acting on node "<<spNode<<endln;
#endif
SP_Constraint *theSPConstraint = thePattern->removeSP_Constraint(theSP->getTag());
if (theSPConstraint != 0) {
#ifdef MMTDEBUG
opserr<<"deleting SP_Constraint "<<theSP->getTag()<<endln;
#endif
delete theSPConstraint;
}
}
}
}
// we also have to remove any sp constraints from the domain that do not belong to load patterns (support fixity)
SP_ConstraintIter &theSPs = theDomain->getSPs();
SP_Constraint *theSP;
// ID theSPTags(0,12);
// int cnt=0;
while ((theSP = theSPs()) != 0) {
int spNode = theSP->getNodeTag();
if (spNode == theNodeTag) {
// theSPTags[cnt] = theSP->getTag();
// cnt++;
#ifdef MMTDEBUG
opserr<<"identified SP_Constraint "<<theSP->getTag()<<" acting on node "<<spNode<<endln;
#endif
SP_Constraint *theSPConstraint = theDomain->removeSP_Constraint(theSP->getTag());
if (theSPConstraint != 0) {
#ifdef MMTDEBUG
opserr<<"deleting SP_Constraint "<<theSP->getTag()<<endln;
#endif
delete theSPConstraint;
}
}
}
// for (int i=0; i<cnt; i++) {
// SP_Constraint *theSPconstraint = theDomain->removeSP_Constraint(theSPTags[i]);
// if (theSPconstraint != 0) {
// delete theSPconstraint;
// }
// }
if (theNode != 0) {
// delete theNode;
/////////////////// M.Talaat : Again, avoid recorder trouble
theNode->revertToStart();
}
RemoveRecorder::remNodeList[numRemNodes] = theNode->getTag();
// RemoveRecorder::remNodes[numRemNodes] = theNode;
// RemoveRecorder::numRemNodes ++;
Node **newRemNodes = new Node *[numRemNodes+1];
for (int ii=0; ii<numRemNodes; ii++)
newRemNodes[ii] = remNodes[ii];
newRemNodes[numRemNodes] = theNode;
if (remNodes != 0)
delete [] remNodes;
remNodes = newRemNodes;
numRemNodes++;
// now give us some notice of what happened
if (fileName != 0)
theFile<<timeStamp<<" Node "<<theNode->getTag()<<"\n";
if (echoTimeFlag == true)
opserr<<"Node "<<theNode->getTag()<<" removed, Time/Load Factor = " <<timeStamp<<endln;
return 0;
}
int
ActorSubdomain::run(void)
{
static Vector theVect(4);
static Vector theVect1(1);
bool exitYet = false;
int res = 0;
while (exitYet == false) {
int action;
res = this->recvID(msgData);
if (res != 0) {
opserr << "ActorSubdomain::run - error receiving msgData\n";
exitYet = true;
action = ShadowActorSubdomain_DIE;
} else {
action = msgData(0);
}
bool change;
int theType, theOtherType, tag, dbTag, loadPatternTag;
int startTag, endTag, axisDirn, numSP, i, numMode, dof;
Element *theEle;
Node *theNod;
SP_Constraint *theSP;
MP_Constraint *theMP;
LoadPattern *theLoadPattern;
NodalLoad *theNodalLoad;
ElementalLoad *theElementalLoad;
DomainDecompositionAnalysis *theDDAnalysis;
const Matrix *theMatrix;
const Vector *theVector;
Matrix *theM;
Vector *theV;
ID *theI, *theNodeTags, *theEleTags;
PartitionedModelBuilder *theBuilder;
IncrementalIntegrator *theIntegrator;
EquiSolnAlgo *theAlgorithm;
LinearSOE *theSOE;
EigenSOE *theEigenSOE;
ConvergenceTest *theTest;
Recorder *theRecorder;
bool res, generalized, findSmallest;
double doubleRes;
int intRes;
NodeResponseType nodeResponseType;
Parameter *theParameter;
int argc;
char **argv;
char *allResponseArgs;
char *currentLoc;
int argLength, msgLength;
Message theMessage;
const ID *theID;
// opserr << "ActorSubdomain ACTION: " << action << endln;
switch (action) {
case ShadowActorSubdomain_setTag:
tag = msgData(1); // subdomain tag
this->setTag(tag);
this->Actor::setCommitTag(tag);
break;
case ShadowActorSubdomain_analysisStep:
this->recvVector(theVect);
this->analysisStep(theVect(0));
break;
case ShadowActorSubdomain_eigenAnalysis:
numMode = msgData(1);
if (msgData(2) == 0)
generalized = true;
else
generalized = false;
if (msgData(3) == 0)
findSmallest = true;
else
findSmallest = false;
this->eigenAnalysis(numMode, generalized, findSmallest);
break;
/*
case ShadowActorSubdomain_buildSubdomain:
theType = msgData(1);
tag = msgData(3); // subdomain tag
this->setTag(tag);
tag = msgData(2); // numSubdomains
theBuilder = theBroker->getPtrNewPartitionedModelBuilder(*this,
theType);
this->recvObject(*theBuilder);
this->buildSubdomain(tag, *theBuilder);
break;
*/
case ShadowActorSubdomain_getRemoteData:
theID = &(this->getExternalNodes());
msgData(0) = theID->Size();
msgData(1) = this->getNumDOF();
this->sendID(msgData);
if (theID->Size() != 0)
this->sendID(*theID);
break;
case ShadowActorSubdomain_getCost:
theVect(0) = this->getCost(); // have to use [] for Sun's CC!
this->sendVector(theVect);
break;
case ShadowActorSubdomain_addElement:
theType = msgData(1);
dbTag = msgData(2);
theEle = theBroker->getNewElement(theType);
if (theEle != 0) {
theEle->setDbTag(dbTag);
this->recvObject(*theEle);
bool result = this->addElement(theEle);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
/*
this->recvID(msgData);
opserr << "ActorSubdomain::addElement() : " << msgData;
msgData(0) = 1;
msgData(1) = 2;
msgData(2) = 3;
msgData(3) = 4;
this->sendID(msgData);
*/
break;
case ShadowActorSubdomain_hasNode:
theType = msgData(1);
res = this->hasNode(theType);
if (res == true)
msgData(0) = 0;
else
msgData(0) = -1;
this->sendID(msgData);
break;
case ShadowActorSubdomain_hasElement:
theType = msgData(1);
res = this->hasElement(theType);
if (res == true)
msgData(0) = 0;
else
msgData(0) = -1;
this->sendID(msgData);
break;
case ShadowActorSubdomain_addNode:
theType = msgData(1);
dbTag = msgData(2);
theNod = theBroker->getNewNode(theType);
if (theNod != 0) {
theNod->setDbTag(dbTag);
this->recvObject(*theNod);
bool result = this->addNode(theNod);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
// opserr << "ActorSubdomain::add node: " << *theNod;
break;
case ShadowActorSubdomain_addExternalNode:
theType = msgData(1);
dbTag = msgData(2);
theNod = theBroker->getNewNode(theType);
if (theNod != 0) {
theNod->setDbTag(dbTag);
this->recvObject(*theNod);
bool result = this->Subdomain::addExternalNode(theNod);
delete theNod;
/*
Node *dummy = new Node(*theNod);
delete theNod;
opserr << *dummy;
opserr << dummy->getMass();
*/
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_addSP_Constraint:
theType = msgData(1);
dbTag = msgData(2);
theSP = theBroker->getNewSP(theType);
if (theSP != 0) {
theSP->setDbTag(dbTag);
this->recvObject(*theSP);
bool result = this->addSP_Constraint(theSP);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_addSP_ConstraintAXIS:
axisDirn = msgData(1);
theI = new ID(msgData(2));
theV = new Vector(2);
SP_Constraint_SetNextTag(msgData(3));
endTag = 0;
this->recvID(*theI);
this->recvVector(*theV);
msgData(0) = 0;
numSP = this->addSP_Constraint(axisDirn, (*theV)(0), *theI, (*theV)(1));
endTag = SP_Constraint_GetNextTag();
msgData(1) = numSP;
msgData(2) = endTag;
this->domainChange();
this->sendID(msgData);
delete theV;
delete theI;
/* DONT BOTHER SENDING
if (numSP > 0) {
theI = new ID(numSP);
for (i = 0; i<numSP; i++) {
theSP = this->getSP_Constraint(i+startTag);
(*theI)(i) = theSP->getClassTag();
}
this->sendID(*theI);
opserr << "Actor: sent: " << *theI;
for (i = 0; i<numSP; i++) {
theSP = this->getSP_Constraint(i+startTag);
if (theSP != 0)
this->sendObject(*theSP);
else
opserr << "ActorSubdomain::addSP_AXIS :: PROBLEMS\n";
}
delete theI;
}
opserr << "ActorSubdomain::addSP_AXIS :: DONE\n";
*/
break;
case ShadowActorSubdomain_addMP_Constraint:
theType = msgData(1);
dbTag = msgData(2);
theMP = theBroker->getNewMP(theType);
if (theMP != 0) {
theMP->setDbTag(dbTag);
this->recvObject(*theMP);
bool result = this->addMP_Constraint(theMP);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_addLoadPattern:
theType = msgData(1);
dbTag = msgData(2);
theLoadPattern = theBroker->getNewLoadPattern(theType);
if (theLoadPattern != 0) {
theLoadPattern->setDbTag(dbTag);
this->recvObject(*theLoadPattern);
bool result = this->addLoadPattern(theLoadPattern);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_addNodalLoadToPattern:
theType = msgData(1);
dbTag = msgData(2);
loadPatternTag = msgData(3);
theNodalLoad = theBroker->getNewNodalLoad(theType);
if (theNodalLoad != 0) {
theNodalLoad->setDbTag(dbTag);
this->recvObject(*theNodalLoad);
bool result = this->addNodalLoad(theNodalLoad, loadPatternTag);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_addElementalLoadToPattern:
theType = msgData(1);
dbTag = msgData(2);
loadPatternTag = msgData(3);
theElementalLoad = theBroker->getNewElementalLoad(theType);
if (theElementalLoad != 0) {
theElementalLoad->setDbTag(dbTag);
this->recvObject(*theElementalLoad);
bool result = this->addElementalLoad(theElementalLoad,
loadPatternTag);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_addSP_ConstraintToPattern:
theType = msgData(1);
dbTag = msgData(2);
loadPatternTag = msgData(3);
theSP = theBroker->getNewSP(theType);
if (theSP != 0) {
theSP->setDbTag(dbTag);
this->recvObject(*theSP);
bool result = this->addSP_Constraint(theSP, loadPatternTag);
if (result == true)
msgData(0) = 0;
else
msgData(0) = -1;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_removeElement:
tag = msgData(1);
theEle = this->removeElement(tag);
if (theEle != 0)
msgData(0) = theEle->getClassTag();
else
msgData(0) = -1;
this->sendID(msgData);
if (theEle != 0) {
this->sendObject(*theEle);
delete theEle;
}
msgData(0) = 0;
break;
case ShadowActorSubdomain_removeNode:
tag = msgData(1);
theNod = this->removeNode(tag);
if (theNod != 0)
msgData(0) = theNod->getClassTag();
else
msgData(0) = -1;
this->sendID(msgData);
if (theNod != 0) {
this->sendObject(*theNod);
delete theNod;
}
msgData(0) = 0;
break;
case ShadowActorSubdomain_removeSP_Constraint:
tag = msgData(1);
theSP = this->removeSP_Constraint(tag);
break;
case ShadowActorSubdomain_removeSP_ConstraintNoTag:
tag = msgData(1);
dof = msgData(2);
loadPatternTag = msgData(3);
msgData(0) = this->removeSP_Constraint(tag, dof, loadPatternTag);
this->sendID(msgData);
break;
case ShadowActorSubdomain_removeMP_Constraint:
tag = msgData(1);
theMP = this->removeMP_Constraint(tag);
break;
case ShadowActorSubdomain_removeLoadPattern:
tag = msgData(1);
theLoadPattern = this->removeLoadPattern(tag);
break;
case ShadowActorSubdomain_removeNodalLoadFromPattern:
tag = msgData(1);
theType = msgData(2);
theNodalLoad = this->removeNodalLoad(tag, theType);
break;
case ShadowActorSubdomain_removeElementalLoadFromPattern:
tag = msgData(1);
theType = msgData(2);
theElementalLoad = this->removeElementalLoad(tag, theType);
break;
case ShadowActorSubdomain_removeSP_ConstraintFromPattern:
tag = msgData(1);
theType = msgData(2);
theSP = this->removeSP_Constraint(tag, theType);
break;
case ShadowActorSubdomain_getElement:
tag = msgData(1);
theEle = this->getElement(tag);
if (theEle != 0)
msgData(0) = theEle->getClassTag();
else
msgData(0) = -1;
this->sendID(msgData);
if (theEle != 0) {
this->sendObject(*theEle);
}
msgData(0) = 0;
break;
case ShadowActorSubdomain_getNode:
tag = msgData(1);
theNod = this->getNode(tag);
if (theNod != 0)
msgData(0) = theNod->getClassTag();
else
msgData(0) = -1;
this->sendID(msgData);
if (theNod != 0) {
this->sendObject(*theNod);
}
msgData(0) = 0;
break;
case ShadowActorSubdomain_Print:
this->Print(opserr, msgData(3));
this->sendID(msgData);
break;
case ShadowActorSubdomain_PrintNodeAndEle:
theNodeTags = 0;
theEleTags = 0;
if (msgData(1) != 0) {
theNodeTags = new ID(msgData(1));
this->recvID(*theNodeTags);
}
if (msgData(2) != 0) {
theEleTags = new ID(msgData(2));
this->recvID(*theEleTags);
}
this->Print(opserr, theNodeTags, theEleTags, msgData(3));
if (theNodeTags != 0)
delete theNodeTags;
if (theEleTags != 0)
delete theEleTags;
this->sendID(msgData);
break;
case ShadowActorSubdomain_applyLoad:
this->recvVector(theVect);
this->applyLoad(theVect(0));
break;
case ShadowActorSubdomain_setCommittedTime:
this->recvVector(theVect);
this->setCurrentTime(theVect(0));
this->setCommittedTime(theVect(0));
break;
case ShadowActorSubdomain_setLoadConstant:
this->setLoadConstant();
break;
case ShadowActorSubdomain_update:
this->update();
break;
case ShadowActorSubdomain_updateTimeDt:
this->updateTimeDt();
break;
case ShadowActorSubdomain_computeNodalResponse:
tag = msgData(1);
if (lastResponse == 0)
lastResponse = new Vector(tag);
else if (lastResponse->Size() != tag) {
delete lastResponse;
lastResponse = new Vector(tag);
}
this->recvVector(*lastResponse);
this->computeNodalResponse();
break;
case ShadowActorSubdomain_record:
this->record();
break;
case ShadowActorSubdomain_commit:
this->commit();
break;
case ShadowActorSubdomain_revertToLastCommit:
this->revertToLastCommit();
break;
case ShadowActorSubdomain_revertToStart:
this->revertToStart();
this->sendID(msgData);
break;
case ShadowActorSubdomain_addRecorder:
theType = msgData(1);
theRecorder = theBroker->getPtrNewRecorder(theType);
if (theRecorder != 0) {
this->recvObject(*theRecorder);
this->addRecorder(*theRecorder);
}
break;
case ShadowActorSubdomain_removeRecorders:
this->removeRecorders();
this->barrierCheck(1);
break;
case ShadowActorSubdomain_removeRecorder:
theType = msgData(1);
this->removeRecorder(theType);
break;
case ShadowActorSubdomain_wipeAnalysis:
this->wipeAnalysis();
break;
case ShadowActorSubdomain_setDomainDecompAnalysis:
theType = msgData(1);
theDDAnalysis =
theBroker->getNewDomainDecompAnalysis(theType, *this);
if (theDDAnalysis != 0) {
this->recvObject(*theDDAnalysis);
this->setDomainDecompAnalysis(*theDDAnalysis);
msgData(0) = 0;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_setAnalysisAlgorithm:
theType = msgData(1);
theAlgorithm = theBroker->getNewEquiSolnAlgo(theType);
if (theAlgorithm != 0) {
this->recvObject(*theAlgorithm);
this->setAnalysisAlgorithm(*theAlgorithm);
msgData(0) = 0;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_setAnalysisIntegrator:
theType = msgData(1);
theIntegrator = theBroker->getNewIncrementalIntegrator(theType);
if (theIntegrator != 0) {
this->recvObject(*theIntegrator);
this->setAnalysisIntegrator(*theIntegrator);
msgData(0) = 0;
} else
msgData(0) = -1;
this->sendID(msgData);
break;
case ShadowActorSubdomain_setAnalysisLinearSOE:
theType = msgData(1);
theSOE = theBroker->getNewLinearSOE(theType);
if (theSOE != 0) {
this->recvObject(*theSOE);
this->setAnalysisLinearSOE(*theSOE);
msgData(0) = 0;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_setAnalysisEigenSOE:
theType = msgData(1);
theEigenSOE = theBroker->getNewEigenSOE(theType);
if (theEigenSOE != 0) {
this->recvObject(*theEigenSOE);
this->setAnalysisEigenSOE(*theEigenSOE);
msgData(0) = 0;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_setAnalysisConvergenceTest:
theType = msgData(1);
theTest = theBroker->getNewConvergenceTest(theType);
if (theTest != 0) {
this->recvObject(*theTest);
this->setAnalysisConvergenceTest(*theTest);
msgData(0) = 0;
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_domainChange:
this->domainChange();
tag = this->getNumDOF();
if (tag != 0) {
if (lastResponse == 0)
lastResponse = new Vector(tag);
else if (lastResponse->Size() != tag) {
delete lastResponse;
lastResponse = new Vector(tag);
}
}
break;
case ShadowActorSubdomain_getDomainChangeFlag:
change = this->getDomainChangeFlag();
if (change == true)
msgData(0) = 0;
else
msgData(0) = 1;
this->sendID(msgData);
break;
case ShadowActorSubdomain_clearAnalysis:
// this->clearAnalysis();
break;
/*
case 50:
const Matrix *theMatrix1 = &(this->getStiff());
this->sendMatrix(*theMatrix1);
break;
case 51:
const Matrix *theMatrix2 = &(this->getDamp());
this->sendMatrix(*theMatrix2);
break;
case 52:
const Matrix *theMatrix3 = &(this->getMass());
this->sendMatrix(*theMatrix3);
break;
*/
case ShadowActorSubdomain_getTang:
theMatrix = &(this->getTang());
this->sendMatrix(*theMatrix);
break;
case ShadowActorSubdomain_getResistingForce:
theVector = &(this->getResistingForce());
this->sendVector(*theVector);
break;
case ShadowActorSubdomain_computeTang:
tag = msgData(1);
this->setTag(tag);
this->computeTang();
break;
case ShadowActorSubdomain_computeResidual:
this->computeResidual();
break;
case ShadowActorSubdomain_clearAll:
this->clearAll();
this->sendID(msgData);
break;
case ShadowActorSubdomain_getNodeDisp:
tag = msgData(1); // nodeTag
dbTag = msgData(2); // dof
doubleRes = this->getNodeDisp(tag, dbTag, intRes);
msgData(0) = intRes;
this->sendID(msgData);
if (intRes == 0) {
theV = new Vector(1);
(*theV)(0) = doubleRes;
this->sendVector(*theV);
delete theV;
}
break;
case ShadowActorSubdomain_setMass:
tag = msgData(1); // nodeTag
dbTag = msgData(2); // noRows
theOtherType = msgData(3); // noRows
theM = new Matrix(dbTag, theOtherType);
this->recvMatrix(*theM);
intRes = this->setMass(*theM, tag);
delete theM;
msgData(0) = intRes;
this->sendID(msgData);
break;
case ShadowActorSubdomain_getNodeResponse:
tag = msgData(1); // nodeTag
nodeResponseType = (NodeResponseType)msgData(2);
theVector = this->getNodeResponse(tag, nodeResponseType);
if (theVector == 0)
msgData(0) = 0;
else {
msgData(0) = 1;
msgData(1) = theVector->Size();
}
this->sendID(msgData);
if (theVector != 0)
this->sendVector(*theVector);
break;
case ShadowActorSubdomain_getElementResponse:
tag = msgData(1); // eleTag
argc = msgData(2);
msgLength = msgData(3);
if (msgLength == 0) {
opserr << "ElementRecorder::recvSelf() - 0 sized string for responses\n";
return -1;
}
allResponseArgs = new char[msgLength];
if (allResponseArgs == 0) {
opserr << "ElementRecorder::recvSelf() - out of memory\n";
return -1;
}
theMessage.setData(allResponseArgs, msgLength);
if (this->recvMessage(theMessage) < 0) {
opserr << "ElementRecorder::recvSelf() - failed to recv message\n";
return -1;
}
//
// now break this single array into many
//
argv = new char *[argc];
if (argv == 0) {
opserr << "ElementRecorder::recvSelf() - out of memory\n";
return -1;
}
currentLoc = allResponseArgs;
for (int j=0; j<argc; j++) {
argv[j] = currentLoc;
argLength = strlen(currentLoc)+1;
currentLoc += argLength;
}
theVector = this->getElementResponse(tag, (const char**)argv, argc);
delete [] argv;
delete [] allResponseArgs;
if (theVector == 0)
msgData(0) = 0;
else {
msgData(0) = 1;
msgData(1) = theVector->Size();
}
this->sendID(msgData);
if (theVector != 0)
this->sendVector(*theVector);
break;
case ShadowActorSubdomain_calculateNodalReactions:
if (msgData(0) == 0)
this->calculateNodalReactions(true);
else
this->calculateNodalReactions(false);
break;
case ShadowActorSubdomain_setRayleighDampingFactors:
theV = new Vector(4);
this->recvVector(*theV);
intRes = this->Subdomain::setRayleighDampingFactors((*theV)(0), (*theV)(1), (*theV)(2), (*theV)
(3));
delete theV;
break;
case ShadowActorSubdomain_addParameter:
theType = msgData(1);
dbTag = msgData(2);
theParameter = theBroker->getParameter(theType);
if (theParameter != 0) {
theParameter->setDbTag(dbTag);
this->recvObject(*theParameter);
//bool result = true;
bool result = this->addParameter(theParameter);
if (result == true)
msgData(0) = 0;
else {
opserr << "Actor::addParameter - FAILED\n";
msgData(0) = -1;
}
} else
msgData(0) = -1;
break;
case ShadowActorSubdomain_removeParameter:
theType = msgData(1);
this->removeParameter(theType);
this->sendID(msgData);
break;
case ShadowActorSubdomain_updateParameterINT:
theType = msgData(1); // tag
dbTag = msgData(2); // value
msgData(0) = this->Domain::updateParameter(theType, dbTag);
this->sendID(msgData);
break;
case ShadowActorSubdomain_updateParameterDOUBLE:
theType = msgData(1); // tag
this->recvVector(theVect1);
msgData(0) = this->Domain::updateParameter(theType, theVect1(0));
this->sendID(msgData);
break;
case ShadowActorSubdomain_DIE:
exitYet = true;
break;
default:
opserr << "ActorSubdomain::invalid action " << action << "received\n";
msgData(0) = -1;
}
// opserr << "DONE ACTION: " << action << endln;
}
// this->sendID(msgData);
return 0;
}
