Exemple #1
0
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