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, BlockMatrix
Dragan,
On Wednesday, September 14, 2016 at 7:34:46 AM UTC-4, dragniko...@gmail.com
wrote:
> 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. I
Hello Wolfgang,
It is ok, you can apply patch yourself.
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 t
Dragan,
Here is the patch. I added one more change in functions Vector::reinit:
if (omit_zeroing_entries == false)
*this = static_cast(0);
The line above fails if the Number type is not primitive. If it is a class it
can't cast. So it could be:
if (omit_zeroing_entries == false)
*
My program uses Parallel::distributed::triangulation.
If I want to output a data array associated with cell, which one of the
following is the correct way of doing it ? And why ?
(This cell data type vector is to be viewed in Paraview.)
|
intindex =0
for(;cell!=endc,++cell,++index)
if(cell->is