auxiliary.h

Go to the documentation of this file.

#ifndef AUXILIARY_H_
#define AUXILIARY_H_

#include <complex.h>
#include <stdio.h>
#include <stdbool.h>
#include <fftw3.h>

//constants
#define PI 3.1415926535897932384626433832795029L

#define TWO_PI (2.0 * PI)

#define FOUR_PI (4.0 * PI)

#define MU0 1.256637061435917e-6

#define EPS0 8.854187817620e-12

#define RHO_CU 1.689e-8


double*
linspace (double a, double b, size_t n, double u[]);

double*
logspace (double a, double b, size_t n, double u[]);

double
wave_length (double f, double sigma, double ep, double mur);

int
alipio_soil (_Complex double* sigma, _Complex double* epsr, double sigma0,
               _Complex double s, double h, double g, double eps_ratio);

/* Calculates the soil parameters \f$ \sigma(s) \f$ and \f$ \varepsilon_r(s) \f$
based on the Smith-Longmire model as presented in [1].

[1] D. Cavka, N. Mora, F. Rachidi, A comparison of frequency-dependent soil
models: application to the analysis of grounding systems, IEEE Trans.
Electromagn. Compat. 56 (February (1)) (2014) 177–187.

@param sigma pointer to where the conductivity \$f \sigma(s) \$f is written in S/m
@param epsr pointer to where the relative permitivitty \f$ \varepsilon_r(s) \f$ is written
@param sigma0 value of the soil conductivity in low frequency \f$ \sigma_0 \f$ in S/m
@param s complex frequency \f$ s = c + j\omega \f$ of interest in rad/s
@param erinf parameter \f$ \varepsilon_\infty'\f$
@return 0 on success
*/
int
smith_longmire_soil (_Complex double* sigma, _Complex double* epsr, double sigma0,
                     _Complex double s, double erinf);

bool
equal_points (const double* point_1, const double* point_2);

bool
equal_points_tol (const double* point_1, const double* point_2, double tol);

double
vector_length (const double start_point[3], const double end_point[3]);

extern double
snrm2_ (int* n, double* x, int* incx);

int
complex_matrix_file (size_t m, size_t n, const _Complex double* a, int lda, FILE* fp);

int
double_matrix_file (size_t m, size_t n, const double* a, int lda, FILE* fp);

extern int
zbesi_ (double* zr, double* zi, double* fnu, int* kode, int* n, double* cyr,
        double* cyi, int* nz, int* ierr);

void
print_matrix (char *desc, int m, int n, const _Complex double* a, int lda);

void
print_matrix_row (char *desc, int m, int n, const _Complex double* a, int lda);

int
matrix_copy (const _Complex double* source, _Complex double* target,
                   size_t lds, size_t ldt, size_t nline, size_t ncol);

int
transpose_copy (const _Complex double* source, _Complex double* target,
                size_t lds, size_t ldt, size_t nline, size_t ncol);

int
pc_copy (const _Complex double* source, _Complex double* target,
         size_t lds, size_t ldt, size_t nline, size_t ncol);

int
pl_copy (const _Complex double* source, _Complex double* target,
         size_t lds, size_t ldt, size_t nline, size_t ncol);

int
pcl_copy (const _Complex double* source, _Complex double* target,
          size_t lds, size_t ldt, size_t nline, size_t ncol);


int
laplace_trans (double* f, _Complex double* g, _Complex double* s, double tmax, size_t nt);

enum
INLT_Filter {
        FILTER_NONE = 0,
      FILTER_BLACKMAN = 1,
    FILTER_HANNING = 2,
    FILTER_LANCZOS = 3,
    FILTER_RIESZ = 4
};

int
inv_laplace_trans (double* f, _Complex double* g, _Complex double* s, double tmax,
                     size_t nt, int filter);

double
heidler (double t, double imax, double tau1, double tau2, int n);

#endif /* AUXILIARY_H_ */