Elmer Discussion Forum

Hello, my PhD thesis is now available on: https://tel.archives-ouvertes.fr/tel-01305764 The Pdf is available here : https://tel.archives-ouvertes.fr/tel-01305764v3/document We are now able to detect accurately steady bifurcation in the case of 2D and 3D incompressible flows using ELMER. The text is ...

Hello, i managed to compile Elmer with a self-complied version of MUMPS_v5 using openBlas for MutliThread ability. It is not straigth forward but it's working... still have some stuff to check. I compiled Mumps as a ".so" shared library (as in a debian package). I added some information in...

Hello,

Have you tested this procedure using the new CMake installation?

Thanks.

Hi peter, you are right, i went first into this "solve with linear constraint" with your advice. That's why the code looks the same. As i didn't fully understand the solver and the fact that i use MUMPS as linear solver, i tried to create a minimal version (at least this is what i thought....

Hello, depending on the indices, AddToMatrixElement can be either normal or really slow. I want to create an augmented system with a vector at the last column of the MATRIX_LIST CollectionMatrix and the transposed vector at the last line like : | Stiff(1,1) Stiff(1,2) | Vec(1) | | Stiff(2,1) Stiff(2...

Hi, actually the -fPIC for MUMPS_5.0 compilation looks like removing some errors during compilation. I create a toolchain in order to make MUMPs_5.0 being found with my self-compile MUMPS_5.0: SET(MUMPSROOT /lab/limatb/app/MUMPS/MUMPS_5.0.0_FILES_fPIC_OpenBLAS) # - - - - - MUMPS INCLUDE - - - - - - ...

Thanks a lot for the trick, i'll try this asap.

Hi peter, You might also add Gmsh to ElmerGrid info for this 827 element at elmergrid/src/femfilein.c https://github.com/ElmerCSC/elmerfem/blob/devel/elmergrid/src/femfilein.c#L3065 static void GmshToElmerIndx(int elemtype,int *topology) ... int order827[]={0,1,2,3,4,5,6,7,8,11,13,9,10,12,14,15,16,1...

Thanks for the info, i'll try this ASAP, and let you know.

One should take care of permutation when everythin is ok..."Newind(1:n) = (/ (i, i = 1, n) /) ! nothing change" ! Write out the mesh ..... ..... IF( WriteData ) THEN IF( BinaryOutput ) WRITE( VtuUnit ) k DO i = ElemFirst, ElemLast IF(.NOT. ActiveElem(i) ) CYCLE CurrentElement => Model % El...