void NeoHookeanMaterialNonLinear<Mesh>::updateJacobianMatrix ( const vector_Type&       disp,
                                                               const dataPtr_Type&      dataMaterial,
                                                               const mapMarkerVolumesPtr_Type mapsMarkerVolumes,
                                                               const displayerPtr_Type& displayer )
{
    this->M_jacobian.reset (new matrix_Type (*this->M_localMap) );

    displayer->leaderPrint (" \n*********************************\n  ");
    updateNonLinearJacobianTerms (this->M_jacobian, disp, dataMaterial, mapsMarkerVolumes, displayer);
    displayer->leaderPrint (" \n*********************************\n  ");
    std::cout << std::endl;

    this->M_jacobian->globalAssemble();
}
void VenantKirchhoffMaterialLinear<Mesh>::updateJacobianMatrix(const vector_Type& disp,
                                                               const dataPtr_Type& dataMaterial,
                                                               const mapMarkerVolumesPtr_Type mapsMarkerVolumes,
                                                               const displayerPtr_Type& displayer)
{
    //displayer->leaderPrint(" \n*********************************\n  ");
    displayer->leaderPrint("  S-  Updating the Jacobian Matrix (constant, Linear Elastic)\n");
    //displayer->leaderPrint(" \n*********************************\n  ");

    //displayer->leaderPrint(" \n*********************************\n  ");
    updateNonLinearJacobianTerms(this->M_jacobian,disp,dataMaterial,mapsMarkerVolumes, displayer);
    //displayer->leaderPrint(" \n*********************************\n  ");

}