[Gmsh] Bug with high-order adaptive visualization

Marc Duflot marc.duflot at cenaero.be
Mon Sep 6 11:43:26 CEST 2010


Hi,

I think there is a bug in the high-order adaptive visualization
algorithm. No adaptation is performed if all the values are integers. It
is illustrated with this simple example:

--8<---

x1 = 1.01;
x2 = 1.00;  // No adaptation if it is an integer value

View "Quad" {
SQ(
0.0,0.0,0.0,
1.0,0.0,0.0,
1.0,1.0,0.0,
0.0,1.0,0.0
){
0.0, 0.0, x1, 0.0
};
SQ(
2.0,0.0,0.0,
3.0,0.0,0.0,
3.0,1.0,0.0,
2.0,1.0,0.0
){
0.0, 0.0, x2, 0.0
};
INTERPOLATION_SCHEME
{
  {1/4,-1/4,-1/4, 1/4},
  {1/4, 1/4,-1/4,-1/4},
  {1/4, 1/4, 1/4, 1/4},
  {1/4,-1/4, 1/4,-1/4}
}
{ {0,0,0},
  {1,0,0},
  {0,1,0},
  {1,1,0}
};
};

View[0].AdaptVisualizationGrid = 1;
View[0].MaxRecursionLevel = 5;
View[0].TargetError = 0;

--8<---


By the way, there is a mistake in the example in section 8.1 of the
manual. Lines 1 and 4 in the first brace of the interpolation scheme
need to be swapped, like I did in the example above.

Thanks,
Marc

-- 
Dr Marc Duflot -- CENAERO
Bâtiment ÉOLE,  1er étage
Rue des Frères Wright, 29
B-6041 Gosselies, Belgium

Tel. : +32 (0)71 91.93.64
Fax. : +32 (0)71 91.93.31

http://research.cenaero.be/people/duflot/