Lines Matching defs:ablock
4 #define ablock(m,n) ablock[(m-1)+5*(n-1)]
13 void matvec_sub(double ablock[],double avec[],double bvec[]) {
19 // subtracts bvec=bvec - ablock*avec
24 // $ - lhs(i,1,ablock,ia,ja,ka,acell)*
26 bvec(1) = bvec(1) - ablock(1,1)*avec(1)
27 - ablock(1,2)*avec(2)
28 - ablock(1,3)*avec(3)
29 - ablock(1,4)*avec(4)
30 - ablock(1,5)*avec(5);
31 bvec(2) = bvec(2) - ablock(2,1)*avec(1)
32 - ablock(2,2)*avec(2)
33 - ablock(2,3)*avec(3)
34 - ablock(2,4)*avec(4)
35 - ablock(2,5)*avec(5);
36 bvec(3) = bvec(3) - ablock(3,1)*avec(1)
37 - ablock(3,2)*avec(2)
38 - ablock(3,3)*avec(3)
39 - ablock(3,4)*avec(4)
40 - ablock(3,5)*avec(5);
41 bvec(4) = bvec(4) - ablock(4,1)*avec(1)
42 - ablock(4,2)*avec(2)
43 - ablock(4,3)*avec(3)
44 - ablock(4,4)*avec(4)
45 - ablock(4,5)*avec(5);
46 bvec(5) = bvec(5) - ablock(5,1)*avec(1)
47 - ablock(5,2)*avec(2)
48 - ablock(5,3)*avec(3)
49 - ablock(5,4)*avec(4)
50 - ablock(5,5)*avec(5);
59 void matmul_sub(double ablock[], double bblock[], double cblock[]) {
69 cblock(1,1) = cblock(1,1) - ablock(1,1)*bblock(1,1)
70 - ablock(1,2)*bblock(2,1)
71 - ablock(1,3)*bblock(3,1)
72 - ablock(1,4)*bblock(4,1)
73 - ablock(1,5)*bblock(5,1);
74 cblock(2,1) = cblock(2,1) - ablock(2,1)*bblock(1,1)
75 - ablock(2,2)*bblock(2,1)
76 - ablock(2,3)*bblock(3,1)
77 - ablock(2,4)*bblock(4,1)
78 - ablock(2,5)*bblock(5,1);
79 cblock(3,1) = cblock(3,1) - ablock(3,1)*bblock(1,1)
80 - ablock(3,2)*bblock(2,1)
81 - ablock(3,3)*bblock(3,1)
82 - ablock(3,4)*bblock(4,1)
83 - ablock(3,5)*bblock(5,1);
84 cblock(4,1) = cblock(4,1) - ablock(4,1)*bblock(1,1)
85 - ablock(4,2)*bblock(2,1)
86 - ablock(4,3)*bblock(3,1)
87 - ablock(4,4)*bblock(4,1)
88 - ablock(4,5)*bblock(5,1);
89 cblock(5,1) = cblock(5,1) - ablock(5,1)*bblock(1,1)
90 - ablock(5,2)*bblock(2,1)
91 - ablock(5,3)*bblock(3,1)
92 - ablock(5,4)*bblock(4,1)
93 - ablock(5,5)*bblock(5,1);
94 cblock(1,2) = cblock(1,2) - ablock(1,1)*bblock(1,2)
95 - ablock(1,2)*bblock(2,2)
96 - ablock(1,3)*bblock(3,2)
97 - ablock(1,4)*bblock(4,2)
98 - ablock(1,5)*bblock(5,2);
99 cblock(2,2) = cblock(2,2) - ablock(2,1)*bblock(1,2)
100 - ablock(2,2)*bblock(2,2)
101 - ablock(2,3)*bblock(3,2)
102 - ablock(2,4)*bblock(4,2)
103 - ablock(2,5)*bblock(5,2);
104 cblock(3,2) = cblock(3,2) - ablock(3,1)*bblock(1,2)
105 - ablock(3,2)*bblock(2,2)
106 - ablock(3,3)*bblock(3,2)
107 - ablock(3,4)*bblock(4,2)
108 - ablock(3,5)*bblock(5,2);
109 cblock(4,2) = cblock(4,2) - ablock(4,1)*bblock(1,2)
110 - ablock(4,2)*bblock(2,2)
111 - ablock(4,3)*bblock(3,2)
112 - ablock(4,4)*bblock(4,2)
113 - ablock(4,5)*bblock(5,2);
114 cblock(5,2) = cblock(5,2) - ablock(5,1)*bblock(1,2)
115 - ablock(5,2)*bblock(2,2)
116 - ablock(5,3)*bblock(3,2)
117 - ablock(5,4)*bblock(4,2)
118 - ablock(5,5)*bblock(5,2);
119 cblock(1,3) = cblock(1,3) - ablock(1,1)*bblock(1,3)
120 - ablock(1,2)*bblock(2,3)
121 - ablock(1,3)*bblock(3,3)
122 - ablock(1,4)*bblock(4,3)
123 - ablock(1,5)*bblock(5,3);
124 cblock(2,3) = cblock(2,3) - ablock(2,1)*bblock(1,3)
125 - ablock(2,2)*bblock(2,3)
126 - ablock(2,3)*bblock(3,3)
127 - ablock(2,4)*bblock(4,3)
128 - ablock(2,5)*bblock(5,3);
129 cblock(3,3) = cblock(3,3) - ablock(3,1)*bblock(1,3)
130 - ablock(3,2)*bblock(2,3)
131 - ablock(3,3)*bblock(3,3)
132 - ablock(3,4)*bblock(4,3)
133 - ablock(3,5)*bblock(5,3);
134 cblock(4,3) = cblock(4,3) - ablock(4,1)*bblock(1,3)
135 - ablock(4,2)*bblock(2,3)
136 - ablock(4,3)*bblock(3,3)
137 - ablock(4,4)*bblock(4,3)
138 - ablock(4,5)*bblock(5,3);
139 cblock(5,3) = cblock(5,3) - ablock(5,1)*bblock(1,3)
140 - ablock(5,2)*bblock(2,3)
141 - ablock(5,3)*bblock(3,3)
142 - ablock(5,4)*bblock(4,3)
143 - ablock(5,5)*bblock(5,3);
144 cblock(1,4) = cblock(1,4) - ablock(1,1)*bblock(1,4)
145 - ablock(1,2)*bblock(2,4)
146 - ablock(1,3)*bblock(3,4)
147 - ablock(1,4)*bblock(4,4)
148 - ablock(1,5)*bblock(5,4);
149 cblock(2,4) = cblock(2,4) - ablock(2,1)*bblock(1,4)
150 - ablock(2,2)*bblock(2,4)
151 - ablock(2,3)*bblock(3,4)
152 - ablock(2,4)*bblock(4,4)
153 - ablock(2,5)*bblock(5,4);
154 cblock(3,4) = cblock(3,4) - ablock(3,1)*bblock(1,4)
155 - ablock(3,2)*bblock(2,4)
156 - ablock(3,3)*bblock(3,4)
157 - ablock(3,4)*bblock(4,4)
158 - ablock(3,5)*bblock(5,4);
159 cblock(4,4) = cblock(4,4) - ablock(4,1)*bblock(1,4)
160 - ablock(4,2)*bblock(2,4)
161 - ablock(4,3)*bblock(3,4)
162 - ablock(4,4)*bblock(4,4)
163 - ablock(4,5)*bblock(5,4);
164 cblock(5,4) = cblock(5,4) - ablock(5,1)*bblock(1,4)
165 - ablock(5,2)*bblock(2,4)
166 - ablock(5,3)*bblock(3,4)
167 - ablock(5,4)*bblock(4,4)
168 - ablock(5,5)*bblock(5,4);
169 cblock(1,5) = cblock(1,5) - ablock(1,1)*bblock(1,5)
170 - ablock(1,2)*bblock(2,5)
171 - ablock(1,3)*bblock(3,5)
172 - ablock(1,4)*bblock(4,5)
173 - ablock(1,5)*bblock(5,5);
174 cblock(2,5) = cblock(2,5) - ablock(2,1)*bblock(1,5)
175 - ablock(2,2)*bblock(2,5)
176 - ablock(2,3)*bblock(3,5)
177 - ablock(2,4)*bblock(4,5)
178 - ablock(2,5)*bblock(5,5);
179 cblock(3,5) = cblock(3,5) - ablock(3,1)*bblock(1,5)
180 - ablock(3,2)*bblock(2,5)
181 - ablock(3,3)*bblock(3,5)
182 - ablock(3,4)*bblock(4,5)
183 - ablock(3,5)*bblock(5,5);
184 cblock(4,5) = cblock(4,5) - ablock(4,1)*bblock(1,5)
185 - ablock(4,2)*bblock(2,5)
186 - ablock(4,3)*bblock(3,5)
187 - ablock(4,4)*bblock(4,5)
188 - ablock(4,5)*bblock(5,5);
189 cblock(5,5) = cblock(5,5) - ablock(5,1)*bblock(1,5)
190 - ablock(5,2)*bblock(2,5)
191 - ablock(5,3)*bblock(3,5)
192 - ablock(5,4)*bblock(4,5)
193 - ablock(5,5)*bblock(5,5);