Re: [Libmesh-users] oscillations boundary conditions example code 9

2010-12-08 Thread John Peterson
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

Re: [Libmesh-users] oscillations boundary conditions example code 9

2010-12-08 Thread John Peterson
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