Hello Stephane
thanks a lot for the vtk reader. It works perfectly !
I tried recently to use Robin type boundary conditions
for simple academic problems like:
vector n = (100,10,10);
vector a = (0,0,0);
vector b = (1,1,1);
mesh M = structured(n,a,b);
solve(u) in M
memory(matrix=none),
method(type=penalty)
{
pde(u)
-div(grad(u)) = 0;
dnu(u) + u = 2. on M xmax;
u=0. on M xmin;
};
save(vtk,"temperature",u,M);
it computes the right solution. But as soon
as I try to multiply "u" by something, I get
a syntax error.
dnu(u) + 3*u = 2. on M xmax;
maybe I missed something with the way the
boundary conditions have to be expressed ?
thanks for your advice
Best regards
Benoit
_______________________________________________
ff3d-users mailing list
[email protected]
http://lists.nongnu.org/mailman/listinfo/ff3d-users
