- HXb = Hm * Xb[iens]
- if "SimulatedObservationAtBackground" in self._parameters["StoreSupplementaryCalculations"]:
- self.StoredVariables["SimulatedObservationAtBackground"].store( numpy.ravel(HXb) )
- d = EnsembleY[:,iens] - HXb
- if "Innovation" in self._parameters["StoreSupplementaryCalculations"]:
- self.StoredVariables["Innovation"].store( numpy.ravel(d) )
- Xa = Xb[iens] + K*d
+ HXb = numpy.ravel(numpy.dot(Hm, Xb[iens]))
+ if self._toStore("SimulatedObservationAtBackground"):
+ self.StoredVariables["SimulatedObservationAtBackground"].store( HXb )
+ d = numpy.ravel(EnsembleY[:,iens]) - HXb
+ if self._toStore("Innovation"):
+ self.StoredVariables["Innovation"].store( d )
+ Xa = numpy.ravel(Xb[iens]) + numpy.dot(K, d)