This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

How does dolfin handle box integration schemes?

–1 vote

Hi,

This is problably a naive question but I was wondering how dolfin handle box integration cell schemes. I know this is useful for some convection-diffusion problems where current (flux) conservation is strictly required. If someone has an example that would be useful.

Thanks,
Chaffra

asked Jan 13, 2014 by chaffra FEniCS User (1,830 points)
edited Jan 13, 2014 by chaffra

I'm not sure what do you mean. If you're asking for quadrature, it is handled by FFC in ffc/ffc/quadrature_schemes.py. Gauss quadrature is normally used.

commented Jan 14, 2014 by Jan Blechta FEniCS Expert (51,420 points)

Please have a look at the first paragraph at this link http://www.iue.tuwien.ac.at/phd/klima/node23.html . Wondering if this type of discretization is implemented or can be implemented in fenics.

commented Jan 15, 2014 by chaffra FEniCS User (1,830 points)

I don't think that FEniCS is designed so that this can be possible. But maybe I am wrong.

commented Jan 16, 2014 by Jan Blechta FEniCS Expert (51,420 points)
...