inputOutput_example.cpp

Annotated source file

Here is the full file examples/inputOutput_example.cpp. Clicking on a function or class name will lead you to its reference documentation.

 
// For detailed options about visualization look at gsView.cpp example.
 
 
 
#include <iostream>
#include <gismo.h>
 
using namespace gismo;
 
 
int main(int argc, char* argv[])
{
 
    std::string input("curves3d/bspline3d_curve_01.xml");
    std::string output("");
    bool bezier=false; 
 
    gsCmdLine cmd("Tutorial Input Output");
    cmd.addPlainString("filename", "G+Smo input geometry file.", input);
    cmd.addString("o", "output", "Name of the output file", output);
    cmd.addSwitch("b","bezier", "Output using Bezier elements", bezier);
    try { cmd.getValues(argc,argv); } catch (int rv) { return rv; }
 
 
    std::string ext;
    ext = (bezier) ? ".vtu" : ".vtp";
 
    if (!gsFileManager::fileExists(input))
    {
        gsWarn << "The file cannot be found!\n";
        return EXIT_FAILURE;
    }
 
    gsInfo << "Read file \"" << input << "\"\n";
 
    gsFileData<> fileData(input);
 
    gsGeometry<>::uPtr pGeom;
    if (fileData.has< gsGeometry<> >())
    {
        pGeom = fileData.getFirst< gsGeometry<> >();
    }
    else
    {
        gsWarn << "Input file doesn't have a geometry.\n";
        return EXIT_FAILURE;
    }
 
    gsInfo << "The file contains: \n" << *pGeom << "\n";
 
    if ( output.empty() )
    {
        gsInfo << "Call program with option -o <basename> to write data to files\n";
        gsInfo << "<basename>Paraview"+ext+", <basename>Paraview.pvd, <basename>.xml, <basename>ControlPoints.csv\n";
        return EXIT_SUCCESS;
    }
 
    // writing a ParaView file
    const std::string out = output + "Paraview";
    if (!fileData.has< gsMultiPatch<> >()) 
    {
        if (!bezier)
        {
            gsWriteParaview(*pGeom, out);
        }
        else
        {
            gsMultiPatch<> mPatch;
            mPatch.addPatch(*pGeom);
            gsWriteParaviewBezier(mPatch, out);
        }
        gsInfo << "Wrote ParaView files: " << out << ext << ", " << out << ".pvd\n";
    }
    else
    {
        gsMultiPatch<> mPatch;
        fileData.getFirst(mPatch);
        if (!bezier)
        {
            gsWriteParaview(mPatch, out);
        }
        else
        {
            gsWriteParaviewBezier(mPatch, out);
        }
        gsInfo << "Wrote ParaView file: " << out << ".pvd\n";
    }
 
    // writing a G+Smo .xml file
    gsFileData<> fd;
    fd << *pGeom;
    // output is a string. The extention .xml is added automatically
    fd.save(output);
    gsInfo << "Wrote G+Smo file:     " << output << ".xml \n";
 
    // writing a .csv file
    if (pGeom->targetDim()==2)
        gsWriteCsv(output+".csv", pGeom->coefs(), {"X","Y"});
    else if (pGeom->targetDim()==3)
        gsWriteCsv(output+".csv", pGeom->coefs(), {"X","Y","Z"});
    gsInfo << "Wrote CSV file:       " << output+".csv \n";
 
    return EXIT_SUCCESS;
}