-
Notifications
You must be signed in to change notification settings - Fork 6
/
Matrix_operations_multipy.c
71 lines (58 loc) · 1.43 KB
/
Matrix_operations_multipy.c
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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
/*
Copyright (C) 1995 The GeoFramework Consortium
This file is part of Ellipsis3D.
Ellipsis3D is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License, version 2,
as published by the Free Software Foundation.
Ellipsis3D is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
Author:
Louis Moresi <[email protected]>
*/
#include <math.h>
#include "element_definitions.h"
#include "global_defs.h"
/* Procedure which multiplies the m*n matrix A by the n*p matrix B
The result is the m*p matrix C*/
void mul_A_B(
standard_precision A[7][16],
standard_precision B[16][7],
standard_precision C[7][7]
)
{
int i,j,k ;
for(i=1;i<7;i++){
for(j=1;j<7;j++){
C[i][j] = 0.0 ;
for(k=1;k<16;k++){
C[i][j] += A[i][k] * B[k][j] ;
}
}
}
return ;
}
/* Procedure which multiplies the transpose of the m*n matrix A by the m*p matrix B
The result is the n*p matrix C*/
void mul_At_B(
standard_precision A[16][7],
standard_precision B[16][16],
standard_precision C[7][16]
)
{
int i,j,k ;
for(i=1;i<7;i++){
for(j=1;j<16;j++){
C[i][j] = 0.0 ;
}
}
for(i=1;i<7;i++){
for(j=1;j<16;j++){
for(k=1;k<16;k++){
C[i][j] += A[k][i] * B[k][j] ;
}
}
}
return ;
}