[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/