Error in parallel computation of H-AMLS #2
Description
When executing the H-AMLS method with multiple threads, in some cases it has been observed that the program is aborted due to a memory error. The error is probably caused when executing the function "solve_lower_right_symmetric" in scr/THAMLS.cc with which the transformed mass matrix in H-AMLS is computed. The first guess is, that when computing with multiple threads the usage of the subroutine "EIGEN_MISC::delete_upper_block_triangular()" in "solve_lower_right_symmetric()" may cause that error when in the next step of the routine the HLIBpro function "solve_lower_right()" is started for further processing. Hence, the function "EIGEN_MISC::delete_upper_block_triangular()" should be modified that it is inline with the parallel HLIBpro arithmetic function "solve_lower_right()". Using an HLIBpro version older than v2.7 this error did not arise.