/*--------------------------------------------------------------------------*/
SCICOS_BLOCKS_IMPEXP void m_frequ(scicos_block *block, int flag)
{
double *mat = NULL;
double *Dt = NULL;
double *off = NULL;
SCSINT32_COP *icount = NULL;
double t = 0.0;
time_counter_t** work = (time_counter_t**) block->work;
time_counter_t *counter = NULL;
int m = 0;
mat = GetRealOparPtrs(block, 1);
Dt = GetRealOparPtrs(block, 2);
off = GetRealOparPtrs(block, 3);
icount = Getint32OparPtrs(block, 4);
m = GetOparSize(block, 1, 1);
switch (flag)
{
case 4 : /* the workspace is used to store discrete counter value */
{
if ((*work = (time_counter_t*) scicos_malloc(sizeof(time_counter_t) * 2)) == NULL)
{
set_block_error(-16);
return;
}
counter = *work;
*counter = *icount;
(*(counter + 1)) = 0;
break;
}
/* event date computation */
case 3 :
{
counter = *work;
t = get_scicos_time();
*counter += (int)mat[*(counter + 1)]; /*increase counter*/
block->evout[(int)mat[*(counter + 1) + m] - 1] = *off + ((double) * counter / (*Dt)) - t;
(*(counter + 1))++;
*(counter + 1) = *(counter + 1) % m;
break;
}
/* finish */
case 5 :
{
scicos_free(*work); /*free the workspace*/
break;
}
default :
break;
}
}
/*--------------------------------------------------------------------------*/
SCICOS_BLOCKS_IMPEXP void bit_set_32(scicos_block *block, int flag)
{
int n, m, i;
SCSINT32_COP *opar;
SCSINT32_COP *u, *y;
opar = Getint32OparPtrs(block, 1);
u = Getint32InPortPtrs(block, 1);
y = Getint32OutPortPtrs(block, 1);
n = GetInPortCols(block, 1);
m = GetInPortRows(block, 1);
for (i = 0; i < m * n; i++)
{
*(y + i) = ((*(u + i)) | (*opar));
}
}
/*--------------------------------------------------------------------------*/
SCICOS_BLOCKS_IMPEXP void bit_clear_32(scicos_block *block, int flag)
{
int n = 0, m = 0, i = 0;
SCSINT32_COP *opar = NULL;
SCSINT32_COP *u = NULL, *y = NULL;
opar = Getint32OparPtrs(block, 1);
u = Getint32InPortPtrs(block, 1);
y = Getint32OutPortPtrs(block, 1);
m = GetInPortRows(block, 1);
n = GetInPortCols(block, 1);
for (i = 0; i < m * n; i++)
{
*(y + i) = ((*(u + i)) & (*opar));
}
}
/*--------------------------------------------------------------------------*/
SCICOS_BLOCKS_IMPEXP void gainblk_i32s(scicos_block *block,int flag)
{
if ((flag==1)|(flag==6)){
int i = 0,j = 0,l = 0,ji = 0,jl = 0,il = 0;
long *u = NULL,*y = NULL;
int mu = 0,ny = 0,my = 0,mo = 0,no = 0;
long *opar = NULL;
double k = 0.,D = 0.,C = 0.;
mo=GetOparSize(block,1,1);
no=GetOparSize(block,1,2);
mu=GetInPortRows(block,1);
my=GetOutPortRows(block,1);
ny=GetOutPortCols(block,1);
u=Getint32InPortPtrs(block,1);
y=Getint32OutPortPtrs(block,1);
opar=Getint32OparPtrs(block,1);
k=pow(2,32)/2;
if (mo*no==1){
for (i=0;i<ny*mu;++i){
D=(double)(opar[0])*(double)(u[i]);
if (D>=k) D=k-1;
else if (D<-k) D=-k;
y[i]=(long)D;
}
}else{
for (l=0;l<ny;l++)
{for (j=0;j<my;j++)
{D=0;
jl=j+l*my;
for (i=0;i<mu;i++)
{ji=j+i*my;
il=i+l*mu;
C=(double)(opar[ji])*(double)(u[il]);
D=D + C;}
if (D>=k) D=k-1;
else if (D<-k) D=-k;
y[jl]=(long)D;
}
}
}
}
}
