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

Last change on this file since 14656 was 14656, checked in by Eric.Larour, 12 years ago

CHG: implementation of revised ISSM toolkit.
We now have the following objects in the ISSM toolkit:
IssmMat and IssmVec: these are the wrappers to all our toolkit objects. These are hooked up to the
src/c/objects/matrix/Matrix.h and Vector.h objects.
We need of course enums that go with them, which map into Petsc constructs, such as MpiDenseEnum, DenseEnum,
etc ...
The toolkit now implements a MatDense matrix, and a future MatMpiDense matrix, as well as a SeqVec and
future MpiVec vector.
There is also an abstract class, called IssmAbsMat and IssmAbsVec, from which all our matrix and vector objects,
except for IssmMat and IssmVec, derive.
Updated all the wrappers and modules to use these new objects.
The toolkit options database is derived from the .toolkit file which is read at the beginning of any run. Very similar
to what Petsc does with its options database. Created a static class to hold this options database, in src/c/classes/ToolkitOptions.h
very similar to our static class holding the IssmComm.
File size: 1.4 KB
RevLine 
[1105]1/*! \file ContourToNodesx.c
2 */
3
4#include "./ContourToNodesx.h"
5
[14656]6int ContourToNodesx(IssmSeqVec<IssmPDouble>** pflags,double* x, double* y, int nods, Contour<IssmPDouble>** contours,int numcontours,int edgevalue){
[1105]7
8 int i;
9
10 /*Contour:*/
[13220]11 Contour<IssmPDouble>* contouri=NULL;
[8303]12 int numnodes;
[1105]13 double* xc=NULL;
14 double* yc=NULL;
15
16 /*output: */
[14656]17 IssmSeqVec<IssmPDouble>* flags=NULL;
18 flags=new IssmSeqVec<IssmPDouble>(nods);
[1105]19
20 /*Loop through all contours: */
21 for (i=0;i<numcontours;i++){
22 contouri=*(contours+i);
[8303]23 numnodes=contouri->nods;
[1105]24 xc=contouri->x;
25 yc=contouri->y;
[8303]26 IsInPoly(flags,xc,yc,numnodes,x,y,0,nods,edgevalue);
[1105]27 }
28
29 /*Assemble vector: */
[11695]30 flags->Assemble();
[1105]31
32 /*Assign output pointers: */
33 *pflags=flags;
[11401]34 return 1;
[1105]35}
[12121]36
[14656]37int ContourToNodesx(IssmSeqVec<IssmPDouble>** pflags,double* x, double* y, int nods, DataSet* contours, int edgevalue){
[12121]38
39 /*Contour:*/
[13220]40 Contour<IssmPDouble>* contouri=NULL;
[12121]41 double* xc=NULL;
42 double* yc=NULL;
43
44 /*output: */
[14656]45 IssmSeqVec<IssmPDouble>* flags=NULL;
46 flags=new IssmSeqVec<IssmPDouble>(nods);
[12121]47
48 /*Loop through all contours: */
49 if(contours){
[13762]50 for(int i=0;i<contours->Size();i++){
[13220]51 Contour<IssmPDouble>* contour=(Contour<IssmPDouble>*)contours->GetObjectByOffset(i);
[12121]52 IsInPoly(flags,contour->x,contour->y,contour->nods,x,y,0,nods,edgevalue);
53 }
54 }
55
56 /*Assemble vector: */
57 flags->Assemble();
58
59 /*Assign output pointers: */
60 *pflags=flags;
61 return 1;
62}
Note: See TracBrowser for help on using the repository browser.