dotest

FreeFem++可以生成高质量的有限元网格。可以用于流体力学

rm -f *[0s].mesh
g=
quad=
if [ "$1" = "-g" ] ;then g=-g;shift; fi
if [ "$1" = "-q" ] ;then quad="-2q -coef 2 -thetaquad 10" ;shift; fi
if [ "$1" = "-g" ] ;then g=-g;shift; fi

if [ -z "$bamg" -o ! -x "$bamg" ] ;then bamg=../../bamg$g; fi
echo " to close the graphic enter f on graphic window until they exist "
echo "Mesh  square " 
echo ${bamg} -g  square_g.mesh -o square_0.mesh $quad
${bamg} -g  square_g.mesh -o square_0.mesh $quad
if [  $? -ne 0 -o ! -s square_0.mesh ] ; then echo "the file square_0.mesh  is empty or d'ont exist";exit 1 ; fi

echo ${bamg} -g  square_g.mesh -o square_0.mesh $quad
${bamg} -splitpbedge  -g  square_g.mesh -o square_s.mesh  $quad
if [  $? -ne 0 -o ! -s square_0.mesh ] ; then echo "the file square_s.mesh  is empty or d'ont exist";exit 1 ; fi

echo "Mesh  circle  " 
${bamg} -g  circle_g.mesh -o circle_0.mesh $quad
if [ $? -ne 0 -o  ! -s circle_0.mesh ] ; then echo "the file circle_0.mesh  is empty or d'ont exist";exit 1 ; fi

echo "Mesh  a octogone  " 
${bamg} -g  octogone_g.mesh -o octogone_0.mesh $quad
if [ ! -s octogone_0.mesh ] ; then echo "the file octogone_0.mesh  is empty or d'ont exist";exit 1 ; fi

echo "Mesh  a  square with raff  " 
${bamg} -g  square_raf_g.mesh -o square_raf_0.mesh $quad
if [ $? -ne 0 -o ! -s  square_raf_0.mesh ] ; then echo "the file  square_raf_0.mesh  is empty or d'ont exist";exit 1 ; fi

echo "Mesh  a  square with raff  split pb edge" 
${bamg} -splitpbedge -g  square_raf_g.mesh -o square_raf_s.mesh $quad
if [ $? -ne 0 -o ! -s  square_raf_0.mesh ] ; then echo "the file  square_raf_0.mesh  is empty or d'ont exist";exit 1 ; fi

echo "Mesh  a  square with anisotropique mesh  " 
${bamg} -g  square_raf_g.mesh -o square_raf_ani_0.mesh -M square_raf_ani_g.mtr $quad
if [ $? -ne 0 -o ! -s  square_raf_ani_0.mesh ] ; then echo "the file  square_raf_ani_0.mesh  is empty or d'ont exist";exit 1 ; fi

echo "Mesh  a  square with anisotropique mesh with -splitpbedge " 
${bamg} -splitpbedge -g  square_raf_g.mesh -o square_raf_ani_s.mesh -M square_raf_ani_g.mtr $quad
if [ $? -ne 0 -o ! -s  square_raf_ani_s.mesh ] ; then echo "the file  square_raf_ani_s.mesh  is empty or d'ont exist";exit 1 ; fi

${bamg} -splitpbedge -g  l_g.mesh -o l_0.mesh  $quad
if [ $? -ne 0 -o ! -s  l_0.mesh ] ; then echo "the file   l_0.mesh  is empty or d'ont exist";exit 1 ; fi