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

2016-09-16 Thread dragnikolic
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

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

2016-09-16 Thread Wolfgang Bangerth
On 09/15/2016 05:12 AM, dragniko...@gmail.com wrote: 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 double only (argument

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

2016-09-15 Thread dragnikolic
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

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

2016-09-14 Thread Wolfgang Bangerth
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,

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

2016-09-14 Thread Bruno Turcksin
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.

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

2016-09-14 Thread dragnikolic
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

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

2016-09-14 Thread Wolfgang Bangerth
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)

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

2016-09-12 Thread dragnikolic
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

2016-09-10 Thread Wolfgang Bangerth
On 09/10/2016 11:22 AM, dragniko...@gmail.com wrote: Yes, no problem. I see those files changed a lot since the 8.4.1 release. Anyway, I see the two spots are the same as in 8.4.1. I'll send you the patch. Fantastic, thanks! W. --

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

2016-09-10 Thread Wolfgang Bangerth
On 09/10/2016 10:33 AM, dragniko...@gmail.com wrote: 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:

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

2016-09-10 Thread dragnikolic
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