sempy 0.0.18 (experimental)

A Python implementation of the spectral element method

Getting started

The quickest way to get started using the Sempy package is to install it with easy_install This can be achieved by e.g. typing:

$ easy_install sempy

on the command line. This requires an internet connection to work. Another possibility is to download the latest version as a zip or tar.gz file and then install it with:

$ python setup.py install

It might be necessary to use super user privileges on the above commands. For instance:

$ sudo easy_install sempy

Aslo, it might be necessary to run:

$ sudo easy_install --upgrade sempy

on a regular basis to obtain updated versions.

An example script

Solve the Poisson problem with a simple script:

import sempy

X = sempy.Space( filename = 'square', n = 4, dim = 2 )
A = sempy.operators.Laplacian( X ).matrix

f = sempy.Function( X, basis_coeff = 1.0 )
b = sempy.operators.Mass( X ).action_local( f.basis_coeff )
u = sempy.Function( X, basis_coeff = 0.0 )

[v, flag] = sempy.linsolvers.Krylov().solve( A, b, u.glob() )
u.basis_coeff = X.mapping_q( v )
u.plot()