<reponame>TareqZaman/pyamg
"""Test C/F splitting and simple aggregation vis."""

import tempfile
import os

from scipy.sparse import csr_matrix
from numpy import array, ones, uint32

from numpy.testing import TestCase

from pyamg.vis import vis_aggregate_groups


class TestVis(TestCase):
    def setUp(self):
        self.file_name = tempfile.mktemp()

    def tearDown(self):
        os.remove(self.file_name)

    def test_1(self):
        V = array([[0.0, 0.0],
                   [1.0, 0.0],
                   [2.0, 0.0],
                   [0.0, 1.0],
                   [1.0, 1.0],
                   [2.0, 1.0],
                   [0.0, 2.0],
                   [1.0, 2.0],
                   [2.0, 2.0],
                   [0.0, 3.0],
                   [1.0, 3.0],
                   [2.0, 3.0]])

        E2V = array([[0, 4, 3],
                     [0, 1, 4],
                     [1, 5, 4],
                     [1, 2, 5],
                     [3, 7, 6],
                     [3, 4, 7],
                     [4, 8, 7],
                     [4, 5, 8],
                     [6, 10, 9],
                     [6, 7, 10],
                     [7, 11, 10],
                     [7, 8, 11]], dtype=uint32)

        row = array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11])
        col = array([1, 0, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1])
        data = ones((1, 12), dtype=uint32).ravel()

        AggOp = csr_matrix((data, (row, col)), shape=(12, 2))

        vis_aggregate_groups(V=V, E2V=E2V, AggOp=AggOp, mesh_type='tri',
                             fname=self.file_name)

    def test_2(self):
        V = array([[0.0, 0.0],
                   [1.0, 0.0],
                   [2.0, 0.0],
                   [0.0, 1.0],
                   [1.0, 1.0],
                   [2.0, 1.0],
                   [3.0, 1.0],
                   [0.0, 2.0],
                   [1.0, 2.0],
                   [2.0, 2.0],
                   [3.0, 2.0],
                   [4.0, 2.0],
                   [0.0, 3.0],
                   [1.0, 3.0],
                   [2.0, 3.0],
                   [3.0, 3.0],
                   [4.0, 3.0],
                   [5.0, 3.0]])

        E2V = array([[0, 4, 3],
                     [0, 1, 4],
                     [1, 5, 4],
                     [1, 2, 5],
                     [2, 6, 5],
                     [3, 8, 7],
                     [3, 4, 8],
                     [4, 9, 8],
                     [4, 5, 9],
                     [5, 10, 9],
                     [5, 6, 10],
                     [6, 11, 10],
                     [7, 13, 12],
                     [7, 8, 13],
                     [8, 14, 13],
                     [8, 9, 14],
                     [9, 15, 14],
                     [9, 10, 15],
                     [10, 16, 15],
                     [10, 11, 16],
                     [11, 17, 16]], dtype=uint32)

        row = array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
                     10, 11, 12, 13, 14, 15, 16, 17])
        col = array([0, 1, 3, 0, 1, 1, 3, 0, 0,
                     1, 3, 4, 0, 0, 0, 2, 4, 4])
        data = ones((1, 18), dtype=uint32).ravel()

        AggOp = csr_matrix((data, (row, col)), shape=(18, 5))

        vis_aggregate_groups(V=V, E2V=E2V, AggOp=AggOp, mesh_type='tri',
                             fname=self.file_name)
