Skip to content

Commit

Permalink
Remove VHTMeshMorph for output because it doesn't really work now.
Browse files Browse the repository at this point in the history
As currently implemented, the coordinates stored on the mesh (not just
the ones managed by the FS) would have to change. I would like to make
the mesh only handle topology (at the file level, it is already that way
for analysis), in which case the coordinates could be written from a
bonified FS.
  • Loading branch information
jedbrown committed Jun 19, 2011
1 parent f0c3350 commit 0605df4
Showing 1 changed file with 0 additions and 20 deletions.
20 changes: 0 additions & 20 deletions src/fs/tests/vht.c
Original file line number Diff line number Diff line change
Expand Up @@ -553,24 +553,6 @@ static void VHTMeshMorph_NonDimensionalizeSI(void *vctx,double *coords)
coords[1] = scale * y;
coords[2] = scale * z;
}
static void VHTMeshMorph_NonDimensionalize(void *vctx,double *coords)
{
VHT vht = vctx;
dReal scale = dUnitNonDimensionalize(vht->scase->utable.Length,1.0); // nondimensionalize 1 meter
double x = coords[0],y = coords[1],z = coords[2];
coords[0] = scale * x;
coords[1] = scale * y;
coords[2] = scale * z;
}
static void VHTMeshMorph_Dimensionalize(void *vctx,double *coords)
{
VHT vht = vctx;
dReal scale = dUnitDimensionalize(vht->scase->utable.Length,1.0); // dimensionalize 1 length unit
double x = coords[0],y = coords[1],z = coords[2];
coords[0] = scale * x;
coords[1] = scale * y;
coords[2] = scale * z;
}
static dErr VHTSetFromOptions(VHT vht)
{
char scasename[256] = "Exact0";
Expand Down Expand Up @@ -2514,15 +2496,13 @@ int main(int argc,char *argv[])
err = PetscViewerFileSetName(view,dhmfilename);dCHK(err);
err = PetscViewerFileSetMode(view,FILE_MODE_WRITE);dCHK(err);
err = dFSGetMesh(vht->fsu,&mesh);dCHK(err);dCHK(err);
err = dMeshMorph(mesh,VHTMeshMorph_Dimensionalize,vht);dCHK(err);
err = VHTExtractGlobalSplit(vht,X,&Xu,&Xp,&Xe);dCHK(err);
err = dFSDirichletProject(vht->fsu,Xu,dFS_INHOMOGENEOUS);dCHK(err);
err = dFSDirichletProject(vht->fsp,Xp,dFS_INHOMOGENEOUS);dCHK(err);
err = dFSDirichletProject(vht->fse,Xe,dFS_INHOMOGENEOUS);dCHK(err);
err = VecView(Xu,view);dCHK(err);
err = VecView(Xp,view);dCHK(err);
err = VecView(Xe,view);dCHK(err);
err = dMeshMorph(mesh,VHTMeshMorph_NonDimensionalize,vht);dCHK(err);
err = PetscViewerDestroy(&view);dCHK(err);
}

Expand Down

0 comments on commit 0605df4

Please sign in to comment.