Search found 4 matches

by Victor
09 Oct 2020, 17:35
Forum: Bug reports
Topic: Referencing a pointer not allocated
Replies: 0
Views: 5305

Referencing a pointer not allocated

Running AdvDiffFluxes in VtuOutputSolver.F90:1653

InvFieldPerm => InvNodePerm

The pointer InvNodePerm is not allocated.
What is the fix, maybe NULLIFY InvFieldPerm?
Please advise
by Victor
09 Oct 2020, 17:32
Forum: Bug reports
Topic: zero**zero in ElementDescription.F90 [FIXED]
Replies: 2
Views: 3132

Re: zero**zero in ElementDescription.F90 [FIXED]

There is at least another place in ElementDescription.F90 computing 0**0 at line 1008.
But by now I know how to fix it.
However running AdvDiffFluxes I run into another bug, that I am describing in another issue
concerning referencing a pointer not allocated.
by Victor
28 Sep 2020, 08:50
Forum: Bug reports
Topic: zero**zero in ElementDescription.F90 [FIXED]
Replies: 2
Views: 3132

zero**zero in ElementDescription.F90 [FIXED]

Running several test cases, one of them is AdvDiffFluxes, at ElementDescription.F90 lines 789 and 826 s = s + Coeff(i) * u**p(i) both u and p(i) are zero, thus trying to compute zero**zero. Easy to check by putting if(u==0.and.p(i)==0) call Fatal('ElementDescription','Attempted computation of 0**0')...
by Victor
26 Sep 2020, 12:16
Forum: Bug reports
Topic: Loss of precision in SaveScalars.F90 [FIXED]
Replies: 1
Views: 2827

Loss of precision in SaveScalars.F90 [FIXED]

On elmerfem 8.4, running test case mgdyn2D_compute_complex_power at SaveScalars.F90:1565

TmpValues(1:n) = Values(1:n)

TmpValues is REAL while Values is DOUBLE PRECISION.
This may cause a loss of precision.
Shouldn't be TmpValues declared DOUBLE PRECISION as well?