Hello.

The time discretization in ff3d is let to the user. You have to write your 
scheme using the language. For instance, writing an implicit Euler scheme to 
solve the heat equation on a mesh 'M', one can write :

-----------------------------------
double dt = 0.1;
double i=0;

function un = 0;

for (i=0; i<10; ++i) {
  sove(u) in M
  {
    test(v)
    int(u*v) + int(dt*grad(u)*grad(v)) = int(un*v);
    u = 1 on 0; // dirichlet on boundary 0
  }
  un = u;
}
-----------------------------------

Best regards,
Stéphane.

_______________________________________________
ff3d-users mailing list
[email protected]
http://lists.nongnu.org/mailman/listinfo/ff3d-users

Reply via email to