diff options
Diffstat (limited to '2.3-1/src/c')
4 files changed, 100 insertions, 0 deletions
diff --git a/2.3-1/src/c/elementaryFunctions/includes/logspace.h b/2.3-1/src/c/elementaryFunctions/includes/logspace.h new file mode 100644 index 00000000..1f73402b --- /dev/null +++ b/2.3-1/src/c/elementaryFunctions/includes/logspace.h @@ -0,0 +1,24 @@ +/* +SCilab2C FOSSEE IIT Bombay + */ + +#ifndef __LOGSPACE_H__ +#define __LOGSPACE_H__ + +#include "dynlib_elementaryfunctions.h" +#ifdef __cplusplus +extern "C" { +#endif + + +EXTERN_ELEMFUNCT void dlogspaces(double low_limit,double up_limit,double range_num, double* out); + +EXTERN_ELEMFUNCT void dlogspacea(double *low_limit,int row,double *up_limit,double range_num, double* out); + +#ifdef __cplusplus +} /* extern "C" */ +#endif + + +#endif /* !__LOGSPACE_H__ */ + diff --git a/2.3-1/src/c/elementaryFunctions/interfaces/int_logspace.h b/2.3-1/src/c/elementaryFunctions/interfaces/int_logspace.h new file mode 100644 index 00000000..34e325ac --- /dev/null +++ b/2.3-1/src/c/elementaryFunctions/interfaces/int_logspace.h @@ -0,0 +1,10 @@ +/* Scilab2C FOSSEE IITB */ + +#ifndef __INT_LOGSPACE_H__ +#define __INT_LOGSPACE_H__ + +#define d0d0d0logspaced2(in1,in2,in3,out) dlogspaces(in1,in2,in3,out) + +#define d2d2d0logspaced2(in1,size1,in2,size2,in3,out) dlogspacea(in1,size1[0],in2,in3,out) + +#endif diff --git a/2.3-1/src/c/elementaryFunctions/logspace/dlogspacea.c b/2.3-1/src/c/elementaryFunctions/logspace/dlogspacea.c new file mode 100644 index 00000000..f8b283c6 --- /dev/null +++ b/2.3-1/src/c/elementaryFunctions/logspace/dlogspacea.c @@ -0,0 +1,40 @@ +/* + Scilab2C FOSSEE IIT Bombay + */ + +#include "logspace.h" +#include<math.h> +void dlogspacea(double *low_limit,int _row,double *up_limit,double range_num,double *out) +{ + int i,j,k; + double temp; + double step_iterate[_row]; // for each row the spacing between two values is different. + for(i=0;i<_row;i++) + { + step_iterate[i] = pow(10,((up_limit[i]-low_limit[i])/(range_num-1))); + + + } + for(j=0;j < _row;j++) + { + out[j] = pow(10,low_limit[j]); // For every row first element is equal to 10 raise to the first value of low_limit array + temp = out[j]; + for(k=1;k < (double)range_num;k++ ) + { + out[(_row*k)+j] = temp*step_iterate[j]; /* Output matrix positions for 3 X 5 matrix are [0 3 6 9 12;1 4 7 10 13;2 5 8 11 14] so (_row*k)+j) used*/ + temp = out[(_row*k)+j]; + if(k == (double)range_num-1 ) + { + out[(_row*k)+j] = pow(10,((double)up_limit[j])); // For every row Last value of output is equal to 10 raise to first value of up_limit array + } + + } + + + } + + + + +} + diff --git a/2.3-1/src/c/elementaryFunctions/logspace/dlogspaces.c b/2.3-1/src/c/elementaryFunctions/logspace/dlogspaces.c new file mode 100644 index 00000000..7970321b --- /dev/null +++ b/2.3-1/src/c/elementaryFunctions/logspace/dlogspaces.c @@ -0,0 +1,26 @@ +/* + Scilab2C FOSSEE IIT Bombay + */ + +#include "logspace.h" +#include<math.h> +void dlogspaces(double low_limit,double up_limit,double range_num,double *out) +{ + int j; + double temp = pow(10,low_limit); + double step_iterate = pow(10,((up_limit-low_limit)/(range_num-1))); + out[0] = pow(10,low_limit); //First value of output is equal to low_limit value + for(j=1; j<(double)range_num; j++) + { + out[j] = temp*step_iterate; + temp = out[j]; + if(j == (double)range_num-1 ) + { + out[j] = pow(10,((double)up_limit)); // Last value of output is equal to up_limit value + } + } + + + +} + |