Startpage >> Main >> PLap

P Lap

Computation of P-Laplace operator; thanks to Frank Schadt ( frankschadt@web.de )

Iter1end

download example: plaplace2.edp or return to main page

verbosity=0.;
real Dx=.02, R=1.;
real p=1.5, gamma=0.1;
border C(t=0.,2.*pi){x=R*cos(t);y=R*sin(t);label=1;};
mesh Th=buildmesh(C(floor(2.*pi*R/Dx)));
fespace Vh(Th,P1);
Vh uh, uh0=0, V=1, vh;
macro Grad(u)[dx(u),dy(u)]//

problem probdup(uh,vh) = int2d(Th)(V*Grad(uh)'*Grad(vh))
         - int2d(Th) ( vh ) + on(1,uh=0);
for (int i=0;i<=1000;i++) {
   probdup;
   V=gamma*V+(1-gamma)*(Grad(uh)'*Grad(uh))^(p/2-1);
plot(uh,cmm="iter="+i+";max="+uh[].max+";p="+p+";gamma="+gamma,fill=true,value=true,dim=3,wait=0,boundary=true);
}

return to main page

Page last modified on January 07, 2015, at 11:24 AM