Hi Edscott,
you are right, you have to recompile the complete Dune stack with Clang
and the flags for the desired saninitzer. When you use a different build
directory, it is easy to switch back to your GCC-built stack.
Hope this helps,
Christoph
Am 22.08.2018 um 18:23 schrieb Edscott Wilson:
>
Hi Birger,
I think you are right. In Dumux 3.0 the volume is scv.volume() though.
fvGeometry.subContVol[scvIdx].volume would be correct in Dumux 2.12 and
earlier.
Timo
Viele Grüße,
Timo
> Am 22.08.2018 um 13:08 schrieb Birger Hagemann
> :
>
> Dear Martin,
>
> The unit for the well
Hi, Christoph,
clang sounds good. I'll give it a try, But will I have to recompile dune
with clang or will it be compatible with gcc generated .a libraries?
Best Regards,
Edscott
El mar., 21 de ago. de 2018 a la(s) 17:10, Christoph Grüninger (
f...@grueninger.de) escribió:
> Hi Edscott,
>
Hi Timo,
Problem occurs exactly the same way with or without -O0. In order to
test whether problem only occurs with gdb, I add line 185 and run.
brookscorey.hh:
184swe = min(max(swe, 0.0), 1.0); // the equation below is only
defined for 0.0 <= sw <= 1.0
185std::cerr<<"swe
Dear Birger,
thanks for this comment, you are right.
You can directly use *scv.volume()* to get the volume.
Regards,
Martin
On 08/22/2018 10:38 AM, Birger Hagemann wrote:
Dear Martin,
The unit for the well index is m³. Consequently the unit for your
expression for values[contiWEEqIdx] is
Dear Martin,
The unit for the well index is m³. Consequently the unit for your expression
for values[contiWEEqIdx] is kg/s. However, the source term has the unit
kg/m³/s. I think you have to divide by the cell volume:
values[contiWEqIdx] = densityW/viscosityW * WI *(pbh -
Dear Ranjeet,
you can implement the Peaceman well model as a solution-dependent source,
using the following function in your problem file
*NumEqVector source(const Element ,**
** const FVElementGeometry&
fvGeometry,**
**