Hi Gerardo,
"does not work" is kind of an useless problem description. What do you
expect? What happens actually?
Having a look into my crystal ball unveils the following hint:
- Adding a std::cout to your function to proof that it is actually
evaluated.
- Your grid must match the conditions. It must be fine enough that some
nodes are inside the described area and that this area is inside the
grid. Please double check. Maybe a std::cout inside the first case would
be helpful, too.
- Have a look at the resulting linear system, especially the right hand
side. Does your source term have an effect?
Bye
Christoph
--
... as Sir Cyril Hinshelwood has observed ...fluid dynamicists
were divided into hydraulic engineers who observed things that
could not be explained and mathematicians who explained things
that could not be observed. -- James Lighthill
*********************************************
CMWR 2014: 10th - 13th June 2014 in Stuttgart
Please visit www.cmwr14.de
*********************************************
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
