[cig-commits] commit: Make staticTopFront point to the right index
Mercurial
hg at geodynamics.org
Tue Apr 17 14:12:58 PDT 2012
changeset: 927:b767b5fbf956
tag: tip
user: Walter Landry <wlandry at caltech.edu>
date: Tue Apr 17 14:12:50 2012 -0700
files: plugins/EulerDeform/EulerDeform_Build.cxx plugins/EulerDeform/EulerDeform_Remesh.cxx
description:
Make staticTopFront point to the right index
diff -r 9fdb89ec9676 -r b767b5fbf956 plugins/EulerDeform/EulerDeform_Build.cxx
--- a/plugins/EulerDeform/EulerDeform_Build.cxx Tue Apr 17 12:31:38 2012 -0700
+++ b/plugins/EulerDeform/EulerDeform_Build.cxx Tue Apr 17 14:12:50 2012 -0700
@@ -316,14 +316,14 @@ void EulerDeform_Build(void* component,
sys->static_left_coord =
Dictionary_GetDouble( uwCtx->dictionary, "minX");
+ if(sys->staticBack && sys->wrapTop
+ && !sys->staticTopBack)
+ sys->static_back_coord =
+ Dictionary_GetDouble( uwCtx->dictionary, "minZ");
+
if(sys->staticFront && sys->wrapTop
&& !sys->staticTopFront)
sys->static_front_coord =
- Dictionary_GetDouble( uwCtx->dictionary, "minZ");
-
- if(sys->staticBack && sys->wrapTop
- && !sys->staticTopBack)
- sys->static_back_coord =
Dictionary_GetDouble( uwCtx->dictionary, "maxZ");
}
}
diff -r 9fdb89ec9676 -r b767b5fbf956 plugins/EulerDeform/EulerDeform_Remesh.cxx
--- a/plugins/EulerDeform/EulerDeform_Remesh.cxx Tue Apr 17 12:31:38 2012 -0700
+++ b/plugins/EulerDeform/EulerDeform_Remesh.cxx Tue Apr 17 14:12:50 2012 -0700
@@ -107,17 +107,8 @@ void EulerDeform_Remesh(TimeIntegrand* c
/* Copy back coords. */
for( ind_i = 0; ind_i < nInds; ind_i++ )
- {
- std::cout << "side "
- << ind_i << " "
- << inds[ind_i] << " "
- << sys->sideCoords[ind_i][0] << " "
- << sys->sideCoords[ind_i][1] << " "
- << sys->sideCoords[ind_i][2] << " "
- << "\n";
memcpy( sys->mesh->verts[inds[ind_i]], sys->sideCoords[ind_i],
nDims * sizeof(double));
- }
FreeObject( tmpIndSet );
FreeArray( sys->sideCoords );
@@ -136,12 +127,12 @@ void EulerDeform_Remesh(TimeIntegrand* c
}
if(sys->staticFront && !sys->staticTopFront)
{
- EulerDeform_Remesh_Corner(sys->mesh,0,1,sys->static_front_coord,2,1,0);
+ EulerDeform_Remesh_Corner(sys->mesh,grid->sizes[2]-1,
+ grid->sizes[2]-2,sys->static_front_coord,2,1,0);
}
if(sys->staticBack && !sys->staticTopBack)
{
- EulerDeform_Remesh_Corner(sys->mesh,grid->sizes[2]-1,
- grid->sizes[2]-2,sys->static_back_coord,2,1,0);
+ EulerDeform_Remesh_Corner(sys->mesh,0,1,sys->static_back_coord,2,1,0);
}
}
@@ -196,25 +187,9 @@ void EulerDeform_Remesh(TimeIntegrand* c
for( n_i = 0; n_i < nDomainNodes; n_i++ )
memcpy( oldCrds[n_i], sys->mesh->verts[n_i], nDims * sizeof(double) );
- for( n_i = 0; n_i < nDomainNodes; n_i++ )
- std::cout << "coord "
- << n_i << " "
- << sys->mesh->verts[n_i][0] << " "
- << sys->mesh->verts[n_i][1] << " "
- << sys->mesh->verts[n_i][2] << " "
- << "\n";
-
/* Remesh the system. */
Stg_Component_Execute( sys->remesher, NULL, True );
Mesh_Sync( sys->mesh );
-
- for( n_i = 0; n_i < nDomainNodes; n_i++ )
- std::cout << "Before "
- << n_i << " "
- << sys->mesh->verts[n_i][0] << " "
- << sys->mesh->verts[n_i][1] << " "
- << sys->mesh->verts[n_i][2] << " "
- << "\n";
/* Swap old coordinates back in temporarily. */
newCrds = sys->mesh->verts;
More information about the CIG-COMMITS
mailing list