lin_alg.h

Go to the documentation of this file.

#ifndef LIN_ALG_H
# define LIN_ALG_H


# ifdef __cplusplus
extern "C" {
# endif // ifdef __cplusplus

# include "types.h"
#include "solver_interface.h"

c_float* vec_copy(const c_float *a,
                  size_t         n);

void prea_vec_copy(const c_float *a,
                   c_float       *b,
                   size_t         n);

void prea_int_vec_copy(const c_int *a,
                       c_int       *b,
                       size_t       n);

void vec_set_scalar(c_float *a,
                    c_float  sc,
                    size_t   n);

void vec_set_scalar_int(c_int *a, 
                        c_int  sc, 
                        size_t n);

void vec_self_mult_scalar(  c_float *a,
                            c_float sc,
                            size_t  n);

void vec_mult_scalar(const c_float  *a,
                     c_float        sc,
                     c_float        *b,
                     size_t         n);

c_float vec_prod(const c_float *a,
                 const c_float *b,
                 size_t         n);

c_float vec_norm_two(   const c_float   *a, 
                        size_t          n);

c_float vec_norm_inf(const c_float *a,
                     size_t         n);

void vec_add_scaled(const c_float *a,
                    const c_float *b,
                    c_float       *c,
                    c_float        sc,
                    size_t         n);

void vec_mult_add_scaled(   c_float         *a, 
                            const c_float   *b, 
                            c_float         sc1, 
                            c_float         sc2, 
                            size_t          n);

void vec_ew_recipr(const c_float *a,
                   c_float       *b,
                   size_t         n);

void vec_ew_max_vec(const c_float *a,
                    const c_float *b,
                    c_float       *c,
                    size_t         n);

void vec_ew_min_vec(const c_float *a,
                    const c_float *b,
                    c_float       *c,
                    size_t         n);

void vec_ew_mid_vec(const c_float *a,
                    const c_float *bmin,
                    const c_float *bmax,
                    c_float       *c,
                    size_t         n);

void vec_ew_prod(const c_float *a,
                 const c_float *b,
                 c_float       *c,
                 size_t         n);

void vec_ew_div(const c_float *a,
                const c_float *b,
                c_float       *c,
                size_t         n);

void vec_ew_sqrt(const c_float *a,
                 c_float       *b,
                 size_t         n);


/* MATRIX FUNCTIONS ----------------------------------------------------------*/

/* Moved to solver_interface.c*/


# ifdef __cplusplus
}
# endif // ifdef __cplusplus

#endif // ifndef LIN_ALG_H