/usr/share/doc/freefem++/examples/examples++/aalaplace-nc.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 | mesh Th=square(10,10);
fespace Vh(Th,P1nc);
Vh u,v;
u=0;
plot(u);
func f= 0;
func g= x*x+y*y*2;
u=0;
int i=0;
// dcl the problem
problem a(u,v,solver=CG,init=i) =
int2d(Th)( dx(u)*dx(v) + dy(u)*dy(v))
+ int2d(Th) ( v*f )
+ on(1,2,3,4,u=g+i);
cout << "-------------------" << i << endl;
i++;
a; // solve the problem
cout << "-------------------" << endl;
plot(u);
|