// Matrix multiplication sample
//   a[n][m] x b[m][s] -> c[n][s]

void MatMlt(int N, int M, int S, double *a, double *b, double *c)
{
    int n,m,s ;
	
    for (s=0 ; s<S ; s++) {
        for (n=0 ; n<N ; n++) {
            c[S*n+s] = 0.0 ;
            for (m=0 ; m<M ; m++) {
                c[S*n+s] += a[M*n+m] * b[S*m+s] ;
            }
        }
    }
}