Method writeMesh

SRC/recorder/VTKHDF_Recorder.cpp:1432–2015

int VTKHDF_Recorder::writeMesh() {
    // 1) Check domain
    if (theDomain == nullptr) {
        opserr << "WARNING: No domain found -- VTKHDF_Recorder::writeMesh\n";
        return -1;
    }

    // Clear old data structures (if you store them as class members)
    theNodeMapping.clear();
    theEleMapping.clear();
    theNodeTags.clear();
    theEleTags.clear();
    theEleClassTags.clear();
    theEleVtkTags.clear();
    theEleVtkOffsets.clear();

    // 2) Gather node info
    NodeIter &theNodes = theDomain->getNodes();
    Node *theNode;
    numNode = 0;  
    maxNDM = 0;   // max spatial dimension
    maxNDF = 0;   // max degrees of freedom

    while ((theNode = theNodes()) != nullptr) {
        int nodeTag = theNode->getTag();

        const Vector &crd = theNode->getCrds();
        if (crd.Size() > maxNDM) {
            maxNDM = crd.Size();
        }

        const Vector &disp = theNode->getTrialDisp();
        if (disp.Size() > maxNDF) {
            maxNDF = disp.Size();
        }

        // Build mapping (nodeTag -> index)
        theNodeMapping[nodeTag] = numNode;
        theNodeTags.push_back(nodeTag);
        numNode++;
    }

    // 3) Gather element info
    ElementIter &theElements = theDomain->getElements();
    Element *theElement;
    numElement = 0;
    numConnectivityIds = 0;
    
    theEleVtkOffsets.push_back(numConnectivityIds);
    while ((theElement = theElements()) != nullptr) {
        int eleTag   = theElement->getTag();
        int classTag = theElement->getClassTag();

        // Check if recognized by our VTK type map
        auto it = vtktypes.find(classTag);
        if (it != vtktypes.end()) {
            int vtkType = it->second;