[deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

Artur, On Wednesday, May 9, 2018 at 8:51:37 PM UTC-4, Artur Safin wrote: > 1) Can FE_Nedelec > work > for unstructured adaptively refined meshes (for linear polynomials)? > 2) There is a new Nedelec FE class (pull 2240 >

[deal.II] Re: Initialising vector solution for mixed formulation problem

Hi Daniel, Thanks for your reply, the error i was getting was: void dealii::VectorTools::interpolate(const Mapping &, const DoFHandlerType &, const Function &, VectorType &, const dealii::ComponentMask &) [dim = 2, spacedim = 2, VectorType = dealii::BlockVector, DoFHandlerType = DoFHandl

Re: [deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

On Thu, May 10, 2018, at 08:10 CDT, Bruno Turcksin wrote: > The current version of FE_Nedelec should only be used for cartesian mesh. I have had no problems on unstructured, locally refined 2D meshes with our FE_Nedelec element so far. (With verified correct results when comparing against ana

Re: [deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

2018-05-10 16:35 GMT-04:00 Matthias Maier : > > We keep repeating that the current element only works on Cartesian > meshes, but I am not sure whether this is actually true. > Many people have reported problems see this thread for example https://groups.google.com/forum/#!searchin/dealii/nedelec%7

Re: [deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

On Thu, May 10, 2018, at 15:41 CDT, Bruno Turcksin wrote: > Many people have reported problems see this thread for example > https://groups.google.com/forum/#!searchin/dealii/nedelec%7Csort:date/dealii/ao-bi_Nh7gc/jHq1qMJlAgAJ > If I remember correctly, J-P couldn't make Nedelec work for his pr

Re: [deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

2018-05-10 16:52 GMT-04:00 Matthias Maier : > with the addition "in 2D" and for "lowest and second order" is true :-) > Yeah you are right, it's always in 3D that people have problems. Bruno -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https:/

[deal.II] Re: Initialising vector solution for mixed formulation problem

Jane, > the error i was getting was: > > void dealii::VectorTools::interpolate(const Mapping &, > const DoFHandlerType &, const Function VectorType::value_type> &, VectorType &, const dealii::ComponentMask &) > [dim = 2, spacedim = 2, VectorType = dealii::BlockVector, > DoFHandlerType = Do

[deal.II] ASPECT 2.0.0 released

We are pleased to announce the release of ASPECT 2.0.0. ASPECT is the Advanced Solver for Problems in Earth's ConvecTion. It uses modern numerical methods such as adaptive mesh refinement, multigrid, and a modular software design to provide a fast, flexible, and extensible mantle convection solver.

[deal.II] Re: Epetra_FEVector not being created as a LinearMap

I solved my problem. In the code I included below, I was using the Vector object from the namespace dealii::LinearAlgebraTrilinos::MPI. When I switch to using Vector objects from the namespace dealii::LinearAlgebra::EpetraWrappers (and add the trilinos_epetra_vector.h include file), the Epetr

Re: [deal.II] Nedelec Elements for unstructured adaptively refined meshes

> On 10 May 2018, at 22:41, Bruno Turcksin wrote: > > If I remember correctly, J-P couldn't make Nedelec work for his problem a few > months ago. Yes, that’s correct. I was working on a magneto-static problem in 3d, and wasn’t able to achieve sensible results. That said, I couldn’t say for ce

Re: [deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

The problem in the thread you cited is for 3D. I am happy to update our disclaimer for the class. I am just saying that our statement: "For the moment, it is safe to use the element on globally refined meshes with consistent orientation of faces." with the addition "in 2D" and for "lowe

Re: [deal.II] Re: Epetra_FEVector not being created as a LinearMap

On 05/11/2018 12:30 PM, josh.hano...@gmail.com wrote: I solved my problem. In the code I included below, I was using the Vector object from the namespace dealii::LinearAlgebraTrilinos::MPI. When I switch to using Vector objects from the namespace dealii::LinearAlgebra::EpetraWrappers (and add t