File: compositebasistest.py

package info (click to toggle)
dune-functions 2.10.0-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 1,544 kB
  • sloc: cpp: 14,241; python: 661; makefile: 3
file content (35 lines) | stat: -rw-r--r-- 1,447 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
# SPDX-FileCopyrightText: Copyright © DUNE Project contributors, see file AUTHORS.md
# SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later

import dune.grid
import dune.functions as functions

import basistest   # The general test suite for function space bases


# Test all currently supported global bases in dimension `dimension`.
# Only the size of the global bases are tested.
def test(dimension):
    lowerLeft = [0] * dimension
    upperRight = [1] * dimension
    elements = [4] * dimension

    grid = dune.grid.structuredGrid(lowerLeft,upperRight,elements)

    # Test a basis consisting of three identical scalar-valued children
    basis = functions.defaultGlobalBasis(grid,
                                         functions.Composite(functions.Lagrange(order=1),
                                                             functions.Lagrange(order=1),
                                                             functions.Lagrange(order=1)))
    basistest.checkBasis(basis)

    # Test a composite of a power node and a Lagrange basis
    basis = functions.defaultGlobalBasis(grid,
                                         functions.Composite(functions.Power(functions.Lagrange(order=2),exponent=dimension),
                                                             functions.Lagrange(order=1)))
    basistest.checkBasis(basis)


# Run tests for grids of dimension 2 and 3
test(2)
test(3)