Python main_integrate3dの例

コード例 #1

def test_main_integrate3d():
    cube = [
        [
            (0, 0, 0),
            (0, 0, 5),
            (0, 5, 0),
            (0, 5, 5),
            (5, 0, 0),
            (5, 0, 5),
            (5, 5, 0),
            (5, 5, 5),
        ],
        [2, 6, 7, 3],
        [3, 7, 5, 1],
        [7, 6, 4, 5],
        [1, 5, 4, 0],
        [3, 1, 0, 2],
        [0, 4, 6, 2],
    ]
    vertices = cube[0]
    faces = cube[1:]
    hp_params = hyperplane_parameters(faces, vertices)
    assert main_integrate3d(1, faces, vertices, hp_params) == -125
    assert main_integrate3d(1, faces, vertices, hp_params, max_degree=1) == {
        1: -125,
        y: Rational(-625, 2),
        z: Rational(-625, 2),
        x: Rational(-625, 2),
    }

コード例 #2

ファイル: test_intpoly.py プロジェクト: asmeurer/sympy

def test_main_integrate3d():
    cube = [[(0, 0, 0), (0, 0, 5), (0, 5, 0), (0, 5, 5), (5, 0, 0),\
             (5, 0, 5), (5, 5, 0), (5, 5, 5)],\
            [2, 6, 7, 3], [3, 7, 5, 1], [7, 6, 4, 5], [1, 5, 4, 0],\
            [3, 1, 0, 2], [0, 4, 6, 2]]
    vertices = cube[0]
    faces = cube[1:]
    hp_params = hyperplane_parameters(faces, vertices)
    assert main_integrate3d(1, faces, vertices, hp_params) == -125
    assert main_integrate3d(1, faces, vertices, hp_params, max_degree=1) == \
        {1: -125, y: -S(625)/2, z: -S(625)/2, x: -S(625)/2}

コード例 #3

def test_main_integrate3d():
    cube = [[(0, 0, 0), (0, 0, 5), (0, 5, 0), (0, 5, 5), (5, 0, 0),\
             (5, 0, 5), (5, 5, 0), (5, 5, 5)],\
            [2, 6, 7, 3], [3, 7, 5, 1], [7, 6, 4, 5], [1, 5, 4, 0],\
            [3, 1, 0, 2], [0, 4, 6, 2]]
    vertices = cube[0]
    faces = cube[1:]
    hp_params = hyperplane_parameters(faces, vertices)
    assert main_integrate3d(1, faces, vertices, hp_params) == -125
    assert main_integrate3d(1, faces, vertices, hp_params, max_degree=1) == \
        {1: -125, y: -S(625)/2, z: -S(625)/2, x: -S(625)/2}