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

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

On 09/23/2016 11:30 AM, dragniko...@gmail.com wrote: Oh yes, that' exactly what I wanted (in the post-processing case). Thanks. And if I want to do it during integration? What do you mean by "during integration"? What are you integrating at that time? Best W. -- ---

[deal.II] Re: Trilinos SparseMatrix mmult bug

Hello Bruno, Sorry for the vagueness, I was not confident if this was an actual bug in deal.ii (I've never encountered one before!). Here I've attached a bare-bones version of my code which reproduces the bug I'm encountering. I basically look at 4 cases: 1) deal.ii's TrilinosWrappers::SparseM

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

To answer my own question: finally I resolved a problem. I've been always getting zeros for the surface integral. It was a typo in the _expression_ I have been using (c is the dof, D is diffusivity): but I used q index for c by mistake. Thanks fo

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

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

On 09/23/2016 11:17 AM, dragniko...@gmail.com wrote: 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)

[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] Re: Calculate surface integral over a boundary

Hi Dragan, What is it that you're trying to evaluate? Although you don't specifically mention it, as you're wanting to do this in post-processing I presume that you're wanting to evaluate something to do with the solution? If so, then you can use the get_function_values

[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] Re: How can I have access to the all element of previous solution?

> > Thank you very much for your kind reply!!! Best, Kyusik. -- 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.I

Re: [deal.II] Re: How can I have access to the all element of previous solution?

Kyusik, 2016-09-23 9:11 GMT-04:00 : > In the lecture 31.65, the professor change the step-6.cc a little to apply > Picard iteration. > > Doing this, he uses previous_solution_gradients[q_index]. > > So, like this, I just used previous_solution[q_index] > > But, he also say that previous_solution_

[deal.II] Re: How can I have access to the all element of previous solution?

First of all, Thank you very much for your reply And, I'm sorry, I should have been more specific about my problem In the lecture 31.65, the professor change the step-6.cc a little to apply Picard iteration. Doing this, he uses previous_solution_gradients[q_index]. So, like this, I just

[deal.II] Re: How can I have access to the all element of previous solution?

Kyusik, On Friday, September 23, 2016 at 2:09:26 AM UTC-4, hanks0...@gmail.com wrote: > > In the lecture, previous_solution_gradients[] is used. But, I also want to > know How I can use just previous_solution[]. > > At first, I used previous_solution[q_index] in assemble_system as in > previous

[deal.II] Re: Trilinos SparseMatrix mmult bug

Ben, On Thursday, September 22, 2016 at 5:41:26 PM UTC-4, Ben Shields wrote: > > Hey there, thanks for the response. I made the recommended change to > trilinos_sparse_matrix.cc, and the corresponding test code returns the > appropriate output. However, I'm still getting an issue in my particula

[deal.II] Re: periodic boundary conditions

Hej, I found the issue ... I forgot the line constraints.distribute ( sol ); after the solution routine ... so it was just mainly a visualization issue. Thanks for the quick response! El jueves, 22 de septiembre de 2016, 19:26:28 (UTC+2), Daniel Arndt escribió: > > Juan, > > [...] >> 1) M