File: m_mult.c

package info (click to toggle)
grass 6.0.2-6
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 40,044 kB
  • ctags: 31,303
  • sloc: ansic: 321,125; tcl: 25,676; sh: 11,176; cpp: 10,098; makefile: 5,025; fortran: 1,846; yacc: 493; lex: 462; perl: 133; sed: 1
file content (26 lines) | stat: -rw-r--r-- 408 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
 
/*  @(#)m_mult.c	2.1  6/26/87  */
#include <stdio.h>
#include "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;
}