Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F78976151
ch-05-ex-10-solution.cpp
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Sat, Aug 24, 02:13
Size
1 KB
Mime Type
text/x-c
Expires
Mon, Aug 26, 02:13 (1 d, 21 h)
Engine
blob
Format
Raw Data
Handle
20130475
Attached To
R9360 PCSC_git
ch-05-ex-10-solution.cpp
View Options
/*
* chapter-05-exercise-10.cpp
*
* Gauss with partial pivoting
*
* run the program as: ./ch-05-ex-10-solution Agauss.dat bgauss.dat
*
* Created on: Oct 11, 2012
* Author: Radu Popescu <radu.popescu@epfl.ch>
*/
#include <cassert>
#include <cmath>
#include <cstdlib>
#include <fstream>
#include <iostream>
#include "ArrayUtil.hpp"
#include "Gauss.hpp"
int main(int argc, char *argv[]) {
if (argc != 3) {
std::cout << "Missing arguments. Please run program as: "
<< " ./ch-05-ex-06-solution Agauss.dat bgauss.dat" << std::endl;
return -1;
}
const char *file_name_A = argv[1];
const char *file_name_b = argv[2];
double **A;
int m, n;
Read(file_name_A, A, m, n);
double *b;
Read(file_name_b, b, n);
double *x;
Allocate(x, n);
// Solve Ax = b
GaussPartialPivoting(A, b, n);
BackwardSolve(A, b, n, x);
std::cout << "The matrix A is: " << std::endl;
Print(A, m, n, std::cout);
std::cout << "The rhs b is: " << std::endl;
Print(b, n, std::cout, true);
std::cout << "The solution x is: " << std::endl;
Print(x, n, std::cout, true);
Delete(A, m);
Delete(b);
Delete(x);
return 0;
}
Event Timeline
Log In to Comment