linalg.h

Go to the documentation of this file.

/* High Performance implementation of the Hybrid Electromagnetic Model
Released under the General Public License 3 (GPLv3).
All parameters' units are in the SI base if omitted.

Routines to build and solve the electrode system that depend on linear
algebra libraries (e.g. LAPACK and BLAS).

Every matrix is stored as a flat array in column major format.
*/
#ifndef LINALG_H_
#define LINALG_H_

#include <complex.h>
#include <string.h>
#include "electrode.h"
#include "blas.h"
#include "lapack.h"
//#include "mkl.h"
//#include "mkl_lapacke.h"

/* ================= Immittance Formulation ================= */
int
fill_incidence_imm (_Complex double* wg, const Electrode* electrodes,
                    size_t num_electrodes, const double* nodes,
                    size_t num_nodes);

int
fill_impedance_imm (_Complex double* wg, const _Complex double* zl,
                    const _Complex double* zt, size_t num_electrodes,
                    size_t num_nodes);

int
solve_immittance (_Complex double* wg, _Complex double* ie, size_t num_electrodes,
                  size_t num_nodes);

/* ================= Admittance Formulation ================= */
int
fill_incidence_adm (_Complex double* a, _Complex double* b,
                    const Electrode* electrodes, size_t num_electrodes,
                    const double* nodes, size_t num_nodes);

int
fill_impedance_adm (_Complex double* yn, _Complex double* zl, _Complex double* zt,
                    _Complex double* a, _Complex double* b, size_t num_electrodes,
                    size_t num_nodes);

int
calculate_yla_ytb (_Complex double* yla, _Complex double* ytb,
                   _Complex double* zl, _Complex double* zt,
                   _Complex double* a, _Complex double* b,
                   size_t num_electrodes, size_t num_nodes);

int
fill_impedance_adm2 (_Complex double* yn, _Complex double* yla,
                     _Complex double* ytb, _Complex double* a,
                     _Complex double* b, size_t num_electrodes,
                     size_t num_nodes);

int
solve_admittance (_Complex double* yn, _Complex double* ic, size_t num_nodes);

#endif /* LINALG_H_ */