use sor as PC for gamg

Author: gtheler

pipe.fee

@@ -1,10 +1,11 @@
 # convergence study of linearized stresses in an infinite pipe
 # with respect to the number of elements in the pipe thickness
-PROBLEM mechanical 3D
+PROBLEM mechanical 3D PC gamg KSP cg
+PETSC_OPTIONS -mg_levels_pc_type sor
 
 t0 = wall_time()     # initial wall time
 
-# read mesh aaccording to shape $1, order $2 and number of elements through thickness $2
+# read mesh according to shape $1, order $2 and number of elements through thickness $2
 READ_MESH pipe-$1-$2-$3.msh 
 
 INCLUDE problem.fee  # read problem parameters

run.sh

@@ -49,7 +49,7 @@ for shape in ust10 uct10 sst10 sct10 ssh20 sch20 ssh27 sch27; do
    if [ ! -e pipe-${shape}-${order}-${n}.msh ]; then
     gmsh -v 0 -3 pipe-${shape}-${order}-${n}.geo
    fi
-   feenox pipe.fee --mumps ${shape} ${order} ${n} | tee -a convergence-${shape}-${order}.dat
+   feenox pipe.fee ${shape} ${order} ${n} | tee -a convergence-${shape}-${order}.dat
 
    # create a geo to draw a picture of the mesh
    # showing how many elements are there with Gmsh