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

Well, as I said deal.II is integrated into another software (DAE Tools, mostly chemical engineering applications i.e. the fluid flow coupled with several scalar transport equations). Now, it works well. Basically, the system assembly is done only once (no adaptive meshes used). The system matri

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

Oh yes, that' exactly what I wanted (in the post-processing case). Thanks. And if I want to do it during integration? On Friday, September 23, 2016 at 6:21:23 PM UTC+1, Wolfgang Bangerth wrote: > > On 09/23/2016 11:17 AM, dragn...@gmail.com wrote: > > > > Well, the actual situation is more comp

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

And one more thing I forgot. Although it is typically a post-processing task I would like to get it during the integration (the problem is transient). On Friday, September 23, 2016 at 5:56:16 PM UTC+1, Jean-Paul Pelteret wrote: > > Hi Dragan, > > What is it that you're trying to evaluate? Althou

[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 spec

[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. Drag

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 small inconsis

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 do

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 t

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 == false

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

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. Dragan On Saturday, September 10, 2016 at 5:42:01 PM UTC+1, Wolfgang Bangerth wrote: > > On 09/10/2016 10:33 AM, dragn...@gmail.com wrote:

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 dealii::Vector::lp_norm(int

[deal.II] Help with Vector template instantiations

Hello, I need to use Vector FullMatrix, BlockVector and BlockSparseMatrix with different datatype than what is already provided and instantiated (float, double...). I have the following problem when instantiating with a different type: the linking fails due to the multiple definitions of the f

[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: >> >> >>

[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 deal

[deal.II] Help with the weak form

Hello guys, I would need a small help with deriving the weak form. The equation without the unnecessary terms is the following: log(c/(1-c)) - k*Laplacian(c) = 0 I am not sure how to treat the term log(c/(1-c)). Thank you. Regards, Dragan -- The deal.II project is located at http://www.dea