Elmer Discussion Forum

Hi, If the solver is never executed there will unfortunately be no stiffness matrix. You could however, try this keyword in Solver section: Linear System Solver Disabled = Logical True There might be other ways, like setting linear system iterations to zero etc. but this one should be clean. Codewis...

Hi I added a few days ago some more output (with high "max output level") so you might get a better idea now what is going on. Generally memory issues may be difficult and if it happens in netcdf library it is hard to tell. In Elmer it is a good idea to always use and check the optional ST...

Hi

Elmer has 2nd order edge elements (Nedelec's 1st kind) invoked by "quadratic approximation = true".

-Peter

Hi, can you solve anything with success? CUDA does not like DLL so that could be the culprint. However, in that case you should not be able to run most cases. There is something fishy since the code should be looking for "ProcName_bulk" and not just "_bulk". What is you setup and...

Maybe look at test case "ElasticBeamHinged3D". -Peter

Hi I think forces are always in cartesian directions. What the code does is it assembles the equations in cartesian coordinates. Then it rotates the entries when glueing them into the global matrix. The Dirichlet BCs are then set in this global matrix and hence there the "normal-tangential disp...

Hi

What line causes the segfault? Is it netcdf reader or Elmer side? If you increase verbosity to 30 what is the last line Elmer writes on std out?

-Peter

...your error probably has nothing to do with this. There are ~1000 tests that don't have _bulk defined but pass. As you see in the code there is no "procname" so somehow you have come here without giving the Procedure. Sharing the full case would indeed be helpful.

-Peter

Hi There have been different flavors of solver modules over the years. The standard solver does quite a bit of things that could be left for library. You could just provide rules for bulk and boundary assembly etc. The code in MainUtils.F90 around line 1250 checks this possibility: AssProc = GetProc...

Hi, One can also use iterative solvers for computing Gebhart factors. Standard "Linear System ..." keywords apply. Just use "radiation: " as a name space. E.g. radiation: Linear System Solver = iterative The computation is heavy since linear system is solved for each row. Direct ...