vtkpfemexportmodule.C

Go to the documentation of this file.

/*
 *
 *                 #####    #####   ######  ######  ###   ###
 *               ##   ##  ##   ##  ##      ##      ## ### ##
 *              ##   ##  ##   ##  ####    ####    ##  #  ##
 *             ##   ##  ##   ##  ##      ##      ##     ##
 *            ##   ##  ##   ##  ##      ##      ##     ##
 *            #####    #####   ##      ######  ##     ##
 *
 *
 *             OOFEM : Object Oriented Finite Element Code
 *
 *               Copyright (C) 1993 - 2025   Borek Patzak
 *
 *
 *
 *       Czech Technical University, Faculty of Civil Engineering,
 *   Department of Structural Mechanics, 166 29 Prague, Czech Republic
 *
 *  This library is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU Lesser General Public
 *  License as published by the Free Software Foundation; either
 *  version 2.1 of the License, or (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *  Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public
 *  License along with this library; if not, write to the Free Software
 *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */
 
#include "vtkpfemexportmodule.h"
#include "element.h"
#include "gausspoint.h"
#include "timestep.h"
#include "engngm.h"
#include "node.h"
#include "dof.h"
#include "materialinterface.h"
#include "mathfem.h"
#include "cltypes.h"
#include "material.h"
#include "classfactory.h"
#include "crosssection.h"
#include "unknownnumberingscheme.h"
 
#ifdef __PFEM_MODULE
 #include "pfem/pfemparticle.h"
#endif
 
#include <string>
#include <sstream>
#include <ctime>
 
namespace oofem {
REGISTER_ExportModule(VTKPFEMXMLExportModule)
 
 
// Begin VTKPFEMXMLExportModule
 
VTKPFEMXMLExportModule::VTKPFEMXMLExportModule(int n, EngngModel *e) : VTKBaseExportModule(n, e) {}
 
VTKPFEMXMLExportModule::~VTKPFEMXMLExportModule() { }
 
void
VTKPFEMXMLExportModule::initializeFrom(InputRecord &ir)
{
    ExportModule::initializeFrom(ir);
}
 
std::string
VTKPFEMXMLExportModule::giveOutputFileName(TimeStep *tStep)
{
    return this->giveOutputBaseFileName(tStep) + ".vtu";
}
 
 
std::ofstream
VTKPFEMXMLExportModule::giveOutputStream(TimeStep *tStep)
{
    std::string fileName = giveOutputFileName(tStep);
    std::ofstream streamF;
 
    if ( pythonExport ) {
        streamF = std::ofstream(NULL_DEVICE);//do not write anything
    } else {
        streamF = std::ofstream(fileName);
    }
 
    if ( !streamF.good() ) {
        OOFEM_ERROR("failed to open file %s", fileName.c_str() );
    }
 
    streamF.fill('0');//zero padding
    return streamF;
}
 
 
void
VTKPFEMXMLExportModule::doOutput(TimeStep *tStep, bool forcedOutput)
{
    if ( !( testTimeStepOutput(tStep) || forcedOutput ) ) {
        return;
    }
    
    this->fileStream = this->giveOutputStream(tStep);
    struct tm *current;
    time_t now;
    time(& now);
    current = localtime(& now);
 
    // Write output: VTK header
    this->fileStream << "<!-- TimeStep " << tStep->giveTargetTime() * timeScale << " Computed " << current->tm_year + 1900 << "-" << setw(2) << current->tm_mon + 1 << "-" << setw(2) << current->tm_mday << " at " << current->tm_hour << ":" << current->tm_min << ":" << setw(2) << current->tm_sec << " -->\n";
    this->fileStream << "<VTKFile type=\"UnstructuredGrid\" version=\"0.1\" byte_order=\"LittleEndian\">\n";
    this->fileStream << "<UnstructuredGrid>\n";
 
#ifdef __PFEM_MODULE
    // write out the particles (nodes exported as vertices = VTK_VERTEX)
    Domain *d  = emodel->giveDomain(1);
    int nnode = d->giveNumberOfDofManagers();
 
    int nActiveNode = 0;
    for ( int inode = 1; inode <= nnode; inode++ ) {
        PFEMParticle *particle = dynamic_cast< PFEMParticle * >( d->giveNode(inode) );
        if ( particle ) {
            if ( particle->isActive() ) {
                nActiveNode++;
            }
        }
    }
 
    DofManager *node;
    this->fileStream << "<Piece NumberOfPoints=\"" << nActiveNode << "\" NumberOfCells=\"" << nActiveNode << "\">\n";
    this->fileStream << "<Points>\n <DataArray type=\"Float64\" NumberOfComponents=\"3\" format=\"ascii\"> ";
 
    for ( int inode = 1; inode <= nnode; inode++ ) {
        node = d->giveNode(inode);
        PFEMParticle *particle = dynamic_cast< PFEMParticle * >( node );
        if ( particle ) {
            if ( particle->isActive() ) {
                const FloatArray &coords = node->giveCoordinates();
                for ( int i = 1; i <= coords.giveSize(); i++ ) {
                    this->fileStream << scientific << coords.at(i) << " ";
                }
 
                for ( int i = coords.giveSize() + 1; i <= 3; i++ ) {
                    this->fileStream << scientific << 0.0 << " ";
                }
            }
        }
    }
 
    this->fileStream << "</DataArray>\n</Points>\n";
 
 
    // output the cells connectivity data
    this->fileStream << "<Cells>\n";
    this->fileStream << " <DataArray type=\"Int32\" Name=\"connectivity\" format=\"ascii\"> ";
 
    for ( int ielem = 1; ielem <= nActiveNode; ielem++ ) {
        this->fileStream << ielem - 1 << " ";
    }
 
    this->fileStream << "</DataArray>\n";
 
    // output the offsets (index of individual element data in connectivity array)
    this->fileStream << " <DataArray type=\"Int32\" Name=\"offsets\" format=\"ascii\"> ";
 
    for ( int ielem = 1; ielem <= nActiveNode; ielem++ ) {
        this->fileStream << ielem << " ";
    }
    this->fileStream << "</DataArray>\n";
 
 
    // output cell (element) types
    this->fileStream << " <DataArray type=\"UInt8\" Name=\"types\" format=\"ascii\"> ";
    for ( int ielem = 1; ielem <= nActiveNode; ielem++ ) {
        this->fileStream << 1 << " ";
    }
 
    this->fileStream << "</DataArray>\n";
    this->fileStream << "</Cells>\n";
    this->fileStream << "</Piece>\n";
#endif //__PFEM_MODULE
 
    // Finalize the output:
    this->fileStream << "</UnstructuredGrid>\n</VTKFile>";
    if(this->fileStream){
        this->fileStream.close();
    }
 
 }
 
} // end namespace oofem