[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[ff3d-dev] freefem3d unexpected error
From: |
olivier Gipouloux |
Subject: |
[ff3d-dev] freefem3d unexpected error |
Date: |
Fri, 10 Mar 2006 12:36:55 +0100 |
User-agent: |
Mozilla Thunderbird 1.0.2 (X11/20050317) |
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
- [ff3d-dev] freefem3d unexpected error,
olivier Gipouloux <=