GRASS GIS 7 Programmer's Manual
7.4.4(2019)-exported
cmmult.c
Go to the documentation of this file.
1
/* cmmult.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 <stdlib.h>
9
#include "
ccmath.h
"
10
void
cmmult
(
Cpx
* cm,
Cpx
* a,
Cpx
*
b
,
int
n,
int
m,
int
l
)
11
{
12
Cpx
z, *q0, *p, *q;
13
14
int
i, j, k;
15
16
q0 = (
Cpx
*) calloc(m,
sizeof
(
Cpx
));
17
for
(i = 0; i <
l
; ++i, ++cm) {
18
for
(k = 0, p = b + i; k < m; p +=
l
)
19
q0[k++] = *p;
20
for
(j = 0, p = a, q = cm; j < n; ++j, q +=
l
) {
21
for
(k = 0, z.
re
= z.
im
= 0.; k < m; ++k, ++p) {
22
z.
re
+= p->
re
* q0[k].
re
- p->
im
* q0[k].
im
;
23
z.
im
+= p->
im
* q0[k].
re
+ p->
re
* q0[k].
im
;
24
}
25
*q = z;
26
}
27
}
28
free(q0);
29
}
l
double l
Definition:
driver/set_window.c:5
ccmath.h
b
double b
Definition:
driver/set_window.c:5
complex::re
double re
Definition:
ccmath.h:38
complex::im
double im
Definition:
ccmath.h:38
cmmult
void cmmult(Cpx *cm, Cpx *a, Cpx *b, int n, int m, int l)
Definition:
cmmult.c:10
complex
Definition:
ccmath.h:38
external
ccmath
cmmult.c
Generated on Tue Jan 15 2019 17:11:33 for GRASS GIS 7 Programmer's Manual by
1.8.13