1 //--------------------------------------------------//
3 //--------------------------------------------------//
4 // if(min(J.internalField()) > 0)
6 Info << "Move solid mesh using inverse distance interpolation" << endl;
11 // Create point interpolation
12 volPointInterpolation pointInterpolation(mesh);
16 pMesh.boundary().size(),
17 calculatedFvPatchVectorField::typeName
20 pointVectorField pointDU
29 dimensionedVector("zero", dimLength, vector::zero),
33 // Calculate mesh points displacement
34 pointInterpolation.interpolate(DU, pointDU);
36 //- correct edge interpolation
37 //- this is the stuff from edgeCorrectedVolPointInterpolation but
38 //- that class no longer works
39 # include "performEdgeCorrectedVolPointInterpolation.H"
41 const vectorField& pointDUI =
42 pointDU.internalField();
44 //- see the effect of correctBCs
47 vectorField newPoints = mesh.allPoints();
49 forAll (pointDUI, pointI)
51 newPoints[pointI] += pointDUI[pointI];
54 twoDPointCorrector twoDCorrector(mesh);
55 twoDCorrector.correctPoints(newPoints);
56 mesh.movePoints(newPoints);
62 // FatalErrorIn(args.executable())
63 // << "Negative Jacobian"
64 // << exit(FatalError);