Changeset 13762
- Timestamp:
- 10/19/12 17:31:29 (12 years ago)
- Location:
- issm/trunk-jpl/src/c/modules
- Files:
-
- 48 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/modules/AverageOntoPartitionx/AverageOntoPartitionx.cpp
r13622 r13762 20 20 void AverageOntoPartitionx(double** paverage, Elements* elements, Nodes* nodes, Vertices* vertices, Loads* loads, Materials* materials, Parameters* parameters,double* vertex_response){ 21 21 22 int i,j,k;23 22 int dummy; 24 25 23 int npart; 26 24 double *qmu_part = NULL; … … 48 46 49 47 /*loop on each element, and add contribution of the element to the partition (surface weighted average): */ 50 for(i =0;i<elements->Size();i++){48 for(int i=0;i<elements->Size();i++){ 51 49 Element* element=(Element*)elements->GetObjectByOffset(i); 52 50 element->AverageOntoPartition(partition_contributions,partition_areas,vertex_response,qmu_part); -
issm/trunk-jpl/src/c/modules/Bamgx/Bamgx.cpp
r13255 r13762 20 20 21 21 /*intermediary*/ 22 int i; 22 23 int noerr=1; 23 int i,j,num;24 24 double costheta=2; 25 25 double hminaniso=1e-100; -
issm/trunk-jpl/src/c/modules/ConstraintsStatex/RiftConstraintsState.cpp
r13622 r13762 45 45 int num_unstable_constraints=0; 46 46 int converged=0; 47 int potential;48 47 49 48 RiftConstrain(&num_unstable_constraints,loads,configuration_type); -
issm/trunk-jpl/src/c/modules/ContourToMeshx/ContourToMeshx.cpp
r13622 r13762 11 11 12 12 int ContourToMeshx(SeqVec<double>** pin_nod,SeqVec<double>** pin_elem, double* index, double* x, double* y,DataSet* contours,char* interptype,int nel,int nods, int edgevalue) { 13 14 int noerr=1;15 int i;16 int m,n;17 13 18 14 /*Contour:*/ … … 52 48 /*Take care of the case where an element interpolation has been requested: */ 53 49 if ((strcmp(interptype,"element")==0) || (strcmp(interptype,"element and node")==0)){ 54 for (n=0;n<nel;n++){50 for(int n=0;n<nel;n++){ 55 51 if ( (in_nod_serial[ (int)*(index+3*n+0) -1] == 1) && (in_nod_serial[ (int)*(index+3*n+1) -1] == 1) && (in_nod_serial[ (int)*(index+3*n+2) -1] == 1) ){ 56 52 value=1; in_elem->SetValue(n,value,INS_VAL); … … 69 65 xDelete<double>(in_nod_serial); 70 66 71 return noerr;67 return 1; 72 68 } -
issm/trunk-jpl/src/c/modules/ContourToNodesx/ContourToNodesx.cpp
r13622 r13762 40 40 int ContourToNodesx(SeqVec<IssmPDouble>** pflags,double* x, double* y, int nods, DataSet* contours, int edgevalue){ 41 41 42 int i;43 int m,n;44 45 42 /*Contour:*/ 46 43 Contour<IssmPDouble>* contouri=NULL; … … 56 53 /*Loop through all contours: */ 57 54 if(contours){ 58 for (i=0;i<contours->Size();i++){55 for(int i=0;i<contours->Size();i++){ 59 56 Contour<IssmPDouble>* contour=(Contour<IssmPDouble>*)contours->GetObjectByOffset(i); 60 57 IsInPoly(flags,contour->x,contour->y,contour->nods,x,y,0,nods,edgevalue); -
issm/trunk-jpl/src/c/modules/HoleFillerx/HoleFillerx.cpp
r13622 r13762 11 11 int HoleFillerx(double** pimageout,double* image, int lines,int samps,int smooth){ 12 12 13 FILE *fp1;14 unsigned long filesize;15 13 long infinit; 16 14 double temp; 17 long i, j, k, l, ii, jj , bytes;15 long i, j, k, l, ii, jj; 18 16 long iii, jjj; 19 long test, prevlines, prevpix, test1, test2, thresh,counter; 20 long bell, cycles, leftovers; 17 long test; 21 18 float howlong; 22 19 float nsteps, ssteps, wsteps, esteps; … … 38 35 double* image3=NULL; 39 36 double* image4=NULL; 40 41 /*Added to original routine: */42 int imageoutsize;43 37 44 38 /*^^^^^^^^^^^^^ Remove pixels close to the holes ^^^^^^^^^^^^^*/ -
issm/trunk-jpl/src/c/modules/InterpFromMeshToMesh2dx/InterpFromMeshToMesh2dx.cpp
r13622 r13762 30 30 double areacoord[3]; 31 31 double aa,bb; 32 double data_value;33 32 Icoor2 dete[3]; 34 33 -
issm/trunk-jpl/src/c/modules/IoModelToConstraintsx/IoModelToConstraintsx.cpp
r13622 r13762 14 14 /*intermediary: */ 15 15 int i,j; 16 bool transient = false; 17 18 FILE *fid = NULL; 19 int code = 0; 20 int vector_layout = 0; 21 int counter; 22 int nods; 23 IssmDouble* times=NULL; 24 IssmDouble* values=NULL; 25 bool spcpresent=false; 26 int count=0; 27 int numberofvertices; 16 bool transient = false; 17 FILE *fid = NULL; 18 int code = 0; 19 int vector_layout = 0; 20 IssmDouble *times = NULL; 21 IssmDouble *values = NULL; 22 bool spcpresent = false; 23 int count = 0; 24 int numberofvertices; 28 25 29 26 /*variables being fetched: */ -
issm/trunk-jpl/src/c/modules/Kml2Expx/Kml2Expx.cpp
r13622 r13762 25 25 int sgn,double cm,double sp){ 26 26 27 int i ,iret=0;27 int iret=0; 28 28 double *lat=NULL,*lon=NULL; 29 29 -
issm/trunk-jpl/src/c/modules/MeshProfileIntersectionx/ElementSegment.cpp
r13622 r13762 9 9 * If they do, create a Segment object with the intersection, and add to segments_dataset dataset: */ 10 10 11 int i;12 double alpha;13 11 double alpha1,alpha2; 14 12 double beta1,beta2; -
issm/trunk-jpl/src/c/modules/MeshProfileIntersectionx/MeshProfileIntersectionx.cpp
r13622 r13762 7 7 8 8 int i,j,k; 9 int m,n;10 9 11 10 /*Contour:*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateConstraintsBalancethickness.cpp
r13622 r13762 10 10 void CreateConstraintsBalancethickness(Constraints** pconstraints, IoModel* iomodel){ 11 11 12 /*Fetch parameters: */ 12 13 int stabilization; 13 14 /*Fetch parameters: */15 14 iomodel->Constant(&stabilization,BalancethicknessStabilizationEnum); 16 15 17 /*Output*/18 Constraints* constraints = NULL;19 20 16 /*Recover pointer: */ 21 constraints=*pconstraints;17 Constraints* constraints=*pconstraints; 22 18 23 19 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateLoadsBalancethickness.cpp
r13622 r13762 14 14 15 15 /*Intermediary*/ 16 int i;17 16 int element; 18 17 int stabilization; … … 23 22 iomodel->Constant(&numberofedges,MeshNumberofedgesEnum); 24 23 25 /*Output*/26 Loads* loads = NULL;27 28 24 /*Recover pointer: */ 29 loads=*ploads;25 Loads* loads=*ploads; 30 26 31 27 /*Create loads if they do not exist yet*/ … … 39 35 40 36 /*First load data:*/ 41 for (i =0;i<numberofedges;i++){37 for (int i=0;i<numberofedges;i++){ 42 38 43 39 /*Get left and right elements*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Balancethickness/CreateNodesBalancethickness.cpp
r13073 r13762 26 26 int stabilization; 27 27 28 /*DataSets: */29 Nodes* nodes = NULL;30 31 28 /*Fetch parameters: */ 32 29 iomodel->Constant(&dim,MeshDimensionEnum); … … 36 33 37 34 /*Recover pointer: */ 38 nodes=*pnodes;35 Nodes* nodes=*pnodes; 39 36 40 37 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateConstraintsBedSlope.cpp
r13622 r13762 13 13 void CreateConstraintsBedSlope(Constraints** pconstraints, IoModel* iomodel){ 14 14 15 /*Output*/16 Constraints* constraints = NULL;17 18 15 /*Recover pointer: */ 19 constraints=*pconstraints;16 Constraints* constraints=*pconstraints; 20 17 21 18 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateLoadsBedSlope.cpp
r12832 r13762 13 13 void CreateLoadsBedSlope(Loads** ploads, IoModel* iomodel){ 14 14 15 /*DataSet*/16 Loads* loads = NULL;17 18 15 /*Recover pointer: */ 19 loads=*ploads;16 Loads* loads=*ploads; 20 17 21 18 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/BedSlope/CreateNodesBedSlope.cpp
r13622 r13762 18 18 int i; 19 19 bool continuous_galerkin=true; 20 int numberofvertices;21 22 /*DataSets: */23 Nodes* nodes = NULL;24 20 25 21 /*Fetch parameters: */ 22 int numberofvertices; 26 23 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); 27 24 28 25 /*Recover pointer: */ 29 nodes=*pnodes;26 Nodes* nodes=*pnodes; 30 27 31 28 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateLoadsDiagnosticHoriz.cpp
r13622 r13762 14 14 15 15 /*DataSets*/ 16 Loads *loads = NULL;17 16 Icefront *icefront = NULL; 18 17 Riftfront *riftfront = NULL; … … 48 47 49 48 /*Recover pointer: */ 50 loads=*ploads;49 Loads* loads=*ploads; 51 50 52 51 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateNodesDiagnosticHoriz.cpp
r13622 r13762 21 21 bool isstokes,isl1l2,ismacayealpattyn; 22 22 23 /*DataSets: */24 Nodes* nodes = NULL;25 26 23 /*Fetch parameters: */ 27 24 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); … … 31 28 32 29 /*Recover pointer: */ 33 nodes=*pnodes;30 Nodes* nodes=*pnodes; 34 31 35 32 /*First create nodes*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateLoadsDiagnosticHutter.cpp
r13622 r13762 13 13 void CreateLoadsDiagnosticHutter(Loads** ploads, IoModel* iomodel){ 14 14 15 /*DataSet*/16 Loads* loads=NULL;17 18 15 /*Recover pointer: */ 19 loads=*ploads;16 Loads* loads=*ploads; 20 17 21 18 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHutter/CreateNodesDiagnosticHutter.cpp
r12832 r13762 21 21 bool ishutter; 22 22 23 /*DataSets: */24 Nodes* nodes = NULL;25 26 23 /*Fetch parameters: */ 27 24 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); … … 29 26 30 27 /*Recover pointer: */ 31 nodes=*pnodes;28 Nodes* nodes=*pnodes; 32 29 33 30 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateConstraintsDiagnosticVert.cpp
r13073 r13762 20 20 int numberofvertices; 21 21 22 /*Output*/23 Constraints* constraints = NULL;24 25 22 /*Fetch parameters: */ 26 23 iomodel->Constant(&dim,MeshDimensionEnum); … … 29 26 30 27 /*Recover pointer: */ 31 constraints=*pconstraints;28 Constraints* constraints=*pconstraints; 32 29 33 30 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateLoadsDiagnosticVert.cpp
r12832 r13762 13 13 void CreateLoadsDiagnosticVert(Loads** ploads, IoModel* iomodel){ 14 14 15 /*DataSet*/16 Loads* loads=NULL;17 18 15 /*Recover pointer: */ 19 loads=*ploads;16 Loads* loads=*ploads; 20 17 21 18 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticVert/CreateNodesDiagnosticVert.cpp
r13622 r13762 21 21 int numberofvertices; 22 22 23 /*DataSets: */24 Nodes* nodes = NULL;25 26 23 /*Fetch parameters: */ 27 24 iomodel->Constant(&dim,MeshDimensionEnum); … … 29 26 30 27 /*Recover pointer: */ 31 nodes=*pnodes;28 Nodes* nodes=*pnodes; 32 29 33 30 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateConstraintsEnthalpy.cpp
r13622 r13762 28 28 IssmDouble* times=NULL; 29 29 IssmDouble* values=NULL; 30 Constraints* constraints = NULL;31 30 32 31 /*Fetch parameters: */ … … 37 36 38 37 /*Recover pointer: */ 39 constraints=*pconstraints;38 Constraints* constraints=*pconstraints; 40 39 41 40 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateLoadsEnthalpy.cpp
r13622 r13762 13 13 void CreateLoadsEnthalpy(Loads** ploads, IoModel* iomodel){ 14 14 15 /*DataSet*/16 Loads* loads=NULL;17 18 15 /*Recover pointer: */ 19 loads=*ploads;16 Loads* loads=*ploads; 20 17 21 18 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Enthalpy/CreateNodesEnthalpy.cpp
r13622 r13762 20 20 int numberofvertices; 21 21 22 /*DataSets: */23 Nodes* nodes = NULL;24 25 22 /*Fetch parameters: */ 26 23 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); 27 24 28 25 /*Recover pointer: */ 29 nodes=*pnodes;26 Nodes* nodes=*pnodes; 30 27 31 28 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Hydrology/CreateConstraintsHydrology.cpp
r13622 r13762 14 14 void CreateConstraintsHydrology(Constraints** pconstraints, IoModel* iomodel){ 15 15 16 /*Output*/17 Constraints *constraints = NULL;18 19 16 /*Recover pointer: */ 20 constraints=*pconstraints;17 Constraints* constraints=*pconstraints; 21 18 22 19 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Hydrology/CreateLoadsHydrology.cpp
r12832 r13762 13 13 void CreateLoadsHydrology(Loads** ploads, IoModel* iomodel){ 14 14 15 /*Intermediary*/16 int i;17 18 /*DataSet*/19 Loads* loads = NULL;20 21 15 /*Recover pointer: */ 22 loads=*ploads;16 Loads* loads=*ploads; 23 17 24 18 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Hydrology/CreateNodesHydrology.cpp
r13622 r13762 20 20 int numberofvertices; 21 21 22 /*DataSets: */23 Nodes* nodes = NULL;24 25 22 /*Fetch parameters: */ 26 23 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); 27 24 28 25 /*Recover pointer: */ 29 nodes=*pnodes;26 Nodes* nodes=*pnodes; 30 27 31 28 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateConstraintsMelting.cpp
r13622 r13762 13 13 void CreateConstraintsMelting(Constraints** pconstraints, IoModel* iomodel){ 14 14 15 /*Intermediary*/16 int i;17 int count;18 19 /*Intermediary*/20 Constraints* constraints = NULL;21 22 15 /*Recover pointer: */ 23 constraints=*pconstraints;16 Constraints* constraints=*pconstraints; 24 17 25 18 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateLoadsMelting.cpp
r13073 r13762 14 14 15 15 /*Intermediary*/ 16 int i; 17 int dim; 18 int numberofvertices; 19 20 /*DataSet*/ 21 Loads* loads = NULL; 22 23 /*Fetch parameters: */ 16 int dim; 17 int numberofvertices; 24 18 iomodel->Constant(&dim,MeshDimensionEnum); 25 19 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); … … 29 23 30 24 /*Recover pointer: */ 31 loads=*ploads;25 Loads* loads=*ploads; 32 26 33 27 /*Create loads if they do not exist yet*/ … … 38 32 CreateSingleNodeToElementConnectivity(iomodel); 39 33 40 for (i=0;i<numberofvertices;i++){34 for(int i=0;i<numberofvertices;i++){ 41 35 if((iomodel->my_vertices[i]==1)){ 42 36 if (reCast<int>(iomodel->Data(MeshVertexonbedEnum)[i])){ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Melting/CreateNodesMelting.cpp
r13622 r13762 20 20 int numberofvertices; 21 21 22 /*DataSets: */23 Nodes* nodes = NULL;24 25 22 /*Fetch parameters: */ 26 23 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); 27 24 28 25 /*Recover pointer: */ 29 nodes=*pnodes;26 Nodes* nodes=*pnodes; 30 27 31 28 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateConstraintsPrognostic.cpp
r13622 r13762 10 10 void CreateConstraintsPrognostic(Constraints** pconstraints, IoModel* iomodel){ 11 11 12 /*Fetch parameters: */ 12 13 int stabilization; 13 14 /*Fetch parameters: */15 14 iomodel->Constant(&stabilization,PrognosticStabilizationEnum); 16 15 17 /*Output*/18 Constraints *constraints = NULL;19 20 16 /*Recover pointer: */ 21 constraints=*pconstraints;17 Constraints* constraints=*pconstraints; 22 18 23 19 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateLoadsPrognostic.cpp
r13622 r13762 23 23 int numberofelements; 24 24 25 /*DataSet*/26 Loads* loads = NULL;27 28 25 /*Fetch parameters: */ 29 26 iomodel->Constant(&stabilization,PrognosticStabilizationEnum); … … 32 29 33 30 /*Recover pointer: */ 34 loads=*ploads;31 Loads* loads=*ploads; 35 32 36 33 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Prognostic/CreateNodesPrognostic.cpp
r13056 r13762 26 26 int stabilization; 27 27 28 /*DataSets: */29 Nodes* nodes = NULL;30 31 28 /*Fetch parameters: */ 32 29 iomodel->Constant(&dim,MeshDimensionEnum); … … 36 33 37 34 /*Recover pointer: */ 38 nodes=*pnodes;35 Nodes* nodes=*pnodes; 39 36 40 37 /*Create nodes if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateConstraintsSurfaceSlope.cpp
r13622 r13762 13 13 void CreateConstraintsSurfaceSlope(Constraints** pconstraints, IoModel* iomodel){ 14 14 15 /*Output*/16 Constraints* constraints = NULL;17 18 15 /*Recover pointer: */ 19 constraints=*pconstraints;16 Constraints* constraints=*pconstraints; 20 17 21 18 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateLoadsSurfaceSlope.cpp
r12832 r13762 13 13 void CreateLoadsSurfaceSlope(Loads** ploads, IoModel* iomodel){ 14 14 15 /*DataSet*/16 Loads* loads = NULL;17 18 15 /*Recover pointer: */ 19 loads=*ploads;16 Loads* loads=*ploads; 20 17 21 18 /*Create loads if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/SurfaceSlope/CreateNodesSurfaceSlope.cpp
r13622 r13762 16 16 17 17 /*Intermediary*/ 18 int i;19 18 bool continuous_galerkin=true; 20 19 int numberofvertices; 21 22 /*DataSets: */23 Nodes* nodes = NULL;24 20 25 21 /*Fetch parameters: */ … … 27 23 28 24 /*Recover pointer: */ 29 nodes=*pnodes;25 Nodes* nodes=*pnodes; 30 26 31 27 /*Create nodes if they do not exist yet*/ … … 38 34 iomodel->FetchData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,MaskVertexonwaterEnum); 39 35 40 for (i=0;i<numberofvertices;i++){ 41 36 for(int i=0;i<numberofvertices;i++){ 42 37 if(iomodel->my_vertices[i]){ 43 38 -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateConstraintsThermal.cpp
r13622 r13762 14 14 void CreateConstraintsThermal(Constraints** pconstraints, IoModel* iomodel){ 15 15 16 /*Intermediary*/ 17 int i; 18 int count; 16 /*Fetch parameters: */ 19 17 int dim; 20 21 /*Output*/22 Constraints* constraints = NULL;23 24 /*Fetch parameters: */25 18 iomodel->Constant(&dim,MeshDimensionEnum); 26 19 27 20 /*Recover pointer: */ 28 constraints=*pconstraints;21 Constraints* constraints=*pconstraints; 29 22 30 23 /*Create constraints if they do not exist yet*/ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateLoadsThermal.cpp
r13622 r13762 17 17 int dim; 18 18 int numberofvertices; 19 20 /*DataSet*/21 Loads* loads = NULL;22 19 Pengrid* pengrid = NULL; 23 20 24 21 /*Recover pointer: */ 25 loads=*ploads;22 Loads* loads=*ploads; 26 23 27 24 /*Fetch parameters: */ -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Thermal/CreateNodesThermal.cpp
r13622 r13762 16 16 17 17 /*Intermediary*/ 18 int i;19 18 bool continuous_galerkin=true; 20 int numberofvertices;21 22 /*DataSets: */23 Nodes* nodes = NULL;24 19 25 20 /*Fetch parameters: */ 21 int numberofvertices; 26 22 iomodel->Constant(&numberofvertices,MeshNumberofverticesEnum); 27 23 28 24 /*Recover pointer: */ 29 nodes=*pnodes;25 Nodes* nodes=*pnodes; 30 26 31 27 /*Create nodes if they do not exist yet*/ … … 38 34 iomodel->FetchData(6,MeshVertexonbedEnum,MeshVertexonsurfaceEnum,MaskVertexongroundediceEnum,MaskVertexonfloatingiceEnum,FlowequationVertexEquationEnum,MaskVertexonwaterEnum); 39 35 40 for (i=0;i<numberofvertices;i++){ 41 36 for(int i=0;i<numberofvertices;i++){ 42 37 if(iomodel->my_vertices[i]){ 43 38 /*Add node to nodes dataset: */ -
issm/trunk-jpl/src/c/modules/NodesDofx/NodesDofx.cpp
r13622 r13762 14 14 int noerr=1; 15 15 int found=0; 16 int i;17 16 18 17 /*Do we have any nodes for this analysis type? :*/ -
issm/trunk-jpl/src/c/modules/Orthx/Orthx.cpp
r13622 r13762 11 11 12 12 /*intermediary:*/ 13 IssmDouble norm_ new,norm_old,dot_product;;13 IssmDouble norm_old,dot_product;; 14 14 15 15 /*Initialize output*/ -
issm/trunk-jpl/src/c/modules/Reduceloadx/Reduceloadx.cpp
r13622 r13762 20 20 int global_m,global_n; 21 21 bool fromlocalsize = true; 22 int verbose;23 22 24 23 if(VerboseModule()) _pprintLine_(" Dirichlet lifting applied to load vector"); -
issm/trunk-jpl/src/c/modules/Solverx/Solverx.cpp
r13216 r13762 15 15 16 16 void Solverx(Vector<IssmDouble>** puf, Matrix<IssmDouble>* Kff, Vector<IssmDouble>* pf, Vector<IssmDouble>* uf0,Vector<IssmDouble>* df, Parameters* parameters){ 17 18 /*Intermediary: */19 int analysis_type;20 17 21 18 /*output: */ -
issm/trunk-jpl/src/c/modules/Solverx/SolverxPetsc.cpp
r13622 r13762 36 36 /*Intermediary: */ 37 37 int local_m,local_n,global_m,global_n; 38 int analysis_type;39 38 40 39 /*Solver */ -
issm/trunk-jpl/src/c/modules/TriaSearchx/TriaSearchx.cpp
r13622 r13762 24 24 R2 r; 25 25 I2 I; 26 int i,j,k; 27 int i0,i1,i2; 28 double areacoord[3]; 29 double aa,bb; 30 double data_value; 26 int i; 31 27 Icoor2 dete[3]; 32 int verbose=0;33 28 34 / / read background mesh29 /* read background mesh */ 35 30 Mesh Th(index,x,y,nods,nel); 36 31 Th.CreateSingleVertexToTriangleConnectivity();
Note:
See TracChangeset
for help on using the changeset viewer.