omega_h
omega_h copied to clipboard
i just compiled it , but how to use it.
Dear sir: could you please provide some tutorial to show how to use it. Regards
Hello @ibaned Working with @delalond, I am willing to use omega_h for sub-molecular simulations at Blue Brain Project. Despite a very concise and straightforward code and unit-tests, I have some difficulties to get into it. Do you recommend some simulation projects using omega_h demonstrating how to use your library?
Best
Hi @tristan0x and @delalond, I have invited you to view https://github.com/ibaned/louie which is not yet globally public, but should be one good example. This is a small simulation project that implements the Lagrangian Finite Volume methods for shock hydrodynamics developed by this paper:
Carré, G., Del Pino, S., Després, B., & Labourasse, E. (2009). A cell-centered Lagrangian hydrodynamics scheme on general unstructured meshes in arbitrary dimension. Journal of Computational Physics, 228(14), 5160-5183.
And subsequently by P.H. Maire and others.
There is also a short integration here which uses very high-level APIs to adapt a mesh in Albany:
https://github.com/gahansen/Albany/blob/master/src/adapt/AAdapt_MeshAdapt_Omega_h.cpp
We are planning to open-source our main application very soon, and I'll post here when that is available to look at.
Hello @ibaned, I am a coworker of @tristan0x. Could you please also invite me to view https://github.com/ibaned/louie ?
Many thanks
@smelchio yes, I just did!
Thanks a lot for sharing your example on Lagrangian hydrodynamics. It is indeed interesting even if we are not very familiar with Lagrangian formulation with Finite Volume discretization. To start simple, we are thinking of coding a very simple example, typically solving a Laplacian with Finite Element or Finite Volume discretization. Would you have some guidance or an example to help use? We are especially interested in using shape functions with Omega_h.
Omega_h doesn't have much internal code for Finite Element shape functions, but it has general enough storage to attach degrees of freedom to edges/faces etc. It is more of a mesh data structure than a Finite Element framework, so you need to combine it with some Finite Element code. One example is the Albany integration I pointed out above, but we also have an integration with FEniCS which can be seen here:
https://github.com/ibaned/omega_h/blob/master/example/fenics/cpp/main.cpp
This is solving Poisson's equation with Omega_h adapting the mesh.
Hello @ibaned, I start this post. Could you please also invite me to view https://github.com/ibaned/louie ?
Many thanks
Dear @ztdepztdep
I know it has been a while but some tutorials have been added to this repository since then, available in example
directory.
I have run the examples in the directory. I didn't find the examples about mesh adapt. Could you please upload some examples.