Skip to content

w06-g02-CGsolver-class

Conjugate Gradient Solver (CGsolver) Class

Objective

The overall objective of this assignment is to implement the Conjugate Gradient Method in a class CGsolver. This solver will later be part of the class project Direct Stiffness Method.

The code for this class needs to be split into a header file .hpp and its implementation in the .cpp file. Most of the interface of the CGsolver class is already implemented in tehpc/src/conjugate_gradient_solver.hpp, and is used in a simple test tehpc/tests/CGsolver_test.cpp. You have to implement the method solve of the CGsolver class in tehpc/src/conjugate_gradient_solver.cpp.

Please recycle the code from last week tehpc/examples/conjugate_gradient_function.cpp (method described in details in #2 (closed)).

Implementation

  1. load appropriate modules by performing the following command in your terminal: source tehpc/tools/tehpc_load_modules.sh.

  2. If conjugate_gradient_solver.cpp is not already included at the appropriate place in tehpc/src/CMakeLists.txt, add it. This makes sure that it is compiled.

  3. Add tehpc_add_test(CGsolver_test) in the tehpc/tests/CMakeLists.txt to include the CGsolver_test in compilation.

  4. Compile your code, use the following steps (this is new):

    a. go into tehpc/build/

    b. type: either cmake -DCMAKE_BUILD_TYPE:STRING=Debug -DTEHPC_TESTS:STRING=ON .. or ccmake .. and make sure that build type is Debug and Tests are ON (you only have to do this the very first time).

    c. type: make (do this whenever you want to recompile).

  5. Correct all compilation errors. Use code from tehpc/basics/conjugate_gradient_function.cpp and the provided operator overloading (see tehpc/src/matrix.hpp, tehpc/src/vector.hpp and others) to implement tehpc/src/conjugate_gradient_solver.cpp. Important: you are not allowed to change the CGsolver_test.cpp file.

  6. If the code compiles without errors, run the test:

    a. go into tehpc/build/tests/

    b. type: ./CGsolver_test

  7. Make sure that the code runs without error and check printed message to see if it went to the end. Correct errors if needed. Again, you are not allowed to change the CGsolver_test.cpp files.

  8. run valgrind to check for memory leakage: valgrind --leak-check=full ./CGsolver_test.

  9. correct memory leakage

  10. Make sure that your code is well structured and commented.