#include "EXTERN.h"
#include "perl.h"
#include "XSUB.h"


MODULE = test_c  PACKAGE = test_c

double
do_matrix (double &in_matrix,long max_i,long max_j,OUTLIST out_matrix[i][j])
PROTOTYPE: \@$$
CODE:
  long i,j,a;
  double value,out_matrix[i][j];

  for (i=0;i<=max_i;i++) {
    for (j=(i-1);j>=0;j--) {
      value=0;
      for (a=0;a<=max_j;a++) {
        value=value+in_matrix[i][a]*in_matrix[j][a];
      }
      value=value/max_j;
      out_matrix[i][j]=value;
    }
  }