Method vtu

SRC/recorder/PVDRecorder.cpp:284–356

}

int
PVDRecorder::vtu()
{
    if (theDomain == 0) {
	opserr << "WARNING: failed to get domain -- PVDRecorder::vtu\n";
	return -1;
    }
    // get node ndf
    NodeIter& theNodes = theDomain->getNodes();
    Node* theNode = 0;
    int nodendf = 0;
    while ((theNode = theNodes()) != 0) {
	if(nodendf < theNode->getNumberDOF()) {
	    nodendf = theNode->getNumberDOF();
	}
    }
    if (nodendf < 3) {
	nodendf = 3;
    } else if (nodendf > 3) {
        nodendf = 3;
    }

    // get parts
    this->getParts();

    // get background mesh
    VInt gtags;
    TaggedObjectIter& meshes = OPS_getAllMesh();
    Mesh* mesh = 0;
    while((mesh = dynamic_cast<Mesh*>(meshes())) != 0) {
	ParticleGroup* group = dynamic_cast<ParticleGroup*>(mesh);
	if (group == 0) {
	    continue;
	}
        gtags.push_back(group->getTag());
    }


    // part 0: all nodes
    ID partno(0, (int)parts.size()+(int)gtags.size()+1);
    partno[0] = 0;
    if (this->savePart0(nodendf) < 0) {
        return -1;
    }

    // particle parts
    for (int i=0; i<(int)gtags.size(); ++i) {
        partno[1+i] = 1+i;
        if (this->savePartParticle(1+i, gtags[i],nodendf) < 0) {
            return -1;
        }
    }

    // save other parts
    // int index = 1;
    for(std::map<int,ID>::iterator it=parts.begin(); it!=parts.end(); it++) {
	// int& no = partnum[it->first];
	// if (no == 0) {