Closed jotaraz closed 6 months ago
Well, Julia's scoping rules are different from C/C++...
This code is tested in Example430. It worked for two reasons:
once idofL is set to some value, it is defined also outside the if/end clause, as this does not create a scope block
The example is 1D has one species, so there is only one possible value of idofL
... need to add a multispecies test here, so I keep this open...
In two functions in src/vfvm_assemblydata.jl the
iparam
-for-loop to be moved into some other loops such that all used variables are defined.line 337:
@inline function assemble_res_jac(edge...
theiparam
-for-loop has to be inside theisnodespecies(system, ispec, K)
if-statement, otherwiseidofL
is not defined.line 400:
@inline function assemble_res_jac(bedge...
theiparam
-for-loop has to be insideidofK
for-loop and theisnodespecies(system, ispec, K)
if-statement, otherwiseidofK
andidofL
are not defined.In the first function: