[deal.II] Computing the curl of a solution vector field obtained from Nedelec elements at receiver locations

2017-08-29 Thread Anna Avdeeva
Hi All, I have successfully computed the secondary electric field for 3D case as described in the paper by Grayver, 2014. These field is declared as BlockVector solution, where the block correspond to real and imaginary parts. Now I would like to compute total electric E (secondary plus

Re: [deal.II] Looking for more info on FE_RaviartThomas

2017-08-29 Thread Wolfgang Bangerth
Praveen, I have used FE_RaviartThomasNodal to solve induction equation on Cartesian meshes which worked out very well. Now I would like to solve Maxwell equations on general quad meshes. For this I have to use FE_RaviartThomas, but I do not completely understand this space. I'm not sure

[deal.II] Looking for more info on FE_RaviartThomas

2017-08-29 Thread Praveen C
Dear all I have used FE_RaviartThomasNodal to solve induction equation on Cartesian meshes which worked out very well. Now I would like to solve Maxwell equations on general quad meshes. For this I have to use FE_RaviartThomas, but I do not completely understand this space. I know that the dofs

Re: [deal.II] Re: Anisotropic fe space

2017-08-29 Thread Praveen C
Well, when I asked the question, I did not know how to do it. I then saw from my own answer to Jean-Paul's question that I could use a lower degree FE_RaviartThomasNodal itself to get test functions for moments. Thanks praveen On Mon, Aug 28, 2017 at 9:51 PM, Wolfgang Bangerth

Re: [deal.II] Best way to transfer old solutions in time stepping over mesh refinement?

2017-08-29 Thread Wolfgang Bangerth
On 08/29/2017 02:44 AM, 'Maxi Miller' via deal.II User Group wrote: and then I use the same grid both for the new and the old solution for the fe_values.grad_values(). Is that approach correct? I assume it still introduces some errors, For the record, yes, that's what we typically do --

[deal.II] deal.II Newsletter #3

2017-08-29 Thread Rene Gassmoeller
Hello everyone! This is deal.II newsletter #3. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #4986: added missing header (proposed by asartori86; merged)

[deal.II] Re: Best way to transfer old solutions in time stepping over mesh refinement?

2017-08-29 Thread Jean-Paul Pelteret
Dear Maxi, I'm going to put my foot down and say that, particularly as of late, you're not giving your problems sufficient attention before posting them to the forum. We have a set of guidelines https://groups.google.com/forum/#!topic/dealii/GRZMUTLIm2I (which, to the best of my recollection,

[deal.II] Best way to transfer old solutions in time stepping over mesh refinement?

2017-08-29 Thread 'Maxi Miller' via deal.II User Group
I am trying to implement example 26 using Newton's method and Sacado, in order to get a better feeling for the method. As mentioned before, I had problems with retaining the old solution when remeshing. Thus I rewrote my remeshing sequence as template void HeatEquation::refine_mesh ()

[deal.II] Re: Assembly of matrices for the heat equation when using Sacado generates -nan-values

2017-08-29 Thread Jean-Paul Pelteret
Thanks for giving feedback as to what the source of the issue was. I'm glad that you resolved your problem! On Tuesday, August 29, 2017 at 10:09:57 AM UTC+2, Maxi Miller wrote: > > Found my problem, it was the remeshing. Whenever I was remeshing, I > basically killed the old solution, which

[deal.II] Re: Assembly of matrices for the heat equation when using Sacado generates -nan-values

2017-08-29 Thread 'Maxi Miller' via deal.II User Group
Found my problem, it was the remeshing. Whenever I was remeshing, I basically killed the old solution, which then resulted in nan-values for it. That in turn generated -nan-values for the system matrix. Am Montag, 28. August 2017 09:52:45 UTC+2 schrieb Uwe Köcher: > > looks like your system