def test_NonLinearDS_suboptimalSysID():
"""Test NonLinearDS system identification.
"""
dataFile = os.path.join(TESTBASE, "data/data1.txt")
data, _ = loadDataFromASCIIFile(dataFile)
kpcaP = KPCAParam()
kpcaP._kPar = RBFParam()
kpcaP._kPar._kCen = True
kpcaP._kFun = rbfK
nlds = NonLinearDS(5, kpcaP, False)
nlds.suboptimalSysID(data)
baseNLDSFile = os.path.join(TESTBASE, "data/data1-rbf-kdt-5c-center.pkl")
baseNLDS = pickle.load(open(baseNLDSFile))
err = NonLinearDS.naiveCompare(baseNLDS, nlds)
assert np.allclose(err, 0.0) == True
def test_NonLinearDS_suboptimalSysID():
"""Test NonLinearDS system identification.
"""
dataFile = os.path.join(TESTBASE, "data/data1.txt")
data, _ = loadDataFromASCIIFile(dataFile)
kpcaP = KPCAParam()
kpcaP._kPar = RBFParam()
kpcaP._kPar._kCen = True
kpcaP._kFun = rbfK
nlds = NonLinearDS(5, kpcaP, False)
nlds.suboptimalSysID(data)
baseNLDSFile = os.path.join(TESTBASE, "data/data1-rbf-kdt-5c-center.pkl")
baseNLDS = pickle.load(open(baseNLDSFile))
err = NonLinearDS.naiveCompare(baseNLDS, nlds)
np.testing.assert_almost_equal(err, 0, 2)
