[Getfem-users] Problems with icare

[John Burnell]

Hi

I recently downloaded getfem++ 2.0.1 and compiled it. I tried a couple 
of the test programs and they worked okay. I then tried to run the 
navier-stokes solver icare (in contrib/icare). The program failed with a 
"dimension mismatch" error.

I tried to find the problem with the debugger and found a problem in the 
proper_update routine in icare.h - the indices vector was not being 
cleared so each time proper_update was called, indices grew larger and 
larger.

Fixing that problem gave another "dimension mismatch" error, but I don't 
understand this one. It occurs when on the first call of

getfem::standard_solve(MSM, mixed_dyn, iter);

on line 480 of icare.cc. This then results in a call sequence

mdpb.compute_tangent_matrix  at line 446 of getfem_model_solvers.h
pb.compute_tangent_matrix at line 363 of getfem_model_solvers.h
pb.compute_tangent_matrix at line 363 of getfem_model_solvers.h
((mdbrick_abstract*)sub_bricks[i])->compute_tangent_matrix(MS, i1, j1); 
at line 473 of getfem_modeling.h
((mdbrick_abstract*)sub_bricks[i])->compute_tangent_matrix(MS, i1, j1); 
at line 473 of getfem_modeling.h
do_compute_tangent_matrix(MS, i0, j0); at line 477 of getfem_modeling.h
gmm::copy(get_B(), gmm::sub_matrix(MS.tangent_matrix(), SUBI, SUBJ)); at 
line 1569 of getfem_modeling.h

Which is where the dimension mismatch occurs.  The first matrix has 
dimension of 786x1 and the second 1572x1.

At the moment I don't understand exactly what this sequence of calls is 
doing, so I thought I would ask if anyone knows what is going on?

I have compiled getfem++ under Linux and  Windows with Mingw and get the 
same error.

Thanks in advance for any help
John Burnell


--
==========================================
John G. Burnell
Research Scientist
Applied Mathematics Team
Industrial Research
P O Box 31-310, Lower Hutt
New Zealand.                       
Telephone: 64-4-9313-352       
Fax:       64-4-9313-003        
EMail:  address@hidden


This message has been scanned for viruses and dangerous content by MailScanner, 
and is believed to be clean.