Re: [deal.II] Help with Vector template instantiations

Hello Wolfgang, 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 ==

Re: [deal.II] Help with Vector template instantiations

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

[deal.II] Re: Help with the weak form

Thank you Bruno. Dragan On Friday, September 9, 2016 at 1:07:19 PM UTC+1, Bruno Turcksin wrote: > > Dragon, > > On Friday, September 9, 2016 at 7:02:15 AM UTC-4, dragn...@gmail.com > wrote: >> >> Thank you Denis. >> >> I tried that. The function log(c/(1-c)) looks like: >> >> >>

Re: [deal.II] Help with Vector template instantiations

Hello Wolfgang, Yes, I do explicitly instantiate these templates in my code. The datatype is different from what is already defined in deal.II. Something like: #include #include template class Vector; But these two guys: dealii::Vector::operator=(std::complex) and

[deal.II] Re: Help with the weak form

Thank you Denis. I tried that. The function log(c/(1-c)) looks like: and it is difficult to linearise in the whole range. And the assembled matrices are solved out of

Re: [deal.II] Help with Vector template instantiations

Hello Wolfgang, Ok, good. The commit is fine. I'll see if I can pinpoint the exact cause of segmentation faults. One more question, in fact a small inconsistency that is of my interest: the function VectorTools::interpolate_boundary_values is not templated for the generic Number type but for

[deal.II] Re: Calculate surface integral over a boundary

Hello Jean-Paul, Thank you for the help. Well, the actual situation is more complex, requires a generic solution, but it simplifies into a calculation of total quantity through a boundary. For instance, a simple case: say we have heat or mass transfer problem (or any scalar quantity) and

[deal.II] Calculate surface integral over a boundary

Hello guys, I searched the forum about how to calculate the surface integral over a boundary (during the post-processing) and found one post, where he was directed to the step-12. I am not sure if I can understand from there. Any directions how to do it or some other place to look? Thanks.

Re: [deal.II] Help with Vector template instantiations

True, I've checked in the master and it is templated for the Function too. I use 8.4.1. It's fine then. Thanks. On Friday, September 16, 2016 at 5:39:33 PM UTC+1, Wolfgang Bangerth wrote: > > On 09/15/2016 05:12 AM, dragn...@gmail.com wrote: > > > > One more question, in fact a