vec.h

#pragma once
 
#ifdef __cplusplus
extern "C" {
#endif
 

void vec_copy(double *y, double const * x, int n)
{
    for(int i=0;i<n;i++) y[i]=x[i];
}

void vec_fill(double *y, double x, int n)
{
    for(int i=0;i<n;i++) y[i]=x;
}
 
 
#ifdef __cplusplus
}
#endif