running

acceptance/facebody/sin50xy/accept-sin50xy.sh

@@ -11,33 +11,35 @@ else
 fi
 
 tecplot=-t
-field="-u 5"
+field="-u sin50xy"
 egads="-g square.egads"
 
 function adapt_cycle {
     inproj=$1
     outproj=$2
+    complexity=$3
 
     ${src}/ref_acceptance ${field} ${inproj}.meshb \
 	  ${inproj}.solb
 
     ${src}/ref multiscale ${inproj}.meshb ${inproj}.solb \
-	  1000 ${inproj}-metric.solb
+	  ${complexity} ${inproj}-metric.solb
 
     ${src}/ref adapt ${inproj}.meshb ${egads} -m ${inproj}-metric.solb \
 	  -x ${outproj}.meshb -f ${outproj}.tec
 
     ${src}/ref_acceptance ${field} ${outproj}.meshb \
 	  ${outproj}.solb
+
+    ${src}/ref visualize ${outproj}.meshb ${outproj}.solb ${outproj}.plt
 }
 
 cp square.meshb cycle00.meshb
 
-adapt_cycle cycle00 cycle01
-adapt_cycle cycle01 cycle02
-adapt_cycle cycle02 cycle03
-adapt_cycle cycle03 cycle04
-adapt_cycle cycle04 cycle05
-adapt_cycle cycle05 cycle06
+adapt_cycle cycle00 cycle01 1000
+adapt_cycle cycle01 cycle02 1000
+adapt_cycle cycle02 cycle03 2000
+adapt_cycle cycle03 cycle04 4000
+adapt_cycle cycle04 cycle05 8000