ff3d-users
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[ff3d-users] "No cell has reference that matches" can be resolve, but...


From: Robert Li
Subject: [ff3d-users] "No cell has reference that matches" can be resolve, but...
Date: Mon, 26 Jan 2004 16:13:24 -0800 (PST)

Dear Stephane Del Pino,
  
  I can resolve the problem "No cell has reference
that matches (2)!", when I set the distance(from the
origin) is 0.2, not 0.5,0.4,0.3 in pov file. That is
to say I can not get the actual face of the box(z=-0.5
bottom, z=2 up). I don't know why? 
  (I used FDM method here, because I will use more
than one box, just like  soil layers, I want to
extract different layer boundary, then set boundary
condition). which method of FEM and FDM is possible?


Robert


>>>>>>POV file
box { <-0.5,  -0.5,  -0.5>, < 4, 4, 2  > pigment {
color rgb <1,0,0> } }
plane  { <0,0,-1>, 0.2 pigment { color rgb <1,0,1> }}
plane  { <0,0,1>, 1.7 pigment { color rgb <1,1,0> }}

>>>>>>>>input file 
vector n = (20, 20, 30);
vertex a = (-1, -1, -2);
vertex b = (5, 5, 8);

mesh M = structured(n,a,b);
 
scene S = pov("Tplane.pov");    // the pov-ray file for
the geometry
domain O = domain(S, inside(<1,0,0>) and 
inside(<1,0,1>) and  inside(<1,1,0>));

mesh Omesh = surface (O, M);
mesh s1 = extract(Omesh, 1);
mesh s2 = extract(Omesh, 2);
mesh s3 = extract(Omesh, 3);

solve(u) in O by M  //FDM  
{
   pde(u)
     - div(grad(u)) = 0;
 u=-10 on s3;                   //<1,1,0>;
 u=3 on s2;                             //<1,0,1>;
u=0 on M;
};
save(medit,"TT",u,M,dos);
save(medit,"TT",M,dos);
exec("medit TT");

__________________________________
Do you Yahoo!?
Yahoo! SiteBuilder - Free web site building tool. Try it!
http://webhosting.yahoo.com/ps/sb/




reply via email to

[Prev in Thread] Current Thread [Next in Thread]