File: m_mult.c

package info (click to toggle)
grass 6.4.4-1
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 104,028 kB
  • ctags: 40,409
  • sloc: ansic: 419,980; python: 63,559; tcl: 46,692; cpp: 29,791; sh: 18,564; makefile: 7,000; xml: 3,505; yacc: 561; perl: 559; lex: 480; sed: 70; objc: 7
file content (24 lines) | stat: -rw-r--r-- 394 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
 
/*  @(#)m_mult.c        2.1  6/26/87  */
#include <stdio.h>
#include <grass/libtrans.h>

#define		N	3

/*
 * m_mult: matrix multiplication (return c = a * b)
 *  3x3 matric by 3x1 matric
 */

int m_mult(double a[N][N], double b[N], double c[N])
{
    register int i, j;

    for (i = 0; i < N; i++) {
	c[i] = 0.0;

	for (j = 0; j < N; j++)
	    c[i] += (a[i][j] * b[j]);
    }

    return 1;
}