def test_bq2(): "2-point boundary quadrature" raise SkipTest("not implemented") dx = 1.0 dy = dx Q = PISM.BoundaryQuadrature2(1, 1, 1) for q in [0, 1]: for s in [0, 1, 2, 3]: psi = [Q.germ(s, q, k).val for k in [0, 1]] assert sum(psi) == 1.0, "side = {}, q = {}, psi = {}".format( s, q, psi)