#include <mpi.h>
#include <mpe.h>
#include <iostream>
#include <fstream>
#include <vector>
#include <iomanip>
#include <imatrix.hpp>
#include <mvi_util.hpp>
#include "ilss_par.hpp"
Constants | |
const real | zerotest = 1e6 |
const real | delta = 1e-15 |
const real | eps1 = 1e-15 |
const real | sqrt_01 = 0.31622777 |
const real | limit = 1e20 |
Functions | |
static ivector | null (ivector A) |
static imatrix | transp (imatrix &A) |
real | REL (rvector A, rvector B) |
bool | TOO_BAD (ivector &A) |
void | GUESS_ZEROES (rvector &x0, rvector &x1) |
void | LSS (imatrix &A, ivector &b, ivector &Y, imatrix &MyA, vector< int > &mycolumns, intvector &mycol, bool distribute, int procs, int mypid, int &errc, int &commerrc, rmatrix &R1, rmatrix &R2, imatrix &C, int &FLAGS, ofstream &ausg) |
void | LSS (imatrix &A, ivector &b, ivector &Y, imatrix &MyA, vector< int > &mycolumns, intvector &mycol, bool distribute, int procs, int mypid, int &errc, int &commerrc, ofstream &ausg) |
void | LSS (imatrix &A, ivector &b, ivector &Y, int procs, int mypid, int &errc, int &commerrc, ofstream &ausg) |
void | LSS (imatrix &A, ivector &b, ivector &Y, int &errc, int &commerrc, ofstream &ausg) |
Variables | |
static int | m |
static int | n |
static int | dim |
Parallel Version Copyright (C) 2004-2007 Markus Grimmer
Based on: Serial Version Copyright (C) 2003 Carlos Holbig and Walter Kraemer
Parallel verified solution of interval linear systems using Rump's algorithm.
|
x1 is the new, x0 the old value of an iteration. If a component of x1 has decreased by more than a factor of zerotest, then this component is set to 0. The same is done if the sign of a component has changed.
|
|
Parallel linear system solver (square matrix expected). This version: Wrapper without local variables. Without parameters for distributed data and parallel environment. -> Distribution will be done inside the routine. -> Parallel environment will be initialized inside the routine.
|
|
Parallel linear system solver (square matrix expected). This version: Wrapper without local variables. Without parameters for distributed data -> Distribution will be done inside the routine. Includes parameters for already initialized parallel environment.
|
|
Parallel linear system solver (square matrix expected). This version: Wrapper without local variables. Includes parameters for already distributed data (matrix columns) and already initialized parallel environment
|
|
Parallel linear system solver (square matrix expected). The result of Y is an enclosure of the solution errc = 0: Y is enclosure of the solution errc = 1: no enclosure obtained, bad condition (?) errc = 2: no enclosure obtained, matrix A singular (?) FLAGS = 0: R1,R2,C have not yet been computed FLAGS = 1: only R1 has been computed FLAGS = 2: R1 and corresponding C have been computed FLAGS = 3: R1 and R2 are computedbut not the corresponding C FLAGS = 4: R1, R2 and the corresponding C have been computed This version: Includes parameters for already distributed data (matrix columns) and already initialized parallel environment
|
|
Zero vector |
|
REL computes componentwise the maximum relative error of A w.r.t B. if A[i] and B[i] do not have the same sign or if B[i] = 0, then rel. error = 0 for this component.
|
|
TOO_BAD = accuracy of A is far too bad note: 0 for false, 1 for true;
|
|
Transposition of a matrix
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|