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; }