FESTIM icon indicating copy to clipboard operation
FESTIM copied to clipboard

Derived quantity objects `data` and `t` attributes are not emptied at beginning of simulation

Open RemDelaporteMathurin opened this issue 1 year ago • 1 comments

Describe the bug When running the same simulation twice with derived quantities, their t and data attributes isn't reset to an empty list but instead values are just appended.

To Reproduce Run the following on fenicsx

import festim as F

my_model = F.HydrogenTransportProblem()

my_model.mesh = F.Mesh1D([0, 1, 2, 3, 4])

tungsten = F.Material(D_0=1, E_D=0)

my_subdomain = F.VolumeSubdomain1D(id=1, borders=[0, 4], material=tungsten)
left_surface = F.SurfaceSubdomain1D(id=1, x=0)

my_model.subdomains = [my_subdomain, left_surface]

mobile_H = F.Species("H")
my_model.species = [mobile_H]

my_model.temperature = 400

left_flux = F.SurfaceFlux(field=mobile_H, surface=left_surface)

my_model.exports = [left_flux]

my_model.settings = F.Settings(
    atol=1e-15,
    rtol=1e-15,
    transient=False,
)

for i in range(3):
    my_model.initialise()
    my_model.run()

    print(left_flux.data)

Produces:

[0.0]
[0.0, 0.0]
[0.0, 0.0, 0.0]

Expected behavior left_flux.data should be of length 1 everytime

RemDelaporteMathurin avatar Apr 01 '24 13:04 RemDelaporteMathurin

I noticed this bug also occurs on the main branch!

RemDelaporteMathurin avatar Apr 05 '24 17:04 RemDelaporteMathurin