GRASS Programmer's Manual  6.5.svn(2014)-r66266
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
mmul.c
Go to the documentation of this file.
1 /* mmul.c CCMATH mathematics library source code.
2  *
3  * Copyright (C) 2000 Daniel A. Atkinson All rights reserved.
4  * This code may be redistributed under the terms of the GNU library
5  * public license (LGPL). ( See the lgpl.license file for details.)
6  * ------------------------------------------------------------------------
7  */
8 #include "ccmath.h"
9 void mmul(double *c, double *a, double *b, int n)
10 {
11  double *p, *q, s;
12 
13  int i, j, k;
14 
15  trnm(b, n);
16  for (i = 0; i < n; ++i, a += n) {
17  for (j = 0, q = b; j < n; ++j) {
18  for (k = 0, p = a, s = 0.; k < n; ++k)
19  s += *p++ * *q++;
20  *c++ = s;
21  }
22  }
23  trnm(b, n);
24 }
float b
Definition: named_colr.c:8
void trnm(double *a, int n)
Definition: trnm.c:8
tuple q
Definition: forms.py:2019
void mmul(double *mat, double *a, double *b, int n)
Definition: mmul.c:9
int n
Definition: dataquad.c:291