2021_05_26 v0.1.12 How to obtain a Michell structure? use; max_penalty = 3 Niter = 150 full_penalty_iter = Niter*1.1 ngauss = 1 S[y, x] = abs(s1) + abs(s2)

- V0.1.21 is the last containing code (although commented out) using the classical penalty formulation