xref: /libCEED/python/tests/test-2-elemrestriction.py (revision 0ef725981a32b9079ff6c5100673b913b8f4d7c0)

# Copyright (c) 2017, Lawrence Livermore National Security, LLC. Produced at
# the Lawrence Livermore National Laboratory. LLNL-CODE-734707. All Rights
# reserved. See files LICENSE and NOTICE for details.
#
# This file is part of CEED, a collection of benchmarks, miniapps, software
# libraries and APIs for efficient high-order finite element and spectral
# element discretizations for exascale applications. For more information and
# source code availability see http://github.com/ceed.
#
# The CEED research is supported by the Exascale Computing Project 17-SC-20-SC,
# a collaborative effort of two U.S. Department of Energy organizations (Office
# of Science and the National Nuclear Security Administration) responsible for
# the planning and preparation of a capable exascale ecosystem, including
# software, applications, hardware, advanced system engineering and early
# testbed platforms, in support of the nation's exascale computing imperative.

# @file
# Test Ceed ElemRestriction functionality

import os
import libceed
import numpy as np
import check

# -------------------------------------------------------------------------------
# Test creation, use, and destruction of an element restriction
# -------------------------------------------------------------------------------


def test_200(ceed_resource):
    ceed = libceed.Ceed(ceed_resource)

    ne = 3

    x = ceed.Vector(ne + 1)
    a = np.arange(10, 10 + ne + 1, dtype="float64")
    x.set_array(a, cmode=libceed.USE_POINTER)

    ind = np.zeros(2 * ne, dtype="int32")
    for i in range(ne):
        ind[2 * i + 0] = i
        ind[2 * i + 1] = i + 1
    r = ceed.ElemRestriction(ne, 2, 1, 1, ne + 1, ind,
                             cmode=libceed.USE_POINTER)

    y = ceed.Vector(2 * ne)
    y.set_value(0)

    r.apply(x, y)

    with y.array_read() as y_array:
        for i in range(2 * ne):
            assert 10 + (i + 1) // 2 == y_array[i]

# -------------------------------------------------------------------------------
# Test creation, use, and destruction of a strided element restriction
# -------------------------------------------------------------------------------


def test_201(ceed_resource):
    ceed = libceed.Ceed(ceed_resource)

    ne = 3

    x = ceed.Vector(2 * ne)
    a = np.arange(10, 10 + 2 * ne, dtype="float64")
    x.set_array(a, cmode=libceed.USE_POINTER)

    strides = np.array([1, 2, 2], dtype="int32")
    r = ceed.StridedElemRestriction(ne, 2, 1, 2 * ne, strides)

    y = ceed.Vector(2 * ne)
    y.set_value(0)

    r.apply(x, y)

    with y.array_read() as y_array:
        for i in range(2 * ne):
            assert 10 + i == y_array[i]

# -------------------------------------------------------------------------------
# Test creation and destruction of a blocked element restriction
# -------------------------------------------------------------------------------


def test_202(ceed_resource, capsys):
    ceed = libceed.Ceed(ceed_resource)

    ne = 8
    blksize = 5

    x = ceed.Vector(ne + 1)
    a = np.arange(10, 10 + ne + 1, dtype="float64")
    x.set_array(a, cmode=libceed.USE_POINTER)

    ind = np.zeros(2 * ne, dtype="int32")
    for i in range(ne):
        ind[2 * i + 0] = i
        ind[2 * i + 1] = i + 1
    r = ceed.BlockedElemRestriction(ne, 2, blksize, 1, 1, ne + 1, ind,
                                    cmode=libceed.USE_POINTER)

    y = ceed.Vector(2 * blksize * 2)
    y.set_value(0)

    r.apply(x, y)

    print(y)

    x.set_value(0)
    r.T.apply(y, x)
    print(x)

    stdout, stderr, ref_stdout = check.output(capsys)
    assert not stderr
    assert stdout == ref_stdout

