source: issm/trunk-jpl/src/c/modules/ContourToNodesx/ContourToNodesx.cpp@ 13787

/*! \file  ContourToNodesx.c
 */

#include "./ContourToNodesx.h"

int ContourToNodesx(SeqVec<IssmPDouble>** pflags,double* x, double* y, int nods, Contour<IssmPDouble>** contours,int numcontours,int edgevalue){

        int i;

        /*Contour:*/
        Contour<IssmPDouble>* contouri=NULL;
        int      numnodes;
        double*  xc=NULL;
        double*  yc=NULL;
        double   value;

        /*output: */
        SeqVec<IssmPDouble>* flags=NULL;
        flags=new SeqVec<IssmPDouble>(nods);

        /*Loop through all contours: */
        for (i=0;i<numcontours;i++){
                contouri=*(contours+i);
                numnodes=contouri->nods;
                xc=contouri->x;
                yc=contouri->y;
                IsInPoly(flags,xc,yc,numnodes,x,y,0,nods,edgevalue);
        }

        /*Assemble vector: */
        flags->Assemble();

        /*Assign output pointers: */
        *pflags=flags;

        return 1;
}

int ContourToNodesx(SeqVec<IssmPDouble>** pflags,double* x, double* y, int nods, DataSet* contours, int edgevalue){

        /*Contour:*/
        Contour<IssmPDouble>* contouri=NULL;
        double*  xc=NULL;
        double*  yc=NULL;

        /*output: */
        SeqVec<IssmPDouble>* flags=NULL;
        flags=new SeqVec<IssmPDouble>(nods);

        /*Loop through all contours: */
        if(contours){
                for(int i=0;i<contours->Size();i++){
                        Contour<IssmPDouble>* contour=(Contour<IssmPDouble>*)contours->GetObjectByOffset(i);
                        IsInPoly(flags,contour->x,contour->y,contour->nods,x,y,0,nods,edgevalue);
                }
        }

        /*Assemble vector: */
        flags->Assemble();

        /*Assign output pointers: */
        *pflags=flags;
        return 1;
}