• Christophe Geuzaine's avatar
    · 3f96520c
    Christophe Geuzaine authored
    WARNING: potentially dangerous changes to the way nonlinear systems are assembled and solved:
    
    1) the computation of the residual DofData->res = b-Ax has been moved from SolveJac to GenerateJac
    
    2) the calculation of the full Jacobian matrix has also been moved from SolveJac to GenerateJac
    
    3) DofData->Jac now stores the full Jacobian (as soon as GenerateJac is done), instead of only the 
       JacNL terms ; DofData->A is not perturbed anymore
    
    This only applies if Flag_Only is not set.
    
    Performance used to be a problem with this approach in the past (PETSc 3.3) ; it seems OK now.
    
    Please test and report any breakage.
    
    Christophe & Kevin.
    
    3f96520c
magstadyn_a.pro 14.3 KB