Computational routine
eng

##
cumsum_r

### File content

    #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];}}
    }