Scilab [trunk]: scilab/modules/elementaries_functions/src/c/vfrexp.c Source File

00001 #include <math.h>
00002 #include "machine.h"
00003 
00004 void C2F(vfrexp)(int *n, double *x, int *ix,double *y,int *iy,double *z, int *iz);
00005 
00006 void C2F(vfrexp)(int *n, double *x, int *ix,double *y,int *iy,double *z, int *iz)
00007 {
00008   int ix1,iy1,iz1,i,j;
00009   ix1=0;
00010   iy1=0;
00011   iz1=0;
00012 
00013   if (*ix<0) ix1=-(*n-1)*(*ix);
00014   if (*iy<0) iy1=-(*n-1)*(*iy);
00015   if (*iz<0) iz1=-(*n-1)*(*iz);
00016 
00017   for (i=0;i<*n;i++) {
00018     y[iy1]=frexp(x[ix1],&j);
00019     z[iz1]=j;
00020     iy1 += *iy;
00021     ix1 += *ix;
00022     iz1 += *iz;
00023   }
00024 }