Commit 19a7247e by Pavel Kus Committed by Andreas Marek

### removing forgotten timing output

parent 48931833
 ... ... @@ -736,7 +736,6 @@ ! Compute eigenvectors of the rank-1 modified matrix. ! Parts for multiplying with upper half of Q: call obj%timer%start("strange_loop") do i = 1, ncnt j = idx(idxq1(i+ns)) ! Calculate the j-th eigenvector of the deflated system ... ... @@ -747,7 +746,6 @@ &(obj,tmp,nnzu,ddiff(j)) ev(1:nnzu,i) = zu(1:nnzu) / tmp(1:nnzu) * ev_scale(j) enddo call obj%timer%stop("strange_loop") if(useGPU) then !TODO: it should be enough to copy l_rows x ncnt ... ... @@ -801,7 +799,6 @@ ! Compute eigenvectors of the rank-1 modified matrix. ! Parts for multiplying with lower half of Q: call obj%timer%start("strange_loop") do i = 1, ncnt j = idx(idxq1(i+ns)) ! Calculate the j-th eigenvector of the deflated system ... ... @@ -812,7 +809,6 @@ &(obj,tmp,nnzl,ddiff(j)) ev(1:nnzl,i) = zl(1:nnzl) / tmp(1:nnzl) * ev_scale(j) enddo call obj%timer%stop("strange_loop") if(useGPU) then !TODO the previous loop could be possible to do on device and thus ... ...
 ... ... @@ -278,24 +278,24 @@ thread_width = (l_nev-1)/max_threads + 1 ! number of eigenvectors per OMP thread #if REALCASE == 1 call obj%get("stripe_width_real",stripe_width, error) ! call obj%get("stripe_width_real",stripe_width, error) #ifdef DOUBLE_PRECISION_REAL !stripe_width = 48 ! Must be a multiple of 4 stripe_width = 48 ! Must be a multiple of 4 #else stripe_width = stripe_width * 2 !stripe_width = 96 ! Must be a multiple of 8 !stripe_width = stripe_width * 2 stripe_width = 96 ! Must be a multiple of 8 #endif #endif /* REALCASE */ #if COMPLEXCASE == 1 call obj%get("stripe_width_complex",stripe_width, error) !call obj%get("stripe_width_complex",stripe_width, error) #ifdef DOUBLE_PRECISION_COMPLEX !stripe_width = 48 ! Must be a multiple of 2 stripe_width = 48 ! Must be a multiple of 2 #else stripe_width = stripe_width * 2 !stripe_width = 48 ! Must be a multiple of 4 !stripe_width = stripe_width * 2 stripe_width = 48 ! Must be a multiple of 4 #endif #endif /* COMPLEXCASE */ ... ...
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!