Elmer Discussion Forum

Hey Peter, Thanks for the info. In my case, it's failing in 3D elements, and as far as I can tell, they don't appear to be degenerate. I don't know enough about the transformation between global and local coordinate systems, but is it possible that a given element could be sufficiently distorted fro...

Hi all, A quick search of the forums reveals that the warning "GlobalToLocal: did not converge" has come up in a few instances. I'm experiencing this error myself, and I'm wondering if anyone can help me understand the root cause of this? The backtrace for this error is roughly: GlobalToLo...

Hi, I believe there's a problem with the SaveBoundaryValues routine when restarting from a run in which SaveBoundaryValues was already used, and all the requested Parameters already exist as variables. In this situation, values are never updated and so any VTU, SaveLine or result file output shows t...

Hi, I found the solution to this problem, so I thought I'd post it here in case it helps anyone else. GeneralUtils.src provides the following: CALL FreeQuadrantTree( MeshA % RootQuadrant ) which recursively DEALLOCATEs the Quadrant Tree, freeing the memory back to the system. This solved my memory l...

Hi all, I've not put this under 'bug reports' because it could well be my own code which is at fault. I've been experiencing very significant memory leaks from a long run. The program leaks ~27GB of RAM over 1 week. Using valgrind, I've pinpointed the problem to BuildQuadrantTree, in Interpolation.f...

Hi all, This isn't really a bug, so I'm posting it here instead, in case anyone else comes searching for an answer to quite an obscure issue to do with User Functions. I have found that if one writes a function to prescribe Mesh Update in some way, e.g. Mesh Update 1 = Variable Coordinate 2 Real Pro...

Hi Peter,

Thanks for this, good advice to get me on the right track. I'll have a think about your comments.

Joe

Hi all, I have a problem that I think might be solvable with the level set method, but I thought I'd see if anyone has any advice first. Basically, I have two field variables and I'd like to check if there's anywhere in the mesh where both are greater than zero. At individual nodes, this is trivial;...

Hi Martina, This sounds like a problem I had with a solver I wrote, which involved the VariableAdd function. Basically, the solver called: VariableAdd( variables, mesh, SolverPointer, Name, 1, Values, Perm ) On the assumption that this call was successful, the values were simply updated with: Values...

Hi all, I'm seeking some advice on linking the residual of the limited heat equation to subglacial hydrology at the bed. My model setup has a Dirichlet temperature condition at the bed (Pressure Melting Point) and, if I understand correctly, the heat sinks and sources imposed in satisfying this Diri...