()
| 50 | l = [] |
| 51 | |
| 52 | def CalcError(): |
| 53 | space_flux.Update() |
| 54 | gf_flux.Update() |
| 55 | |
| 56 | flux = a.Flux(u) |
| 57 | gf_flux.Set (flux) |
| 58 | |
| 59 | err = (flux-gf_flux)*(flux-gf_flux) |
| 60 | elerr = Integrate (err, mesh, VOL, element_wise=True) |
| 61 | |
| 62 | maxerr = max(elerr) |
| 63 | l.append ( (v.ndof,sum(elerr)) ) |
| 64 | print ("maxerr = ", maxerr) |
| 65 | |
| 66 | for el in mesh.Elements(): |
| 67 | mesh.SetRefinementFlag(el, elerr[el.nr] > 0.25*maxerr) |
| 68 | |
| 69 | |
| 70 |
no test coverage detected