Dear Yves, Please see the attached files for the case when N=2 and N=3.
Best, Houtan On Wed, Nov 25, 2009 at 4:54 PM, Yves Renard <[email protected]>wrote: > > Dear Houtan, > > Strange. This error means that the acces to a dense matrix is wrong. > Of course, it does not inform of the function where the error is located. > To > have this information, it is necessary to recompile whith the -g option : > > rm test_range_basis.o > make test_range_basis CXXFLAGS="-g" > > run it with gdb : > > gdb test_range_basis > > then : > > run test_range_basis.param > > and when the error occurs : > > where > > > > > Yves. > > > On mercredi 25 novembre 2009, houtan wrote: > > Dear Yves, > > > > With N=2, the output from executing ./test_range_basis > > test_range_basis.param is: > > MESH_TYPE=GT_PK(2,1) > > FEM_TYPE=FEM_PK(2,1) > > INTEGRATION=IM_TRIANGLE(6) > > Selecting Neumann and Dirichlet boundaries > > Number of dof : 40401 > > Number of dof mult : 160801 > > Assembly of the mass matrix > > depart range basis > > Segmentation fault > > > > ------------------------------ > > > > With N=3, the output from executing ./test_range_basis > > test_range_basis.param is: > > MESH_TYPE=GT_PK(3,1) > > FEM_TYPE=FEM_PK(3,3) > > INTEGRATION=IM_TETRAHEDRON(6) > > Selecting Neumann and Dirichlet boundaries > > Number of dof : 97336 > > Number of dof mult : 226981 > > Assembly of the mass matrix > > depart range basis > > terminate called after throwing an instance of 'gmm::gmm_error' > > what(): Error in ../src/gmm/gmm_matrix.h, line 361 : > > out of range > > Abort > > > > Best, > > Houtan > > > > On Wed, Nov 25, 2009 at 2:31 PM, Yves Renard > <[email protected]>wrote: > > > Dear Houtan, > > > > > > Thank you for your contribution. > > > Can you try > > > ./test_range_basis test_range_basis.param > > > in the "tests" directory and report the error message ? > > > (The execution time is rather large, you can change N in > > > test_range_basis.param to N = 2). > > > > > > > > > Yves. > > > > > > On mercredi 25 novembre 2009, houtan wrote: > > > > Here are the steps I followed to build GetFem++ on Mac OS X 10.6.2 > > > > using gcc 4.2.1. There is one issue that remains to be resolved. > > > > > > > > 1) Download and install gfortran (gfortran-42-5646.pkg) from > > > > http://r.research.att.com/tools/ > > > > > > > > 2) ./configure CFLAGS='-arch i686' CXXFLAGS='-arch i686' > > > > CPPFLAGS='-arch i686' FCFLAGS='-arch i686' LDFLAGS='-arch i686' > > > > BLAS_LIBS='-framework Accelerate -lblas' --prefix='/usr/local' > > > > > > > > 3) Add the following to the top of src/getfem_interpolated_fem.cc > > > > #include <sys/types.h> > > > > #ifndef uint > > > > #define uint unsigned int > > > > #endif > > > > > > > > Equivalently, you can replace all instances of uint with unsigned or > > > > unsigned int in src/getfem_interpolated_fem.cc > > > > > > > > 4) make > > > > > > > > 5) make check > > > > > > > > 6) sudo make install > > > > > > > > The problem is that make check fails test/test_range_basis. Do you > know > > > > what the problem may be or do you have any suggestions of how to go > > > > about solving the problem? Every other test passes successfully. > > > > > > > > output: > > > > ./test_range_basis > > > > > > > /var/folders/jE/jEH8M1+hHsKcpoT6laWMak+++TI/-Tmp-/test_range_basis.param_ > > >24 > > > > > > >7 2>&1 failed > > > > FAIL: ../tests/test_range_basis.pl > > > > > > > > Best, > > > > Houtan > > > > > > -- > > > > > > Yves Renard ([email protected]) tel : (33) > 04.72.43.87.08 > > > Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29 > > > 20, rue Albert Einstein > > > 69621 Villeurbanne Cedex, FRANCE > > > http://math.univ-lyon1.fr/~renard<http://math.univ-lyon1.fr/%7Erenard>< > http://math.univ-lyon1.fr/%7Erenard> > > > > > > --------- > > > > -- > > Yves Renard ([email protected]) tel : (33) 04.72.43.87.08 > Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29 > 20, rue Albert Einstein > 69621 Villeurbanne Cedex, FRANCE > http://math.univ-lyon1.fr/~renard <http://math.univ-lyon1.fr/%7Erenard> > > --------- >
(gdb) run test_range_basis.param Starting program: /Users/houtan/Downloads/getfem-4.0.0/tests/test_range_basis test_range_basis.param Reading symbols for shared libraries .+++++.... done MESH_TYPE=GT_PK(2,1) FEM_TYPE=FEM_PK(2,1) INTEGRATION=IM_TRIANGLE(6) Selecting Neumann and Dirichlet boundaries Number of dof : 40401 Number of dof mult : 160801 Assembly of the mass matrix depart range basis terminate called after throwing an instance of 'gmm::gmm_error' what(): Error in ./gmm/gmm_matrix.h, line 361 : out of range Program received signal SIGABRT, Aborted. 0x91882732 in __kill () (gdb) where #0 0x91882732 in __kill () #1 0x91882724 in kill$UNIX2003 () #2 0x9191598d in raise () #3 0x9192ba44 in abort () #4 0x93711fda in __gnu_cxx::__verbose_terminate_handler () #5 0x9371017a in __cxxabiv1::__terminate () #6 0x937101ba in std::terminate () #7 0x937102b8 in __cxa_throw () #8 0x001627e0 in gmm::short_error_throw (file=0x1de774 "./gmm/gmm_matrix.h", line=361, func=0x1de2e8 "", errormsg=0x2140b4 "out of range") at gmm_except.h:75 #9 0x00164b31 in gmm::dense_matrix<double>::operator() (this=0xbfffe170, l=8088064, c=8088064) at gmm_matrix.h:362 #10 0x00035b87 in gmm::range_basis_eff_Gram_Schmidt_dense<gmm::transposed_row_ref<gmm::row_matrix<gmm::rsvector<double> >*> > (b...@0xbfffe5f4, colum...@0xbfffe320, c_ort...@0xbfffe364, EPS=9.9999999999999998e-13) at gmm_range_basis.h:331 #11 0x0003e950 in gmm::range_basis<gmm::transposed_row_ref<gmm::row_matrix<gmm::rsvector<double> >*> > (b...@0xbfffe5f4, colum...@0xbfffe5b4, EPS=9.9999999999999998e-13, skip_init=false) at gmm_range_basis.h:421 #12 0x0003edf9 in gmm::range_basis<gmm::transposed_row_ref<gmm::row_matrix<gmm::rsvector<double> >*> > (b...@0xbfffe5f4, colum...@0xbfffe5b4, EPS=9.9999999999999998e-13) at gmm_range_basis.h:481 #13 0x00004679 in laplacian_problem::assembly (this=0xbfffe6d4) at test_range_basis.cc:198 #14 0x00006dcf in main (argc=2, argv=0xbfffee50) at test_range_basis.cc:254
(gdb) run test_range_basis.param Starting program: /Users/houtan/Downloads/getfem-4.0.0/tests/test_range_basis test_range_basis.param Reading symbols for shared libraries .+++++.... done MESH_TYPE=GT_PK(3,1) FEM_TYPE=FEM_PK(3,3) INTEGRATION=IM_TETRAHEDRON(6) Selecting Neumann and Dirichlet boundaries Number of dof : 97336 Number of dof mult : 226981 Assembly of the mass matrix depart range basis terminate called after throwing an instance of 'gmm::gmm_error' what(): Error in ./gmm/gmm_matrix.h, line 361 : out of range Program received signal SIGABRT, Aborted. 0x91882732 in __kill () (gdb) where #0 0x91882732 in __kill () #1 0x91882724 in kill$UNIX2003 () #2 0x9191598d in raise () #3 0x9192ba44 in abort () #4 0x93711fda in __gnu_cxx::__verbose_terminate_handler () #5 0x9371017a in __cxxabiv1::__terminate () #6 0x937101ba in std::terminate () #7 0x937102b8 in __cxa_throw () #8 0x001627e0 in gmm::short_error_throw (file=0x1de774 "./gmm/gmm_matrix.h", line=361, func=0x1de2e8 "", errormsg=0x2140b4 "out of range") at gmm_except.h:75 #9 0x00164b31 in gmm::dense_matrix<double>::operator() (this=0xbfffe170, l=35446272, c=35446272) at gmm_matrix.h:362 #10 0x00035b87 in gmm::range_basis_eff_Gram_Schmidt_dense<gmm::transposed_row_ref<gmm::row_matrix<gmm::rsvector<double> >*> > (b...@0xbfffe5f4, colum...@0xbfffe320, c_ort...@0xbfffe364, EPS=9.9999999999999998e-13) at gmm_range_basis.h:331 #11 0x0003e950 in gmm::range_basis<gmm::transposed_row_ref<gmm::row_matrix<gmm::rsvector<double> >*> > (b...@0xbfffe5f4, colum...@0xbfffe5b4, EPS=9.9999999999999998e-13, skip_init=false) at gmm_range_basis.h:421 #12 0x0003edf9 in gmm::range_basis<gmm::transposed_row_ref<gmm::row_matrix<gmm::rsvector<double> >*> > (b...@0xbfffe5f4, colum...@0xbfffe5b4, EPS=9.9999999999999998e-13) at gmm_range_basis.h:481 #13 0x00004679 in laplacian_problem::assembly (this=0xbfffe6d4) at test_range_basis.cc:198 #14 0x00006dcf in main (argc=2, argv=0xbfffee50) at test_range_basis.cc:254
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
