vector.h

Go to the documentation of this file.

/* -----------------------------------------------------------------------------------------
  $Id: vector.h 571 2006-11-10 18:30:25Z svacek $
  $Author: svacek $
   ----------------------------------------------------------------------------------------- */
#include <malloc.h>
#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include "wlog.h"
#ifndef __VECTOR_SVH
#define __VECTOR_SVH
/* --------------------------------------------------------------------------------  */
/*@{*/
/** \defgroup VECTOR_OPS Vector Operations */
/** \ingroup LINEAR_ALGEBRA */
/* --------------------------------------------------------------------------------  */
typedef struct
{
  double *a;
  long n;
} vector; /**< Vector of n-th DOUBLEs */

void set_vec(vector *pom,double *a,long n);                                     /**< set link to memory to treat that as vector */
void link_vec(vector *pom,const vector u);                                      /**< set link to memory of the vector u */
void concat2_vec(vector *uv,const vector u,const vector v);                     /**< Concatenate 2 vectors */
void concat3_vec(vector *uv,const vector u,const vector v,const vector w);      /**< Concatenate 3 vectors */
void construct_vec(vector *u,long k);                                           /**< Creates vector u */
void construct_vec_copy(vector *u,const vector v);                              /**< Creates copy of the vector v */
void construct_vec_concat2(vector *uv,const vector u,const vector v);           /**< Construct Concatenation of  2 vectors */
void construct_vec_concat3(vector *uv,const vector u,const vector v,const vector w);  /**< Construct Concatenation of 3 vectors */
void destruct_vec(vector *which);                                               /**< Destruct the vector */
void print_vec(const vector u,char *name);                                      /**< Destruct the vector */
void split_vec(vector *out,const vector in,long start,long end);                /**< Take one part start..end of vector IN and puts it into (alloceted) OUT */
/**< Computes the norm of vector U :
     - l1-norm (k=1), 
     - l2-norm(k=2),
     - l-\f$\infty\f$ norm(k=3) ) 
*/
double norm_vec(const vector u,int k);
/**< Computes the distance of vector u and v as \f$\|u-v\|_k\f$ :
     - l1-norm (k=1), 
     - l2-norm(k=2),
     - l-\f$\infty\f$ norm(k=3) ) 
*/
double distance_vec(const vector u,const vector v,int k);
void zeroize_vec(vector u);                                                     /**< Puts zero into vector */
void load_vec(vector *u,const char *FileName);                                  /**< Load vector from (ASCII) file */
void load_vec_bin(vector *u,const char *FileName);                              /**< Load vector from (BINARY) file */
void save_vec(const vector u,const char *FileName);                             /**< Saves vector to (ASCII) file */
void save_vec_bin(const vector u,const char *FileName);                         /**< Saves vector to (BINARY) file */
void add_vec(vector u,const vector v);                                          /**< Makes u := u+v */
void add_vec_val(vector u,double alpha);                                        /**< Makes u := u+alpha */
void subtract_vec(vector u,const vector v);                                     /**< Makes u := u-v */
void multiply_vec(vector u,double alpha);                                       /**< Makes u := u*alpha */
double product_vec(const vector u,const vector v);                              /**< Returns dot product of \f$ u \cdot v = \sum_i u_i v_i \f$ */
void dot_multiply_vec(vector which,const vector b);                             /**< Makes MATLAB operation componentwise multiple:  which(i)=which(i)*b(i) for all i */
void lincomb_vec(vector u,const vector v,const double alpha,const double beta); /**< Linear combination u:= alpha * u + beta * v */
void copy_vec(vector u,const vector v);                                         /**< Copy vector v to u , i.e. u:=v */
int vec_isnan(vector u);
/*@}*/

#endif

---