25 #define P_ASSIGN2(x,z,p,ax,az,q)	x [p] = 1  argument
33 #define R_ASSEMBLE(x,z,p,ax,az,q)		x [p] += ax [q]  argument
34 #define R_ASSIGN(x,z,p,ax,az,q)			x [p]  = ax [q]  argument
35 #define R_ASSIGN_CONJ(x,z,p,ax,az,q)		x [p]  = ax [q]  argument
36 #define R_ASSIGN_REAL(x,p,ax,q)			x [p]  = ax [q]  argument
38 #define R_IS_NONZERO(ax,az,q)			IS_NONZERO (ax [q])  argument
39 #define R_IS_ZERO(ax,az,q)			IS_ZERO (ax [q])  argument
40 #define R_IS_ONE(ax,az,q)			(ax [q] == 1)  argument
41 #define R_MULT(x,z,p, ax,az,q, bx,bz,r)		x [p]  = ax [q] * bx [r]  argument
42 #define R_MULTADD(x,z,p, ax,az,q, bx,bz,r)	x [p] += ax [q] * bx [r]  argument
43 #define R_MULTSUB(x,z,p, ax,az,q, bx,bz,r)	x [p] -= ax [q] * bx [r]  argument
44 #define R_MULTADDCONJ(x,z,p, ax,az,q, bx,bz,r)	x [p] += ax [q] * bx [r]  argument
45 #define R_MULTSUBCONJ(x,z,p, ax,az,q, bx,bz,r)	x [p] -= ax [q] * bx [r]  argument
46 #define R_ADD(x,z,p, ax,az,q, bx,bz,r)		x [p]  = ax [q] + bx [r]  argument
47 #define R_ADD_REAL(x,p, ax,q, bx,r)		x [p]  = ax [q] + bx [r]  argument
50 #define R_DIV(x,z,p,ax,az,q)			x [p] /= ax [q]  argument
51 #define R_LLDOT(x,p, ax,az,q)			x [p] -= ax [q] * ax [q]  argument
54 #define R_DIV_REAL(x,z,p, ax,az,q, bx,r)	x [p] = ax [q] / bx [r]  argument
55 #define R_MULT_REAL(x,z,p, ax,az,q, bx,r)	x [p] = ax [q] * bx [r]  argument
57 #define R_LDLDOT(x,p, ax,az,q, bx,r)		x [p] -=(ax[q] * ax[q])/ bx[r]  argument
66 #define C_ASSEMBLE(x,z,p,ax,az,q) \  argument
70 #define C_ASSIGN(x,z,p,ax,az,q) \  argument
74 #define C_ASSIGN_REAL(x,p,ax,q)			x [2*(p)]  = ax [2*(q)]  argument
76 #define C_ASSIGN_CONJ(x,z,p,ax,az,q) \  argument
82 #define C_IS_NONZERO(ax,az,q) \  argument
85 #define C_IS_ZERO(ax,az,q) \  argument
88 #define C_IS_ONE(ax,az,q) \  argument
91 #define C_IMAG_IS_NONZERO(ax,az,q)  (IS_NONZERO (ax [2*(q)+1]))  argument
93 #define C_MULT(x,z,p, ax,az,q, bx,bz,r) \  argument
97 #define C_MULTADD(x,z,p, ax,az,q, bx,bz,r) \  argument
101 #define C_MULTSUB(x,z,p, ax,az,q, bx,bz,r) \  argument
106 #define C_MULTADDCONJ(x,z,p, ax,az,q, bx,bz,r) \  argument
111 #define C_MULTSUBCONJ(x,z,p, ax,az,q, bx,bz,r) \  argument
115 #define C_ADD(x,z,p, ax,az,q, bx,bz,r) \  argument
119 #define C_ADD_REAL(x,p, ax,q, bx,r) \  argument
130 #define C_DIV(x,z,p,ax,az,q) \  argument
137 #define C_LLDOT(x,p, ax,az,q) \  argument
142 #define C_DIV_REAL(x,z,p, ax,az,q, bx,r) \  argument
146 #define C_MULT_REAL(x,z,p, ax,az,q, bx,r) \  argument
151 #define C_LDLDOT(x,p, ax,az,q, bx,r) \  argument
161 #define Z_ASSEMBLE(x,z,p,ax,az,q) \  argument
165 #define Z_ASSIGN(x,z,p,ax,az,q) \  argument
169 #define Z_ASSIGN_REAL(x,p,ax,q)			x [p]  = ax [q]  argument
171 #define Z_ASSIGN_CONJ(x,z,p,ax,az,q) \  argument
177 #define Z_IS_NONZERO(ax,az,q) \  argument
180 #define Z_IS_ZERO(ax,az,q) \  argument
183 #define Z_IS_ONE(ax,az,q) \  argument
186 #define Z_IMAG_IS_NONZERO(ax,az,q)  (IS_NONZERO (az [q]))  argument
188 #define Z_MULT(x,z,p, ax,az,q, bx,bz,r) \  argument
192 #define Z_MULTADD(x,z,p, ax,az,q, bx,bz,r) \  argument
196 #define Z_MULTSUB(x,z,p, ax,az,q, bx,bz,r) \  argument
200 #define Z_MULTADDCONJ(x,z,p, ax,az,q, bx,bz,r) \  argument
204 #define Z_MULTSUBCONJ(x,z,p, ax,az,q, bx,bz,r) \  argument
208 #define Z_ADD(x,z,p, ax,az,q, bx,bz,r) \  argument
212 #define Z_ADD_REAL(x,p, ax,q, bx,r) \  argument
223 #define Z_DIV(x,z,p,ax,az,q) \  argument
228 #define Z_LLDOT(x,p, ax,az,q) \  argument
233 #define Z_DIV_REAL(x,z,p, ax,az,q, bx,r) \  argument
237 #define Z_MULT_REAL(x,z,p, ax,az,q, bx,r) \  argument
242 #define Z_LDLDOT(x,p, ax,az,q, bx,r) \  argument