#include "scicos_block4.h"
#include "../machine.h"
void cumsum_r(scicos_block *block,int flag)
{
double *u;
double *y;
int nu,mu,i,j,ij;
mu=GetInPortRows(block,1);
nu=GetInPortCols(block,1);
u=GetRealInPortPtrs(block,1);
y=GetRealOutPortPtrs(block,1);
for(j=0;j<nu;j++)
{i=0;
ij=i+j*mu;
y[ij]=u[ij];
for (i=1;i<mu;i++)
{ij=i+j*mu;
y[ij]=u[ij]+y[ij-1];}}
}