#include <octave/oct.h>
#include "config.h"
#include "oct-obj.h"
#include "oct-map.h"
#include "defun-dld.h"
#include "toplev.h"

DEFUN_DLD (newkron, args, ,"Faster Kronecker Product")
{

  octave_value_list retval;

  Matrix A = args(0).matrix_value();
  Matrix B = args(1).matrix_value();
  
  int ra = A.rows();
  int ca = A.columns();
  int rb = B.rows();
  int cb = B.columns();
  int r = ra*rb;
  int c = ca*cb;
  Matrix tmpmat(rb,cb,0);
  Matrix C(r,c,0);
  for (int i=0; i<ra; i++)
    {
      for (int j=0; j<ca; j++)
	{
	  tmpmat = A(i,j)*B;
	  C = C.insert(tmpmat,i*rb,j*cb);
	}
    }
  retval(0) = octave_value(C);
return retval;

}