`gmsh -2 dual.geo`;
`getdp dual -mesh dual.msh -solve dual`;
`getdp dual -mesh dual.msh -pos err`;

for($ii=0; $ii<10; $ii++) {
	`gmsh -2 dual.geo -bgm errh.pos`;
	`getdp dual -mesh dual.msh -solve dual`;
	`getdp dual -mesh dual.msh -pos err`;
	`getdp dual -mesh dual.msh -solve dual -adapt errp.ada`;
	`getdp dual -mesh dual.msh -pos err`;
}

`getdp dual -mesh dual.msh -solve mwe`;
`getdp dual -mesh dual.msh -pos po`;
