#include <stdio.h>
#include <math.h>
#include "tictoc.h"

void pmatrix(int n,double A[n][n]){
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			printf("%6g ",A[i][j]);
		}
		printf("\n");
	}
}
double fnorm(int n,double A[n][n]){
	double r=0;
	for(int i=0;i<n;i++){
		for(int j=0;j<n;j++){
			r+=A[i][j]*A[i][j];
		}
	}
	return sqrt(r);
}

#define N 10000
double A[N][N];

int main(){
	for(int i=0;i<N;i++){
		for(int j=0;j<N;j++){
			A[i][j]=1.0/(1+i+j);
		}
	}
	// pmatrix(N,A);
	tic();
	double r=fnorm(N,A);
	double t=toc();
	printf("The norm is %.14e\n",r);
	printf("Took %g seconds\n",t);
	return 0;
}