Mika Malinen

I believe creating directly a square (N x N)-matrix, with N the count of edge element DOFs, would be the most natural approach. If one proceeds in a point-wise and...

To obtain a local interpolant having a guaranteed error estimate, I'd use a different approach and start by writing an additional integration loop in order to integrate the local mass...

It seems that the DOFs (the content of the array Asloc after the first BLOCK) are computed correctly. I verified this with some simple source fields which should be contained...

It seems that your code at the moment (https://github.com/ElmerCSC/elmerfem/commit/4c822dc8aac012d1e71a7e05469b37560a86cf84) supposes the variables `AsVar` and `BsVar` to have their value arrays suitable for representing the variable at integration points. However, they...

Hi, This would certainly be a nice feature! I think the linking of a material model should be possible without compiling with the utility elmerf90. The following check may not...

When seeking for procedures, Elmer has used Fortran mangling so that a symbol name has been expected to be in lower case letters and to be appended with one underscore....

The computation of stresses/strains/stress resultants hasn't been implemented in this case. This is of course a serious handicap and we should just find time to implement it ...

> why you guys cannot take ShellMultiSolver and based on it extend your native ShellSolver ? It has been derived from a very old solver (https://github.com/ElmerCSC/elmerfem/blob/devel/fem/src/modules/FacetShellSolve.F90) and doesn't therefore use...

Around lines 2175 it is now defined `cmplx_val = IP % s(t)*detJ*SUM(Wbasis(j,:)*gradv)` which is not feasible if dim == 2. If dim == 2, the value of cmplx_val should be...