/*
 *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 *  Copyright (C) 2008 - INRIA - Arnaud TORSET
 *
 *  This file must be used under the terms of the CeCILL.
 *  This source file is licensed as described in the file COPYING, which
 *  you should have received as part of this distribution.  The terms
 *  are also available at
 *  http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 *
 */


#include "matrixMagnitude.h"

double zmagna(doubleComplex* in, int rows, int cols){
	int i=0,j=0;
	double  out=0, colSum=0;
		
		
	/* Other method : 
		zrowsuma(in,rows,cols,temp);
		out=max(temp,cols);
	but we have to malloc a array */
	
	
	
	if ((rows==1)||(cols==1)){
		for(i=0;i<cols*rows;i++){
			out += zmagns(in[i]);
		}
	}
	else{		
		for(i=0;i<cols;i++){
			colSum = 0;
			for(j=0;j<rows;j++){
				colSum += zmagns(in[i*rows+j]);
			}
			if (colSum>out) out=colSum;
		}
	}
	return out;
}