[Getdp] mesh adaptation

Christophe Geuzaine cag32 at case.edu
Sat Sep 24 03:42:18 CEST 2005


Michal Hajek wrote:
> Hello :)
> 
> I am trying to figure out how to make adaptation of mesh. 
> First I solve simple electrostatics problem and I get correct scalar
> potential on some mesh. Than I compute flux resulting from the potential
> gradient. So far so good. 
> 
> Than I wanted in PostOperation :
> 
> { Name Tar; NameOfPostProcessing MyPostProcessing;
>      Operation {
>             Print [ Flux, OnElementsOf Domain_Elektrika, Dimension 2,
>             Adapt P1, Target 300, Value 3, Format Gmsh, File
>             "refined.msh" ];
> }}
> 
> I thought that I make a very simple approximation of the error - just
> the flux itself. Unfortunately, if Target value is too small (eg.
> smaller than 40 in my case) I receive 
> this error:
> 
> GetDP   : Error     : MIN1D not converged: f(3.25927) = -71350.2
> 
> 
> or if Target value is too large (eg. >40 in my case) I receive another
> error:
> GetDP   : Error     : GSL: function not continuous (fsolver.c, line 40)
> 
> 
> I am using getdp version 1.0.1 compiled with gsl-1.6 on a debian amd64
> mashine. 


Could you try again with getdp 1.1.2?

Christophe


> 
> If someone is interested, I have put .geo, .msh, .pro and solver.par
> files here:
> http://material.karlov.mff.cuni.cz/people/hajek/getdp/
> 
> Thank you for any kind of help :)
> 
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp
> 


-- 
Christophe Geuzaine
Assistant Professor, Case University, Mathematics
http://www.case.edu/artsci/math/geuzaine