Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
On Thursday, 7th July from 1 to 3 pm there will be a maintenance with a short downtime of GitLab.
Open sidebar
elpa
elpa
Commits
cadf98b1
Commit
cadf98b1
authored
Aug 22, 2017
by
Pavel Kus
Browse files
getting rid of CONST_REAL_0_5
parent
7cd787ea
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
src/elpa2/elpa2_tridiag_band_template.F90
View file @
cadf98b1
...
...
@@ -632,12 +632,11 @@
if
(
wantDebug
)
call
obj
%
timer
%
stop
(
"blas"
)
#if REALCASE == 1
x
=
dot_product
(
hv
(
1
:
nc
),
hd
(
1
:
nc
))
*
tau
hd
(
1
:
nc
)
=
hd
(
1
:
nc
)
-
CONST_0_5
*
x
*
hv
(
1
:
nc
)
#endif
#if COMPLEXCASE == 1
x
=
dot_product
(
hv
(
1
:
nc
),
hd
(
1
:
nc
))
*
conjg
(
tau
)
hd
(
1
:
nc
)
=
hd
(
1
:
nc
)
-
0.5
*
x
*
hv
(
1
:
nc
)
#endif
hd
(
1
:
nc
)
=
hd
(
1
:
nc
)
-
0.5_rk
*
x
*
hv
(
1
:
nc
)
if
(
wantDebug
)
call
obj
%
timer
%
start
(
"blas"
)
#if REALCASE == 1
call
PRECISION_SYR2
(
'L'
,
nc
,
-
ONE
,
hd
,
1
,
hv
,
1
,
ab
(
1
,
ns
),
2
*
nb
-1
)
...
...
@@ -1026,12 +1025,11 @@
! Transform diagonal block
#if REALCASE == 1
x
=
dot_product
(
hv
(
1
:
nc
),
hd
(
1
:
nc
))
*
tau
hd
(
1
:
nc
)
=
hd
(
1
:
nc
)
-
CONST_0_5
*
x
*
hv
(
1
:
nc
)
#endif
#if COMPLEXCASE == 1
x
=
dot_product
(
hv
(
1
:
nc
),
hd
(
1
:
nc
))
*
conjg
(
tau
)
hd
(
1
:
nc
)
=
hd
(
1
:
nc
)
-
CONST_REAL_0_5
*
x
*
hv
(
1
:
nc
)
#endif
hd
(
1
:
nc
)
=
hd
(
1
:
nc
)
-
0.5_rk
*
x
*
hv
(
1
:
nc
)
if
(
my_pe
>
0
.and.
iblk
==
1
)
then
! The first column of the diagonal block has to be send to the previous PE
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment