[ff3d-dev] freefem3d unexpected error

[olivier Gipouloux]

hello,

i want to solve a lapalcian type problem in a square 3D box, but with some
special addon: i want to impose a flux jump on  surface cutting the box.

If the jump is constant, i know how to add a bidimentionnal linear form
but if the jump depend on a bilinearform on the surface, a can not add
a bidimensionnal bilinearform to thge tridimentionnal bilinear form (see 
the attachement
bellow)

Is it because it is not possible to solve this proble with freefem3d or 
did i make a mistake?

thank you

address@hidden:~/Documents/ff3d> ./ff3d-2005-12-03 essai.ff
=========================
== Text mode execution ==
=========================
FreeFEM3D computing thread: started
Parsing data

/home/delpino/src/ff3d/language/VariationalFormulaExpression.cpp:56: 
unexpected operator type




--
Olivier Gipouloux
Institut Camille Jordan UMR 5208
Batiment ISTIL
15 Boulevard Latarjet
69622 Villeurbanne Cedex
T� :  (33) 04 72 44 82 47 / Fax : (33) 04 72 43 11 45
https://mcs.univ-lyon1.fr/membres/gipouloux.html

union {
    box {
       <0,0,0>, <1,1,0.5> 
        } 
    box {
       <0,0,0.5>, <1,1,1> 
        } 
    pigment { color rgb <1,0,0> }
       }
intersection {
    box {
       <0,0,0>, <1,1,0.5> 
        } 
    box {
       <0,0,0.5>, <1,1,1> 
        } 
    pigment { color rgb <0,1,0> }
             }

// Definition de la boite //
vector n = (21,21,21);
vector a = (0,0,0);
vector b = (1,1,1);
mesh M = structured(n,a,b);
scene S = pov("essai.pov");
domain Omega = domain(S,inside(<1,0,0>));
domain Sigma = domain(S,inside(<0,1,0>));
mesh Q = surface(Sigma,M);
 solve (u) in Omega by M  krylov(precond=ichol)
          {
           test(v)
              int(dx(u)*dx(v) +dy(u)*dy(v) + dz(u)*dz(v) ) 
// the problematic line is bellow
            + int[Q](dx(u)*dx(v))
// but the line bellow do not give an error.....
           =  int[Q](10.*v);
//          =  0;
              u = 0 on M zmin;
              u = 5. on M zmax;
          }            

save(medit,"essai",M);
save(medit,"essai",u,M);

Olivier.Gipouloux.vcf
Description: Vcard