Modification dans Laplace::solvgradconj

Une seule chose a changer: maintenir hconj[k] nul si k est un indice de vertex Dirichlet. Il faut donc remplacer le test

"if(v[i].where]"  par "if(fabs(u0[i])> precise)"

ou u0 contient la condition de Dirichlet. On rappelle la convention
u0[i] =0 entraine qi n'est pas un point de Dirichlet
On rappelle aussi que les constructeurs sont fait pour que u0 soit dans sol en entree dans solvgradconj. Il faut donc le stocker dans un tableau de travail (une fonction P1) avant de faire des iterations qui change sol.