summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/elementaryFunctions/logspace/dlogspacea.c
blob: f8b283c632cb67c62594c0199b30c07ecec608bc (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
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
         }  

      }
     
 
   }
   
   
     

}