On Wed, Dec 8, 2010 at 10:13 AM, robert wrote:
> after the solve, the initial data are all set to one. I only played with
> the boundary conditions.
So you are solving ex9 (convection-diffusion) with discontinuous data
near the boundaries, and you are surprised there are oscillations?
AFAIK even
On Wed, Dec 8, 2010 at 9:53 AM, robert wrote:
> Hi,
>
> I have tried to modify example code 9. I am using a starting
> concentration of one everywhere but want to set the condition for one
> (or more) boundary to zero.
>
> for (unsigned int qp=0; qp {
> //const Numb