finish L1-normal and area-keep function, work on cylinder/pie/mobius/fandisk/ellipse. Since I use ipopt for L1 term, it seems quite slow. And for sculpture, I should use very small weight for linear equation fitting, actually zero