GenerateJac not working when GetDP compiled with SPARSKIT
When GetDP is compiled with SPARSKIT, the GenerateJac operation always triggers a GetDP error "Wrong arguments in LinAlg_AddMatrixMatrix()". I think this is due to a mistake in the function call at SolvingOperation.cpp:732 :
LinAlg_AddMatrixMatrix(&DofData_P->A, &DofData_P->Jac, &DofData_P->Jac) ;
Arguments 1 and 3 are not the same as they are expected to be, judging from the function definition in LinAlg_SPARSKIT.cpp :
void LinAlg_AddMatrixMatrix(gMatrix *M1, gMatrix *M2, gMatrix *M3)
{
if(M3 == M1)
add_matrix_matrix(&M1->M, &M2->M) ;
else
Message::Error("Wrong arguments in 'LinAlg_AddMatrixMatrix'");
}
I would say this is a simple inversion of the first and second argument.