6. Integration module

This module actually contains the Heun algorithm routines.

Note

The python code is available here : integrator.py .

Example:
>>> from integrator import step
>>> step(y,t,dt)

6.1. Global variables

  • aotensor tensor with the format (int i, int j, int k, float v) in list
  • Li first list of index of tensor
  • Lj second list of index of tensor
  • Lk third list of index of tensor
  • Lv list of tensor values

6.2. Dependencies

>>> import numpy as np
>>> from params_maooam import ndim,f2py
>>> import aotensor as aotensor_mod
>>> if f2py:
>>>     import sparse_mult as mult
>>>     sparse_mul3_f2py = mult.sparse_mult.sparse_mul3

6.3. Functions

  • sparse_mul3
  • tendencies
  • step
integrator.sparse_mul3(arr)

Calculate for each i the sums on j,k of the product

\[tensor(i,j,k)* arr(j) * arr(k)\]

Note

Python-only function

integrator.sparse_mul3_py(arr)[source]

Calculate for each i the sums on j,k of the product

\[tensor(i,j,k)* arr(j) * arr(k)\]

Note

Python-only function

integrator.step(y, t, dt)[source]

RK2 method integration

integrator.tendencies(y)[source]

Calculate the tendencies thanks to the product of the tensor and the vector y