vectorField solidZone2PointsDisplAtSolid =
stress().faceZone2PointDisplacementIncrement(solidZoneIndex());
solidZone2PointsDispl() =
// solidToFluid().pointInterpolate
ggiInterpolator2().slaveToMasterPointInterpolate
(
solidZone2PointsDisplAtSolid
);
residual2Prev() = residual2();
residual2() = solidZone2PointsDispl() - fluid2ZonePointsDispl();
// const scalarField& minLe = minEdgeLength();
// scalar residualNorm = gMax(mag(residual())/(minLe + SMALL));
scalar residual2Norm = ::sqrt(sum(magSqr(residual2())));
这些solidzone2我是想让它代表固体耦合边界二,但是这里一改,上面的成员函数、私有成员函数的相应的地方也要该,我怕一不小心改错了则麼办