def test_Pq(): ''' ''' from eos import Nominal, Spline_eos vt = (data()) eos = Spline_eos(Nominal(), precondition=False) c = eos.get_c() gun = Gun(eos) P, q = eos.Pq_like(*gun.compare(vt, c)) assert P.shape == (50, 50) assert q.shape == (50, ) assert P[23, 23] == P.max() value = 3.286529435e-10 assert close(P.max(), value), 'P.max={0:.9e} != {1:.9e}'.format(P.max(), value) i = 22 assert np.argmin(q) == i, 'argmin={0} != 11'.format(np.argmin(q), i) value = 5.680863201e-01 assert close(-q[i], value), '-q[{2}]={0:.9e} != {1:.9e}'.format(-q[i], value, i) return 0