proteus
proteus copied to clipboard
possible bug in Advection_ASGS df_last initialization
Aron and I found this when reviewing the code. Looks like the ('df',ci,ci) key is missing
https://github.com/erdc-cm/proteus/blob/master/src/SubgridError.py#L76
for ci in range(self.nc):
if self.lag:
self.tau_last.append(numpy.zeros(cq[('u',ci)].shape,'d'))
self.tau.append(numpy.zeros(cq[('u',ci)].shape,'d'))
if cq.has_key(('df',ci,ci)):
self.df_last = copy.deepcopy(cq[('df',ci,ci)])
cq[('df_sge',ci,ci)] = self.df_last
OK, meaning at this stage cq doesn't have ('df',ci,ci) yet so we don't get a proper deepcopy?
@cekees I can't find the file reference. Closing this. Probably the file has been moved (change of address). Close and reopen?