Comments (1)
Hi, what you observe is a totally correct behavior. The key is to understand how L2 elements work. Unlike continuous elements, L2 elements are not precisely defined at boundaries of elements and the notion of trace is used instead, which is not defined uniquely and different common choices exist. Practically, this means "open" bases are used for them typically, like the default Gauss-Legendre, which do not have DOFs at the boundaries. This implies that all DOFs essentially contribute to interpolation of the values (even though we may expect the closest contribute the most). This is the reason why the integration is performed on the whole element.
Projection to H1 might be an option, but keep in mind it is a global procedure, so it is a kind of integral fitting. This might cause artificial dislocation of some values, which might or might not be welcomed depending on the application (i.e., do you expect continuous or discontinuous profiles?). I would try to think differently maybe, I do not know your application exactly, but cannot you reformulate the problem in the weak sense? Like instead of calculation of strong divergence/curl, calculate them weakly? That is how L2 elements are typically used in EM codes 😉
from mfem.
Related Issues (20)
- Boundary Condition E=(1,0,0) for Nedelec elements for tetrahedral meshes HOT 2
- Fast Marching Method with adaption HOT 1
- Parallel solvers for the wave equation HOT 2
- contact case: tribol ,build problem HOT 1
- Boundary Dof (id) and physical location HOT 1
- Is is possible to implement AssemblePABoundary for VectorFEMassIntegrator
- Boundary conditions using the Gradientintegrator in the navier mini app
- Some SWIG Typemap is brokend due to real_t (instead of double) HOT 2
- Preventing silent failures using full assembly on very connected meshes HOT 6
- Removing sundials prefix from executable name is breaking some ninja builds HOT 3
- Which solver to use with MixedVectorWeakDivergenceIntegrator? HOT 2
- Bug in VTK reader or writer? HOT 5
- Possible issue in 3D in det.cpp? HOT 4
- add utf-8 option for MSVC compiler to fix build fail with MSVC and code page 936 HOT 2
- add utf-8 option for MSVC compiler to fix build fail with MSVC and code page 936 HOT 1
- ex5 and ex5p do not work with -pa when using square-disc.mesh HOT 1
- MakeSimplicial corrupt corrdinates with surface mesh HOT 1
- Segmentation Fault Error if Use `AssemblyLevel::FULL` in ex26
- Segmentation fault when run code with cuda in hypre HOT 1
- `FindPointsGSLIB` behavior when `point_pos` is empty on a rank.
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from mfem.