pdsvector.h

/*
 * pdsvector.h
 * 
 * Copyright 2011 Fernando Pujaico Rivera <fernando.pujaico.rivera@gmail.com>
 * 
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 * 
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
 * MA 02110-1301, USA.
 * 
 */

#ifndef __PDSVECTOR_H__
#define __PDSVECTOR_H__

#ifdef __cplusplus
extern "C" {
#endif 

#include <stdio.h>
#include <stdlib.h>
#include <pds/pdsraglobal.h>


#ifndef PDS_OK

        #define PDS_OK 1
#endif

#ifndef PDS_WRONG

        #define PDS_WRONG 0
#endif


typedef struct 
{
        PdsRaReal *V;
        PdsRaNatural Nel;
}PdsVector;

PdsVector *pds_vector_new(PdsRaNatural Nel);


PdsVector *pds_vector_new_array1d(const PdsRaReal *Array1D,PdsRaNatural N);


PdsVector *pds_vector_new_args(PdsRaNatural N,...);


PdsVector *pds_vector_new_vector(const PdsVector *VecSrc);


PdsVector *pds_vector_new_load_data(const char* datafile);


PdsVector *pds_vector_new_load_data_line(const char* datafile,PdsRaNatural line);


PdsVector *pds_vector_new_load_data_col(const char* datafile,PdsRaNatural col);



PdsVector *pds_vector_new_linspace(PdsRaReal xi,PdsRaReal xf,PdsRaNatural N);


PdsVector *pds_vector_new_logspace(PdsRaReal xi,PdsRaReal xf,PdsRaNatural N);


PdsVector *pds_vector_new_rectangular(PdsRaNatural N);


PdsVector *pds_vector_new_hamming(PdsRaNatural N);


PdsVector *pds_vector_new_blackman(PdsRaNatural N,PdsRaReal Alpha);



int pds_vector_init_value(PdsVector *Vector,PdsRaReal Val);


int pds_vector_init_vector(PdsVector *Vector,const PdsVector *VecSrc);




int pds_vector_add_vector(PdsVector *Vector,const PdsVector *VecSrc);


int pds_vector_add_factor_vector(PdsVector *Vector,PdsRaReal Alpha,const PdsVector *VecSrc);


int pds_vector_add_value(PdsVector *Vector,PdsRaReal Val);


int pds_vector_sub_vector(PdsVector *Vector,const PdsVector *VecSrc);


int pds_vector_mul_vector_elements(PdsVector *Vector,const PdsVector *VecSrc);


int pds_vector_mul_value(PdsVector *Vector,PdsRaReal Val);


int pds_vector_div_value(PdsVector *Vector,PdsRaReal Val);


int pds_vector_dot_vector(const PdsVector *Vector1,const PdsVector *Vector2,PdsRaReal *Val);



int pds_vector_printf(const PdsVector *Vector);


int pds_vector_fprintf(const PdsVector *Vector,FILE *fd);


int pds_vector_fscanf(PdsVector *Vector, FILE *fd);


int pds_vector_fwrite(PdsVector *Vector, FILE *fd);


int pds_vector_fread(PdsVector *Vector, FILE *fd);




int pds_vector_norm_vector(const PdsVector *Vector, PdsRaReal *m);


int pds_vector_norm2_vector(const PdsVector *Vector, PdsRaReal *m);


int pds_vector_rms_vector(const PdsVector *Vector, PdsRaReal *m);


int pds_vector_rms_vectors(const PdsVector *VectorA, const PdsVector *VectorB, PdsRaReal *m);


int pds_vector_rms_factor_vectors(const PdsVector *VectorA, const PdsVector *VectorB, PdsRaReal *m);

int pds_vector_mean_vector(const PdsVector *Vector, PdsRaReal *m);


int pds_vector_var_vector(const PdsVector *Vector, PdsRaReal *v);


int pds_vector_std_vector(const PdsVector *Vector, PdsRaReal *v);


int pds_vector_cov_vector (const PdsVector *VectorX,const PdsVector *VectorY, PdsRaReal *c);


int pds_vector_cor_vector (const PdsVector *VectorX,const PdsVector *VectorY, PdsRaReal *c);


PdsVector* pds_vector_new_conv(const PdsVector *Vector1,const PdsVector *Vector2);


int pds_vector_swap_elements(PdsVector *Vector, PdsRaNatural X, PdsRaNatural Y);


int pds_vector_asc_order(PdsVector *Vector);


int pds_vector_desc_order(PdsVector *Vector);


int pds_vector_cmp(const PdsVector *Vector1,const PdsVector *Vector2,PdsRaNatural *n);




int pds_vector_get_value(const PdsVector *Vector, PdsRaNatural x, PdsRaReal *m);


int pds_vector_get_nelements(const PdsVector *Vector,PdsRaReal *Nelements);


int pds_vector_get_min_value(const PdsVector *Vector, PdsRaReal *m);


int pds_vector_get_max_value(const PdsVector *Vector, PdsRaReal *m);



#define PDS_VSET(Vector,id,Val) pds_vector_set_value(Vector,id,Val)

int pds_vector_set_value(PdsVector *Vector, PdsRaNatural x, PdsRaReal m);


int pds_vector_set_bcd(PdsVector *Vector,PdsRaNatural n,PdsRaNatural b0,PdsRaNatural b1);


int pds_vector_set_gray(PdsVector *Vector,PdsRaNatural n,PdsRaNatural b0,PdsRaNatural b1);



int pds_vector_octave_plot(const PdsVector *X,const PdsVector *Y,const char *labelx,const char *labely,const char *octfile,const char *imgfile);


int pds_vector_octave_semilogx(const PdsVector *X,const PdsVector *Y,const char *labelx,const char *labely,const char *octfile,const char *imgfile);


int pds_vector_octave_semilogy(const PdsVector *X,const PdsVector *Y,const char *labelx,const char *labely,const char *octfile,const char *imgfile);



void pds_vector_free(PdsVector *Vector);


void pds_vector_destroy(PdsVector **Vector);


#ifdef __cplusplus
}
#endif 

#endif