[bullmut]

Good Day

I am trying to follow the code for a surfaceDiscplament motion.
To that end i am following the print out call for a single time step from the solver moveDynamicMesh. I am finally at the point where i tracing the steps of the function call

mesh.update();

line 136 in moveDynamicMesh.C
Please bear with me...
This calls the dynamicMotionSolverFvMesh base class which inherits from the dynamicFvMesh class. The function update in this class actually calls the fvMesh class function movePoints

bool Foam::dynamicMotionSolverFvMesh::update()
{
    fvMesh::movePoints(motionPtr_->newPoints());

    if (foundObject<volVectorField>("U"))
    {
        volVectorField& U =
            const_cast<volVectorField&>(lookupObject<volVectorField>("U"));
        U.correctBoundaryConditions();
    }

    return true;
}

But before the fvMesh function can be implemented we need to motionSolver pointer to execute the function newPoints.The motionSolver class is a virtual base class that is used to build specific forms of motion. In this case i am looking at displacementMotionsolver. So now i am looking at the base class motionSolver and the virtual function newPoints exists:

Foam::tmp<Foam::pointField> Foam::motionSolver::newPoints()
{
    solve();
    return curPoints();
}

Now it is not directly replaced by a displacementMotionsolver function however there is a solve() function in the class that inherits from displacementMotionsolver, namely displacementLaplacianFvMotionsolver. Which is ultimately the motion solver i am implementing.
My question is do i assume the final motion solver (displacementLaplacianFvMotionsolver) is being called to execute solve() in which case the code reads:

void Foam::displacementLaplacianFvMotionSolver::solve()
{
    // The points have moved so before interpolation update
    // the motionSolver accordingly
    movePoints(fvMesh_.points());

    diffusivity().correct();
    pointDisplacement_.boundaryField().updateCoeffs();

    Foam::solve
    (
        fvm::laplacian
        (
            diffusivity().operator()(),
            cellDisplacement_,
            "laplacian(diffusivity,cellDisplacement)"
        )
    );
}

which makes little sense to me or do i assume it is the solve function as defined by the base virtual class motionSolver as below:

//- Solve for motion
        virtual void solve() = 0;

In which case the code skips the 3rd code block i posted would just jump straight to the function curPoints.
Any help is greatly appreciated and thank you for even reading this lengthy post