jhale

def lagrange_points(cell, degree):
    """Construct the locations of the equispaced Lagrange nodes on cell.

    :param cell: the :class:`~.reference_elements.ReferenceCell`
    :param degree: the degree of polynomials for which to construct nodes.

    :returns: a rank 2 :class:`~numpy.array` whose rows are the
        coordinates of the nodes.

    The implementation of this function is left as an :ref:`exercise

jhale

Spack for FEniCS on University of Luxembourg clusters

This should also be useful for other clusters where you want to use components (e.g. MPI, compilers) from the module system.

Setting up Spack

Start a session for building

si -N 1 -n 16 -c 1 -t 0-02:00:00 # on iris: -C broadwell or -C skylake

jhale

import numpy as np
import scipy as sp

np.set_printoptions(precision=15)


def maclaurin_series(n):
    m2pi = -2 * np.pi
    ns = np.arange(0, n)
    xs = np.power(m2pi, ns) / sp.special.factorial(ns)

jhale

jhale

jhale

jhale

jhale

jhale

import sys

# Here I make two namespaces, ns_ and A_.
# You can make as many as you want, probably
# one for every KSP/SVD/EPS in your code.
args = [sys.argv[0]] + """
--petsc.ns_svd_type lapack
--petsc.A_mat_type aij
""".split()

jhale

.vagrant/