/usr/share/doc/freefem++/examples/examples++-load/bilapMorley.edp is in freefem++-doc 3.19.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 | load "Morley"
int n=100,nn=n+10;
real[int] xx(nn),yy(nn);
mesh Th=square(40,40); // mesh definition of $\Omega$
fespace Vh(Th,P2Morley); // finite element space
// correct june 21 2007 (Thank to Hui Zhang hzhang@mailst.xjtu.edu.cn who find the error)
macro bilaplacien(u,v) ( dxx(u)*dxx(v)+dyy(u)*dyy(v)+2.*dxy(u)*dxy(v)) // fin macro
real f=1;
Vh [u,ux,uy],[v,vx,vy];
solve bilap([u,ux,uy],[v,vx,vy]) =
int2d(Th)( bilaplacien(u,v) )
- int2d(Th)(f*v)
+ on(1,2,3,4,u=0,ux=0,uy=0)
;
plot(u,cmm="u", wait=1,fill=1);
plot(ux,wait=1,cmm="u_x");
plot(uy,wait=1,cmm="u_y");
for (int i=0;i<=n;i++)
{
xx[i]=real(i)/n;
yy[i]=u(0.5,real(i)/n); // value of uh at point (0.5, i/10.)
}
plot([xx(0:n),yy(0:n)],wait=1);
|