[ff3d-users] Referencing issue

[Arun Viswanathan]

Hi all, I have just got stuck with this code now. I
mailed in a few hours back too. I tried running this
code (both ff3d code and povray code attached below)
with 2 options of domains:
domain O = domain(S, outside(<0,1,0>); AND
domain O = domain(S, outside(<0,1,0>) and
outside(<0.9,0.9,0.9>);

Both give me the same error -- "No cell has reference
that matches (2)!"....Here (2) represents the
reference for the sphere (POV-ray to FEM translation).
The code runs error-free if I remove the 2nd boundary
condition (Robin's condition) but causes the above
errors if I include it. Please give me some
suggestions and guidance as I am unable to see what is
wrong in this code.

Thank you. Best regards,
Arun.


***************ff3d code below**********************

/* The program for the container with 12 drums inside
it. There is a small sphere next to the base of 1 of
the drums to depict the source of leakage from the
drum closest to it */

function D = 0.062;
function C0 = x*y*z;
function al = 0.1;

double dt = 0.1;
double i;

vector n = (10,10,10);
vector a = (0,0,0);
vector b = (2,2,6);
mesh M = structured(n,a,b);
scene S = pov("p2.pov");
domain O = domain(S, outside(<0,1,0>));

for (i=0; i<5; i++)
{
        cout << "\nIteration no. " << i+1 << "\n";
        solve(C) in O by M
        {
                pde(C)
                        C - div(D*dt*grad(C)) = C0;
                dnu(C) = 0 on M;
                al*C + dnu(C) = 0 on <0.9,0.9,0.9>;
        };
        C0 = C;
}
save(medit, "p2-5", M);
save(medit, "p2-5", C, M);

********************povray file below**************

union {
/*      box {
                <0,0,0>, <2,2,6>
                pigment { color rgb <0,0,0> }
        }
*/
        sphere {
                <0.05,0.05,0.5>, 0.05
                pigment { color rgb <0.9,0.9,0.9> }
        }

        cylinder {
                <0.5,0,0.5>, <0.5,1.4,0.5>, 0.4
                pigment { color rgb <0,1,1> }
        }

        cylinder {
                <1.5,0,0.5>, <1.5,1.4,0.5>, 0.4
                pigment { color rgb <0.9,1,1> }
        }

        cylinder {
                <0.5,0,1.5>, <0.5,1.4,1.5>, 0.4
                pigment { color rgb <0.8,0.9,1> }
        }

        cylinder {
                <1.5,0,1.5>, <1.5,1.4,1.5>, 0.4
                pigment { color rgb <0.7,0.8,1> }
        }

        cylinder {
                <0.5,0,2.5>, <0.5,1.4,2.5>, 0.4
                pigment { color rgb <0.6,0.7,1> }
        }

        cylinder {
                <1.5,0,2.5>, <1.5,1.4,2.5>, 0.4
                pigment { color rgb <0.5,0.6,1> }
        }

        cylinder {
                <0.5,0,3.5>, <0.5,1.4,3.5>, 0.4
                 pigment { color rgb <0.4,0.5,1> }
        }

        cylinder {
                 <1.5,0,3.5>, <1.5,1.4,3.5>, 0.4
                 pigment { color rgb <0.3,0.4,1> }
        }

        cylinder {
                 <0.5,0,4.5>, <0.5,1.4,4.5>, 0.4
                 pigment { color rgb <0.2,0.3,1> }
        }

        cylinder {
                 <1.5,0,4.5>, <1.5,1.4,4.5>, 0.4
                 pigment { color rgb <0.1,0.2,1> }
        }

        cylinder {
                 <0.5,0,5.5>, <0.5,1.4,5.5>, 0.4
                 pigment { color rgb <0,0.1,1> }
        }

        cylinder {
                 <1.5,0,5.5>, <1.5,1.4,5.5>, 0.4
                 pigment { color rgb <0,0,1> }
        }
        pigment { color rgb <0,1,0> }
}


        
                
__________________________________
Do you Yahoo!?
Friends.  Fun.  Try the all-new Yahoo! Messenger.
http://messenger.yahoo.com/