source: issm/trunk/src/mex/TriMesh/TriMesh.cpp@ 11995

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

#include "./TriMesh.h"

WRAPPER(TriMesh){
        
        /* input: */
        char*  domainname=NULL;
        double area;
        bool   order;

        /*intermediary: */
        DataSet* domain=NULL;

        /* output: */
        Matrix* index=NULL;
        Vector* x=NULL;
        Vector* y=NULL;
        Matrix* segments=NULL;
        Vector* segmentmarkerlist=NULL;

        /*Boot module: */
        MODULEBOOT();

        /*checks on arguments on the matlab side: */
        CHECKARGUMENTS(NLHS,NRHS,&TriMeshUsage);
        
        /*Fetch data needed for meshing: */
        FetchData(&domainname,DOMAINOUTLINE);
        FetchData(&area,AREA);
        FetchData(&order,ORDER);


        /*Read domain outline: */
        domain=DomainOutlineRead(domainname,false);

        /*call x core: */
        TriMeshx(&index,&x,&y,&segments,&segmentmarkerlist,domain,area,order);
        
        
        /*write outputs: */
        WriteData(INDEX,index);
        WriteData(X,x);
        WriteData(Y,y);
        WriteData(SEGMENTS,segments);
        WriteData(SEGMENTMARKERLIST,segmentmarkerlist);

        /*free ressources: */
        delete domain;
        xdelete(&index);
        xdelete(&x);
        xdelete(&y);
        xdelete(&segments);
        xdelete(&segmentmarkerlist);

        /*end module: */
        MODULEEND();

}

void TriMeshUsage(void) //{{{1
{
        printf("\n");
        printf("   usage: [index,x,y,segments,segmentmarkers]=TriMesh(domainoutlinefilename,area,ordered) \n");
        printf("      where: index,x,y defines a triangulation, segments is an array made \n");
        printf("      of exterior segments to the mesh domain outline, segmentmarkers is an array flagging each segment, \n");
        printf("      outlinefilename an Argus domain outline file, \n");
        printf("      area is the maximum area desired for any element of the resulting mesh, \n");
        printf("      and ordered is a bool that determines whether segments are output in the \n");
        printf("      order they are made by Triangle (ie none), or ordered counter clockwise around the domain outline.\n");
        printf("\n");
}
//}}}