Search found 4812 matches

by raback
07 Mar 2024, 14:42
Forum: External tools
Topic: Problem reading in netcdf file at higher resolution
Replies: 4
Views: 62

Re: Problem reading in netcdf file at higher resolution

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
by raback
07 Mar 2024, 14:37
Forum: ElmerSolver
Topic: Elmer-elmag example throws error
Replies: 21
Views: 336

Re: Elmer-elmag example throws error

...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
by raback
07 Mar 2024, 14:33
Forum: ElmerSolver
Topic: Elmer-elmag example throws error
Replies: 21
Views: 336

Re: Elmer-elmag example throws error

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...
by raback
04 Mar 2024, 11:36
Forum: ElmerSolver
Topic: Solver stops after Gebhard factors
Replies: 13
Views: 253

Re: Solver stops after Gebhard factors

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 ...
by raback
03 Mar 2024, 13:43
Forum: ElmerSolver
Topic: Solver stops after Gebhard factors
Replies: 13
Views: 253

Re: Solver stops after Gebhard factors

Hi The "ViewFactors" binary is called by "ElmerSolver" automatically. However, if the problem is the view factors then just compute those. If the view factors exist and the nodes are the same they are not recomputed unless requested. If you decleare your boundaries open then the ...
by raback
01 Mar 2024, 21:44
Forum: ElmerSolver
Topic: Solver stops after Gebhard factors
Replies: 13
Views: 253

Re: Solver stops after Gebhard factors

This piece of debug code in ViewFactors.F90 let's this case pass. Just checks the normal direction when R>20. ! This fixes the case by "Lhoeltke" on forum. Not a generic one! ! In line ~530 of ViewFactors. #if 1 BLOCK REAL(KIND=dp) :: r0(3), n0(3), rad,dotpr j = 0 k = 0 DO i=1,n Element =>...
by raback
01 Mar 2024, 21:41
Forum: ElmerSolver
Topic: Solver stops after Gebhard factors
Replies: 13
Views: 253

Re: Solver stops after Gebhard factors

Hi This is a strange one. The computation of view factors is one of the older pieces of Elmer and still there seems to be some strange bug. This piece of code has been used alot! I used quite a bit of time today to find the bug but I could not really find the main culprint. What happens is that ~25 ...
by raback
22 Feb 2024, 14:28
Forum: ElmerSolver
Topic: SIGSEGV Fault in parallel, but not serial
Replies: 7
Views: 127

Re: SIGSEGV Fault in parallel, but not serial

Hi Ryan, You jumped quickly to the parallel simulations of Elmer! Maybe you can show some nice color pictures. Do you have a reference for this case? It would be nice addition in "elmer-elmag" repo! Your setting's might not be fully optimal. For example, I think you didn't provide any BC's...
by raback
21 Feb 2024, 18:07
Forum: ElmerSolver
Topic: Elmer-elmag example throws error
Replies: 21
Views: 336

Re: Elmer-elmag example throws error

Hi, Linear solver setting are needed when solving for linear system. Usually this relates only to PDE's but it also relates to some postprocessing solvers where, for example, B=Curl(A) is also solved using the finite element method. No nonlinear system setting are needed for this since the primary s...
by raback
21 Feb 2024, 15:30
Forum: ElmerSolver
Topic: Elmer-elmag example throws error
Replies: 21
Views: 336

Re: Elmer-elmag example throws error

Hi It was previously possible not to define any linear solver settings. Then the code used defaults which happened to be "direct" solver with "banded" method. We thought that it was a very poor default so the user is now obliged to give one (for the CalcFields) solver. Linear Sys...