source: issm/trunk-jpl/src/wrappers/TriMesh/TriMesh.cpp@ 14221

/*
 * TriMesh: mesh a domain using an .exp file
 */

#include "./TriMesh.h"

void TriMeshUsage(void){/*{{{*/
        _printLine_("");
        _printLine_("   usage: [index,x,y,segments,segmentmarkers]=TriMesh(domainoutlinefilename,rifts,area) ");
        _printLine_("      where: index,x,y defines a triangulation, segments is an array made ");
        _printLine_("      of exterior segments to the mesh domain outline, segmentmarkers is an array flagging each segment, ");
        _printLine_("      outlinefilename an Argus domain outline file, ");
        _printLine_("      area is the maximum area desired for any element of the resulting mesh, ");
        _printLine_("");
}/*}}}*/
WRAPPER(TriMesh){
        
        /*intermediary: */
        double  area;
        DataSet *domain = NULL;
        DataSet *rifts  = NULL;

        /* output: */
        SeqMat<int>    *index             = NULL;
        SeqVec<double> *x                 = NULL;
        SeqVec<double> *y                 = NULL;
        SeqMat<int>    *segments          = NULL;
        SeqVec<int>    *segmentmarkerlist = NULL;

        /*Boot module: */
        MODULEBOOT();

        /*checks on arguments: */
        CHECKARGUMENTS(NLHS,NRHS,&TriMeshUsage);
        
        /*Fetch data needed for meshing: */
        FetchData(&domain,DOMAINOUTLINE);
        FetchData(&rifts,RIFTSOUTLINE);
        FetchData(&area,AREA);

        /*call x core: */
        TriMeshx(&index,&x,&y,&segments,&segmentmarkerlist,domain,rifts,area);

        /*write outputs: */
        WriteData(INDEX,index);
        WriteData(X,x);
        WriteData(Y,y);
        WriteData(SEGMENTS,segments);
        WriteData(SEGMENTMARKERLIST,segmentmarkerlist);

        /*free ressources: */
        delete domain;
        delete rifts;
        delete index;
        delete x;
        delete y;
        delete segments;
        delete segmentmarkerlist;

        /*end module: */
        MODULEEND();
}