Mercurial > octave-nkf

comparison liboctave/fCMatrix.h @ 10312:cbc402e64d83

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
untabify liboctave header files
author John W. Eaton <jwe@octave.org>
date Thu, 11 Feb 2010 12:14:48 -0500
parents 9e0ec19df4bc
children 12884915a8e4
comparison
equal deleted inserted replaced
10311:a217e1d74353 10312:cbc402e64d83
154 154
155 FloatComplexColumnVector column (octave_idx_type i) const; 155 FloatComplexColumnVector column (octave_idx_type i) const;
156 156
157 private: 157 private:
158 FloatComplexMatrix tinverse (MatrixType &mattype, octave_idx_type& info, 158 FloatComplexMatrix tinverse (MatrixType &mattype, octave_idx_type& info,
159 float& rcon, int force, int calc_cond) const; 159 float& rcon, int force, int calc_cond) const;
160 160
161 FloatComplexMatrix finverse (MatrixType &mattype, octave_idx_type& info, 161 FloatComplexMatrix finverse (MatrixType &mattype, octave_idx_type& info,
162 float& rcon, int force, int calc_cond) const; 162 float& rcon, int force, int calc_cond) const;
163 163
164 public: 164 public:
165 FloatComplexMatrix inverse (void) const; 165 FloatComplexMatrix inverse (void) const;
166 FloatComplexMatrix inverse (octave_idx_type& info) const; 166 FloatComplexMatrix inverse (octave_idx_type& info) const;
167 FloatComplexMatrix inverse (octave_idx_type& info, float& rcon, int force = 0, 167 FloatComplexMatrix inverse (octave_idx_type& info, float& rcon, int force = 0,
168 int calc_cond = 1) const; 168 int calc_cond = 1) const;
169 169
170 FloatComplexMatrix inverse (MatrixType &mattype) const; 170 FloatComplexMatrix inverse (MatrixType &mattype) const;
171 FloatComplexMatrix inverse (MatrixType &mattype, octave_idx_type& info) const; 171 FloatComplexMatrix inverse (MatrixType &mattype, octave_idx_type& info) const;
172 FloatComplexMatrix inverse (MatrixType &mattype, octave_idx_type& info, 172 FloatComplexMatrix inverse (MatrixType &mattype, octave_idx_type& info,
173 float& rcon, int force = 0, 173 float& rcon, int force = 0,
174 int calc_cond = 1) const; 174 int calc_cond = 1) const;
175 175
176 FloatComplexMatrix pseudo_inverse (float tol = 0.0) const; 176 FloatComplexMatrix pseudo_inverse (float tol = 0.0) const;
177 177
178 FloatComplexMatrix fourier (void) const; 178 FloatComplexMatrix fourier (void) const;
179 FloatComplexMatrix ifourier (void) const; 179 FloatComplexMatrix ifourier (void) const;
191 float rcond (MatrixType &mattype) const; 191 float rcond (MatrixType &mattype) const;
192 192
193 private: 193 private:
194 // Upper triangular matrix solvers 194 // Upper triangular matrix solvers
195 FloatComplexMatrix utsolve (MatrixType &typ, const FloatComplexMatrix& b, 195 FloatComplexMatrix utsolve (MatrixType &typ, const FloatComplexMatrix& b,
196 octave_idx_type& info, float& rcon, 196 octave_idx_type& info, float& rcon,
197 solve_singularity_handler sing_handler, 197 solve_singularity_handler sing_handler,
198 bool calc_cond = false, blas_trans_type transt = blas_no_trans) const; 198 bool calc_cond = false, blas_trans_type transt = blas_no_trans) const;
199 199
200 // Lower triangular matrix solvers 200 // Lower triangular matrix solvers
201 FloatComplexMatrix ltsolve (MatrixType &typ, const FloatComplexMatrix& b, 201 FloatComplexMatrix ltsolve (MatrixType &typ, const FloatComplexMatrix& b,
202 octave_idx_type& info, float& rcon, 202 octave_idx_type& info, float& rcon,
203 solve_singularity_handler sing_handler, 203 solve_singularity_handler sing_handler,
204 bool calc_cond = false, blas_trans_type transt = blas_no_trans) const; 204 bool calc_cond = false, blas_trans_type transt = blas_no_trans) const;
205 205
206 // Full matrix solvers (umfpack/cholesky) 206 // Full matrix solvers (umfpack/cholesky)
207 FloatComplexMatrix fsolve (MatrixType &typ, const FloatComplexMatrix& b, 207 FloatComplexMatrix fsolve (MatrixType &typ, const FloatComplexMatrix& b,
208 octave_idx_type& info, float& rcon, 208 octave_idx_type& info, float& rcon,
209 solve_singularity_handler sing_handler, 209 solve_singularity_handler sing_handler,
210 bool calc_cond = false) const; 210 bool calc_cond = false) const;
211 211
212 public: 212 public:
213 // Generic interface to solver with no probing of type 213 // Generic interface to solver with no probing of type
214 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b) const; 214 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b) const;
215 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b, 215 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b,
216 octave_idx_type& info) const; 216 octave_idx_type& info) const;
217 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b, 217 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b,
218 octave_idx_type& info, float& rcon) const; 218 octave_idx_type& info, float& rcon) const;
219 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b, octave_idx_type& info, 219 FloatComplexMatrix solve (MatrixType &typ, const FloatMatrix& b, octave_idx_type& info,
220 float& rcon, solve_singularity_handler sing_handler, 220 float& rcon, solve_singularity_handler sing_handler,
221 bool singular_fallback = true, 221 bool singular_fallback = true,
222 blas_trans_type transt = blas_no_trans) const; 222 blas_trans_type transt = blas_no_trans) const;
223 223
224 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b) const; 224 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b) const;
225 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b, 225 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b,
226 octave_idx_type& info) const; 226 octave_idx_type& info) const;
227 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b, 227 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b,
228 octave_idx_type& info, float& rcon) const; 228 octave_idx_type& info, float& rcon) const;
229 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b, 229 FloatComplexMatrix solve (MatrixType &typ, const FloatComplexMatrix& b,
230 octave_idx_type& info, float& rcon, 230 octave_idx_type& info, float& rcon,
231 solve_singularity_handler sing_handler, 231 solve_singularity_handler sing_handler,
232 bool singular_fallback = true, 232 bool singular_fallback = true,
233 blas_trans_type transt = blas_no_trans) const; 233 blas_trans_type transt = blas_no_trans) const;
234 234
235 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b) const; 235 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b) const;
236 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b, 236 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b,
237 octave_idx_type& info) const; 237 octave_idx_type& info) const;
238 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b, 238 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b,
239 octave_idx_type& info, float& rcon) const; 239 octave_idx_type& info, float& rcon) const;
240 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b, 240 FloatComplexColumnVector solve (MatrixType &typ, const FloatColumnVector& b,
241 octave_idx_type& info, float& rcon, 241 octave_idx_type& info, float& rcon,
242 solve_singularity_handler sing_handler, 242 solve_singularity_handler sing_handler,
243 blas_trans_type transt = blas_no_trans) const; 243 blas_trans_type transt = blas_no_trans) const;
244 244
245 FloatComplexColumnVector solve (MatrixType &typ, 245 FloatComplexColumnVector solve (MatrixType &typ,
246 const FloatComplexColumnVector& b) const; 246 const FloatComplexColumnVector& b) const;
247 FloatComplexColumnVector solve (MatrixType &typ, const FloatComplexColumnVector& b, 247 FloatComplexColumnVector solve (MatrixType &typ, const FloatComplexColumnVector& b,
248 octave_idx_type& info) const; 248 octave_idx_type& info) const;
249 FloatComplexColumnVector solve (MatrixType &typ, const FloatComplexColumnVector& b, 249 FloatComplexColumnVector solve (MatrixType &typ, const FloatComplexColumnVector& b,
250 octave_idx_type& info, float& rcon) const; 250 octave_idx_type& info, float& rcon) const;
251 FloatComplexColumnVector solve (MatrixType &typ, const FloatComplexColumnVector& b, 251 FloatComplexColumnVector solve (MatrixType &typ, const FloatComplexColumnVector& b,
252 octave_idx_type& info, float& rcon, 252 octave_idx_type& info, float& rcon,
253 solve_singularity_handler sing_handler, 253 solve_singularity_handler sing_handler,
254 blas_trans_type transt = blas_no_trans) const; 254 blas_trans_type transt = blas_no_trans) const;
255 255
256 // Generic interface to solver with probing of type 256 // Generic interface to solver with probing of type
257 FloatComplexMatrix solve (const FloatMatrix& b) const; 257 FloatComplexMatrix solve (const FloatMatrix& b) const;
258 FloatComplexMatrix solve (const FloatMatrix& b, octave_idx_type& info) const; 258 FloatComplexMatrix solve (const FloatMatrix& b, octave_idx_type& info) const;
259 FloatComplexMatrix solve (const FloatMatrix& b, octave_idx_type& info, float& rcon) const; 259 FloatComplexMatrix solve (const FloatMatrix& b, octave_idx_type& info, float& rcon) const;
260 FloatComplexMatrix solve (const FloatMatrix& b, octave_idx_type& info, float& rcon, 260 FloatComplexMatrix solve (const FloatMatrix& b, octave_idx_type& info, float& rcon,
261 solve_singularity_handler sing_handler, 261 solve_singularity_handler sing_handler,
262 blas_trans_type transt = blas_no_trans) const; 262 blas_trans_type transt = blas_no_trans) const;
263 263
264 FloatComplexMatrix solve (const FloatComplexMatrix& b) const; 264 FloatComplexMatrix solve (const FloatComplexMatrix& b) const;
265 FloatComplexMatrix solve (const FloatComplexMatrix& b, octave_idx_type& info) const; 265 FloatComplexMatrix solve (const FloatComplexMatrix& b, octave_idx_type& info) const;
266 FloatComplexMatrix solve (const FloatComplexMatrix& b, octave_idx_type& info, float& rcon) const; 266 FloatComplexMatrix solve (const FloatComplexMatrix& b, octave_idx_type& info, float& rcon) const;
267 FloatComplexMatrix solve (const FloatComplexMatrix& b, octave_idx_type& info, float& rcon, 267 FloatComplexMatrix solve (const FloatComplexMatrix& b, octave_idx_type& info, float& rcon,
268 solve_singularity_handler sing_handler, 268 solve_singularity_handler sing_handler,
269 blas_trans_type transt = blas_no_trans) const; 269 blas_trans_type transt = blas_no_trans) const;
270 270
271 FloatComplexColumnVector solve (const FloatColumnVector& b) const; 271 FloatComplexColumnVector solve (const FloatColumnVector& b) const;
272 FloatComplexColumnVector solve (const FloatColumnVector& b, octave_idx_type& info) const; 272 FloatComplexColumnVector solve (const FloatColumnVector& b, octave_idx_type& info) const;
273 FloatComplexColumnVector solve (const FloatColumnVector& b, octave_idx_type& info, 273 FloatComplexColumnVector solve (const FloatColumnVector& b, octave_idx_type& info,
274 float& rcon) const; 274 float& rcon) const;
275 FloatComplexColumnVector solve (const FloatColumnVector& b, octave_idx_type& info, float& rcon, 275 FloatComplexColumnVector solve (const FloatColumnVector& b, octave_idx_type& info, float& rcon,
276 solve_singularity_handler sing_handler, 276 solve_singularity_handler sing_handler,
277 blas_trans_type transt = blas_no_trans) const; 277 blas_trans_type transt = blas_no_trans) const;
278 278
279 FloatComplexColumnVector solve (const FloatComplexColumnVector& b) const; 279 FloatComplexColumnVector solve (const FloatComplexColumnVector& b) const;
280 FloatComplexColumnVector solve (const FloatComplexColumnVector& b, octave_idx_type& info) const; 280 FloatComplexColumnVector solve (const FloatComplexColumnVector& b, octave_idx_type& info) const;
281 FloatComplexColumnVector solve (const FloatComplexColumnVector& b, octave_idx_type& info, 281 FloatComplexColumnVector solve (const FloatComplexColumnVector& b, octave_idx_type& info,
282 float& rcon) const; 282 float& rcon) const;
283 FloatComplexColumnVector solve (const FloatComplexColumnVector& b, octave_idx_type& info, 283 FloatComplexColumnVector solve (const FloatComplexColumnVector& b, octave_idx_type& info,
284 float& rcon, 284 float& rcon,
285 solve_singularity_handler sing_handler, 285 solve_singularity_handler sing_handler,
286 blas_trans_type transt = blas_no_trans) const; 286 blas_trans_type transt = blas_no_trans) const;
287 287
288 FloatComplexMatrix lssolve (const FloatMatrix& b) const; 288 FloatComplexMatrix lssolve (const FloatMatrix& b) const;
289 FloatComplexMatrix lssolve (const FloatMatrix& b, octave_idx_type& info) const; 289 FloatComplexMatrix lssolve (const FloatMatrix& b, octave_idx_type& info) const;
290 FloatComplexMatrix lssolve (const FloatMatrix& b, octave_idx_type& info, 290 FloatComplexMatrix lssolve (const FloatMatrix& b, octave_idx_type& info,
291 octave_idx_type& rank) const; 291 octave_idx_type& rank) const;
292 FloatComplexMatrix lssolve (const FloatMatrix& b, octave_idx_type& info, 292 FloatComplexMatrix lssolve (const FloatMatrix& b, octave_idx_type& info,
293 octave_idx_type& rank, float& rcon) const; 293 octave_idx_type& rank, float& rcon) const;
294 294
295 FloatComplexMatrix lssolve (const FloatComplexMatrix& b) const; 295 FloatComplexMatrix lssolve (const FloatComplexMatrix& b) const;
296 FloatComplexMatrix lssolve (const FloatComplexMatrix& b, octave_idx_type& info) const; 296 FloatComplexMatrix lssolve (const FloatComplexMatrix& b, octave_idx_type& info) const;
297 FloatComplexMatrix lssolve (const FloatComplexMatrix& b, octave_idx_type& info, 297 FloatComplexMatrix lssolve (const FloatComplexMatrix& b, octave_idx_type& info,
298 octave_idx_type& rank) const; 298 octave_idx_type& rank) const;
299 FloatComplexMatrix lssolve (const FloatComplexMatrix& b, octave_idx_type& info, 299 FloatComplexMatrix lssolve (const FloatComplexMatrix& b, octave_idx_type& info,
300 octave_idx_type& rank, float& rcon) const; 300 octave_idx_type& rank, float& rcon) const;
301 301
302 FloatComplexColumnVector lssolve (const FloatColumnVector& b) const; 302 FloatComplexColumnVector lssolve (const FloatColumnVector& b) const;
303 FloatComplexColumnVector lssolve (const FloatColumnVector& b, 303 FloatComplexColumnVector lssolve (const FloatColumnVector& b,
304 octave_idx_type& info) const; 304 octave_idx_type& info) const;
305 FloatComplexColumnVector lssolve (const FloatColumnVector& b, octave_idx_type& info, 305 FloatComplexColumnVector lssolve (const FloatColumnVector& b, octave_idx_type& info,
306 octave_idx_type& rank) const; 306 octave_idx_type& rank) const;
307 FloatComplexColumnVector lssolve (const FloatColumnVector& b, octave_idx_type& info, 307 FloatComplexColumnVector lssolve (const FloatColumnVector& b, octave_idx_type& info,
308 octave_idx_type& rank, float& rcon) const; 308 octave_idx_type& rank, float& rcon) const;
309 309
310 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b) const; 310 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b) const;
311 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b, 311 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b,
312 octave_idx_type& info) const; 312 octave_idx_type& info) const;
313 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b, 313 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b,
314 octave_idx_type& info, 314 octave_idx_type& info,
315 octave_idx_type& rank) const; 315 octave_idx_type& rank) const;
316 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b, 316 FloatComplexColumnVector lssolve (const FloatComplexColumnVector& b,
317 octave_idx_type& info, 317 octave_idx_type& info,
318 octave_idx_type& rank, float& rcon) const; 318 octave_idx_type& rank, float& rcon) const;
319 319
320 // matrix by diagonal matrix -> matrix operations 320 // matrix by diagonal matrix -> matrix operations
321 321
322 FloatComplexMatrix& operator += (const FloatDiagMatrix& a); 322 FloatComplexMatrix& operator += (const FloatDiagMatrix& a);
323 FloatComplexMatrix& operator -= (const FloatDiagMatrix& a); 323 FloatComplexMatrix& operator -= (const FloatDiagMatrix& a);