Dragan,
It is ok, you can apply patch yourself.
I've committed under your name (which hopefully I got right):
https://github.com/dealii/dealii/pull/3114
I actually also had to specialize many template functions to work with the new
type (those related to matrices and vectors, BlockMatrixBase::add,
SparseMatrix::add and also local_apply_boundary_conditions etc.). There are
some functions that do not work well with the user-defined types like
Vector<>::allocate(), Vector<>::deallocate() and Vector<adouble>::operator
=(). They cause seg. faults. I am not sure why, perhaps because of mem.align
you use (Utilities::System::posix_memalign) and it is fine for primitive
types. Anyway that is something internal for my case.
Hm, good question what is happening there. Feel free to propose individual
patches for each of these cases if you can identify an underlying reason.
I am finishing what I started doing and will write you about an interesting
application of deal.II. It is mostly useful for multi scale modelling in
chemical engineering but is general in nature. It is coupling of deal.II with
an equation-based simulator: basically, using deal.II only to assemble
matrices/rhs (including non-linear FE cases), generating equations and then
solving one or more of these systems together with other differential and
algebraic equations in a large DAE system. All that is done from python
(although implemented in c++). An example could be Lithium-ion batteries:
electrodes are typically made out of a porous material composed of large
numbers of solid particles, and particles are in a electrolyte. Particles
(transport modelled using the FE method) are coupled at the larger electrode
length scale via electrolyte transport (using the finite volume method). Two
software are fully coupled and the boundary conditions can be set using the
equations from the other software etc.
That would make for a fantastic code gallery project!
Best
Wolfgang
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
