def atest_monte_carlo_metropolis_hasting_sampling_with_duplicates(self):
self.set_positions_resampling_set()
setNpRandomSeed()
indices = [0,0,0,3,4,4]
solutions = self.resample(indices)
dautils.metropolisHastingSampling(self.resamplingDrifterSet)
self.assertEqual(self.resamplingDrifterSet.getDrifterPositions().tolist(), solutions)
def atest_monte_carlo_metropolis_hasting_sampling_with_duplicates(self):
self.set_positions_resampling_set()
setNpRandomSeed()
indices = [0,0,0,3,4,4]
solutions = self.resample(indices)
dautils.metropolisHastingSampling(self.resamplingDrifterSet)
self.assertEqual(self.resamplingDrifterSet.getDrifterPositions().tolist(), solutions)
def atest_monte_carlo_metropolis_hasting_sampling(self):
self.set_positions_resampling_set()
setNpRandomSeed()
indices = [0,0,0,3,4,4]
solutions = self.resample(indices)
dautils.metropolisHastingSampling(self.resamplingDrifterSet, self.resamplingVar)
#print self.resamplingDrifterSet.getDrifterPositions().tolist()
assert2DListAlmostEqual(self, self.resamplingDrifterSet.getDrifterPositions().tolist(), solutions, 2, "metropolis hasting sampling, probabilistic duplicates")
def atest_monte_carlo_metropolis_hasting_sampling(self):
self.set_positions_resampling_set()
setNpRandomSeed()
indices = [0,0,0,3,4,4]
solutions = self.resample(indices)
dautils.metropolisHastingSampling(self.resamplingDrifterSet, self.resamplingVar)
#print self.resamplingDrifterSet.getDrifterPositions().tolist()
assert2DListAlmostEqual(self, self.resamplingDrifterSet.getDrifterPositions().tolist(), solutions, 2, "metropolis hasting sampling, probabilistic duplicates")
