Changeset 12004 for issm/branches/trunk-jpl-damage/src/c/shared/Exp/DomainOutlineRead.cpp

Timestamp:

04/16/12 16:06:52 (13 years ago)

Author:

cborstad

Message:

merged trunk-jpl into trunk-jpl-damage through revision 11990

Location:

issm/branches/trunk-jpl-damage

Files:

• . (modified) (1 prop)
• src/c/shared/Exp/DomainOutlineRead.cpp (modified) (5 diffs)

issm/branches/trunk-jpl-damage/src/c/shared/Exp/DomainOutlineRead.cpp

@@ -9 +9 @@
 #include "../Alloc/alloc.h"
 #include "../../include/include.h"
+#include "../../objects/objects.h"
 #include "../Exceptions/exceptions.h"
+#include "../../Container/DataSet.h"
 
 int DomainOutlineRead(int* pnprof,int** pprofnvertices,double*** ppprofx,double*** ppprofy,bool** pclosed,char* domainname,bool whole=true){
 
         
-        /*Error management: */
-        int noerr=1;
+        /*indexing: */
         int i,counter;
 
@@ -39 +40 @@
         if ((fid=fopen(domainname,"r"))==NULL){
                 _error_("%s%s","could not find domain file ",domainname); 
-                noerr=0; goto cleanupandreturn;
         }
 
@@ -95 +95 @@
                 /*Read vertices: */
                 for (i=0;i<n;i++){
-                        fscanf(fid,"%20lf %20lf\n",x+i,y+i);
+                        fscanf(fid,"%20lf %20lf\n",&x[i],&y[i]);
                 }
 
@@ -117 +117 @@
         fclose(fid);
 
-        cleanupandreturn: 
-        /*Free ressources: */
-        if(!noerr){
-                xfree((void**)&x);
-                xfree((void**)&y);
-                n=-1;
-        }
-        
         /*Assign output pointers: */
         *pnprof=nprof;
@@ -132 +124 @@
         if(pclosed)*pclosed=closed;
         else       xfree((void**)&closed);
-        return noerr;
 }
 
+DataSet* DomainOutlineRead(char* domainname,bool whole=true){
+
+        /*indexing: */
+        int i;
+
+        /*intermediary: */
+        int nprof;
+        int* profnvertices=NULL;
+        double** pprofx=NULL;
+        double** pprofy=NULL;
+
+        Contour* contour=NULL;
+
+        /*output: */
+        DataSet* domain=NULL;
+
+        /*get domain outline from intermediary function:*/
+        DomainOutlineRead(&nprof,&profnvertices,&pprofx, &pprofy, NULL,domainname,whole);
+
+        /*now create dataset of contours: */
+        domain=new DataSet(0);
+
+        for(i=0;i<nprof;i++){
+                domain->AddObject(new Contour(i,profnvertices[i],pprofx[i],pprofy[i],1));
+        }
+
+        return domain;
+}