On 4/16/10, Juzar Thingna <
address@hidden> wrote:
> Hello Everyone,
> Here is my input file:
>
> vector n=(70,90,70);
> vector a=(0,7,0);
> vector b=(20,37,20);
> scene S=pov("ferromagnet.pov");
> mesh M=structured(n,a,b);
> domain O =domain(S,inside(<1,0,0>));
> function L=36;
> function wb=0.37824;
> solve(w) in O by M
> krylov(type=bicg)
> {
> pde(w)
> div(L*grad(w))-w=0;
> dnu(w)= wb on M ymin;
> w=0 on M ymax;
> }
> mesh Z = tetrahedrize (O,M);
> save(medit,"w_fm",w,Z);
> save(medit,"w_fm",Z);
> save(raw,"w_fm.dat",w,M);
>
> I want to set the derivative of w along the y-direction to be equal to wb at
> ymin. I thought that setting dnu(w)=wb would do the trick, but apparently I
> don't get the value of the slope at ymin = wb.
> Please could someone guide me as in where I'm going wrong?
>
> The foll is my geometry file:
>
> box {<0,7,0>,<20,37,20> pigment{color rgb <1,0,0>}}
>
> Regards,
> Juzar Thingna
> Center for Computation Science and Engineering
> National University of Singapore
>