Skip to content
Snippets Groups Projects
Commit 37bd2a8e authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

*** empty log message ***

parent e4818098
No related branches found
No related tags found
No related merge requests found
......@@ -540,7 +540,8 @@ class CurvatureField:public Field
options["IField"] = new FieldOptionInt(iField);
options["Delta"] = new FieldOptionDouble(delta);
}
void grad_norm(Field &f,double x,double y,double z, double *g){
void grad_norm(Field &f,double x,double y,double z, double *g)
{
g[0]=f(x+delta/2,y,z)-f(x-delta/2,y,z);
g[1]=f(x,y+delta/2,z)-f(x,y-delta/2,z);
g[2]=f(x,y,z+delta/2)-f(x,y,z-delta/2);
......@@ -591,7 +592,8 @@ class MaxEigenHessianField : public Field
eigenvalues=gsl_vector_alloc(3);
gslmat=gsl_matrix_alloc(3,3);
}
~MaxEigenHessianField(){
~MaxEigenHessianField()
{
gsl_eigen_symm_free(gslwork);
gsl_vector_free(eigenvalues);
gsl_matrix_free(gslmat);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment