// demo
// print matrix
#include <stdio.h>
#include "lob.h"
typedef double R;
void pm(char * s, int m, int n, R w[m][n])
    {int i, j; for(i=0; i<m; ++i) {
    for(j=0; j<n; ++j) printf("%9.6f ", w[i][j]);
    printf(" %s\n", s);}}

int main(){R p[5][7]; {int i, j;
  for(i=0; i<5; ++i) for(j=0; j<7; ++j) p[i][j] = rnd();}
  {pm("b", 5, 7, p); grom(5, 7, p); pm("a", 5, 7, p);}
  {R s[7][5]; {int i, j; for(i=0; i<5; ++i) for(j=0; j<7; ++j)
     s[j][i] = p[i][j];}
  {R w[5][5]; mm(5, 7, 5, p, s, w); pm("t", 5, 5, w);}}
return 0;}