From 2c515cc5c9ab8ee3e6495951ab4fcccfdc63d275 Mon Sep 17 00:00:00 2001
From: ivar <amvdw95@gmail.com>
Date: Thu, 2 Aug 2012 06:46:32 +0000
Subject: [PATCH] rejected samples are now correctly flagged with number 2
---
da/baseclasses/optimizer.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/da/baseclasses/optimizer.py b/da/baseclasses/optimizer.py
index 1b856211..38a6bdea 100755
--- a/da/baseclasses/optimizer.py
+++ b/da/baseclasses/optimizer.py
@@ -132,7 +132,7 @@ class Optimizer(object):
self.obs[:] = np.array(allobs)
self.obs_ids[:] = np.array(allids)
-
+ print self.HX_prime[:,:].shape,np.array(allsimulated).shape
self.HX_prime[:,:] = np.array(allsimulated)
self.Hx[:] = self.HX_prime[:,0]
@@ -354,7 +354,7 @@ class Optimizer(object):
threshold = self.rejection_threshold*np.sqrt(self.R[n,n])
if np.abs(res) > threshold:
msg = 'Rejecting observation %d because residual (%f) exceeds threshold (%f)'%(n,res,threshold) ; logging.debug(msg)
- self.flags[n] == 2
+ self.flags[n] = 2
continue
PHt = 1./(self.nmembers-1)*np.dot(self.X_prime,self.HX_prime[n,:])
--
GitLab