# -------------------------------------------------------------------------------
# Test creation, use, and destruction of a blocked element restriction
# -------------------------------------------------------------------------------


def test_208(ceed_resource, capsys):
    ceed = libceed.Ceed(ceed_resource)

    ne = 8
    blksize = 5

    x = ceed.Vector(ne + 1)
    a = np.arange(10, 10 + ne + 1, dtype="float64")
    x.set_array(a, cmode=libceed.USE_POINTER)

    ind = np.zeros(2 * ne, dtype="int32")
    for i in range(ne):
        ind[2 * i + 0] = i
        ind[2 * i + 1] = i + 1
    r = ceed.BlockedElemRestriction(ne, 2, blksize, 1, 1, ne + 1, ind,
                                    cmode=libceed.USE_POINTER)

    y = ceed.Vector(blksize * 2)
    y.set_value(0)

    r.apply_block(1, x, y)

    print(y)

    x.set_value(0)
    r.T.apply_block(1, y, x)
    print(x)

    stdout, stderr, ref_stdout = check.output(capsys)
    assert not stderr
    assert stdout == ref_stdout

# -------------------------------------------------------------------------------
# Test getting the multiplicity of the indices in an element restriction
# -------------------------------------------------------------------------------


def test_209(ceed_resource):
    ceed = libceed.Ceed(ceed_resource)

    ne = 3

    ind = np.zeros(4 * ne, dtype="int32")
    for i in range(ne):
        ind[4 * i + 0] = i * 3 + 0
        ind[4 * i + 1] = i * 3 + 1
        ind[4 * i + 2] = i * 3 + 2
        ind[4 * i + 3] = i * 3 + 3
    r = ceed.ElemRestriction(ne, 4, 1, 1, 3 * ne + 1, ind,
                             cmode=libceed.USE_POINTER)

    mult = r.get_multiplicity()

    with mult.array_read() as mult_array:
        for i in range(3 * ne + 1):
            val = 1 + (1 if (i > 0 and i < 3 * ne and i % 3 == 0) else 0)
            assert val == mult_array[i]

# -------------------------------------------------------------------------------
# Test creation and view of an element restriction
# -------------------------------------------------------------------------------


def test_210(ceed_resource, capsys):
    ceed = libceed.Ceed(ceed_resource)

    ne = 3

    ind = np.zeros(2 * ne, dtype="int32")
    for i in range(ne):
        ind[2 * i + 0] = i + 0
        ind[2 * i + 1] = i + 1
    r = ceed.ElemRestriction(ne, 2, 1, 1, ne + 1, ind,
                             cmode=libceed.USE_POINTER)

    print(r)

    stdout, stderr, ref_stdout = check.output(capsys)
    assert not stderr
    assert stdout == ref_stdout

# -------------------------------------------------------------------------------
# Test creation and view of a strided element restriction
# -------------------------------------------------------------------------------


def test_211(ceed_resource, capsys):
    ceed = libceed.Ceed(ceed_resource)

    ne = 3

    strides = np.array([1, 2, 2], dtype="int32")
    r = ceed.StridedElemRestriction(ne, 2, 1, ne + 1, strides)

    print(r)

    stdout, stderr, ref_stdout = check.output(capsys)
    assert not stderr
    assert stdout == ref_stdout

# -------------------------------------------------------------------------------
# Test creation and view of a blocked strided element restriction
# -------------------------------------------------------------------------------


def test_212(ceed_resource, capsys):
    ceed = libceed.Ceed(ceed_resource)

    ne = 3

    strides = np.array([1, 2, 2], dtype="int32")
    r = ceed.BlockedStridedElemRestriction(ne, 2, 2, 1, ne + 1, strides)

    print(r)

    stdout, stderr, ref_stdout = check.output(capsys)
    assert not stderr
    assert stdout == ref_stdout

# -------------------------------------------------------------------------------