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

Oups, wrong alloc for SVD. This is far better!

parent f5aa9258
No related branches found
No related tags found
No related merge requests found
...@@ -21,7 +21,8 @@ ...@@ -21,7 +21,8 @@
#define FORMAT_SMS 14 #define FORMAT_SMS 14
#define FORMAT_OPT 15 #define FORMAT_OPT 15
#define FORMAT_PSTEX 16 #define FORMAT_PSTEX 16
#define FORMAT_TEX 17 #define FORMAT_JPEGTEX 17
#define FORMAT_TEX 18
#define CONV_VALUE 0.8 #define CONV_VALUE 0.8
......
// $Id: Utils.cpp,v 1.5 2001-11-08 10:16:41 geuzaine Exp $ // $Id: Utils.cpp,v 1.6 2001-11-19 13:43:17 geuzaine Exp $
#include "Gmsh.h" #include "Gmsh.h"
#include "Numeric.h" #include "Numeric.h"
...@@ -47,7 +47,7 @@ void MeanPlane(List_T *points, Surface *s){ ...@@ -47,7 +47,7 @@ void MeanPlane(List_T *points, Surface *s){
na = 3; na = 3;
U = dmatrix(1,ndata,1,na); U = dmatrix(1,ndata,1,na);
V = dmatrix(1,ndata,1,ndata); V = dmatrix(1,na,1,na);
W = dvector(1,na); W = dvector(1,na);
for (i=0; i<ndata; i++){ for (i=0; i<ndata; i++){
...@@ -79,7 +79,7 @@ void MeanPlane(List_T *points, Surface *s){ ...@@ -79,7 +79,7 @@ void MeanPlane(List_T *points, Surface *s){
norme(res); norme(res);
free_dmatrix(U,1,ndata,1,na); free_dmatrix(U,1,ndata,1,na);
free_dmatrix(V,1,ndata,1,ndata); free_dmatrix(V,1,na,1,na);
free_dvector(W,1,na); free_dvector(W,1,na);
// check coherence of results for non-plane surfaces // check coherence of results for non-plane surfaces
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment