repo.or.cz
/
official-gcc.git
/
blob
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
log
|
graphiclog1
|
graphiclog2
|
commit
|
commitdiff
|
tree
|
refs
|
edit
|
fork
blame
|
history
|
raw
|
HEAD
libgfortran/ChangeLog:
[official-gcc.git]
/
gcc
/
testsuite
/
gfortran.dg
/
graphite
/
id-18.f90
blob
273e670fcee82c575ac0ca40f2f9dbb48016c706
1
MODULE
spherical_harmonics
2
INTEGER
,
PARAMETER
::
dp
=
SELECTED_REAL_KIND
(
14
,
200
)
3
CONTAINS
4
FUNCTION
dlegendre
(
x
,
l
,
m
)
RESULT
(
dplm
)
5
SELECT
CASE
(
l
)
6
CASE
(
0
)
7
dplm
=
0.0
_dp
8
CASE
(
1
)
9
dplm
=
1.0
_dp
10
CASE DEFAULT
11
IF
(
mm
>
0
)
THEN
12
dpmm
= -
m
13
DO
im
=
1
,
mm
14
dpmm
= -
dpmm
15
END DO
16
IF
(
l
==
mm
+
1
)
THEN
17
DO
il
=
mm
+
2
,
l
18
dpll
=
dpmm
19
END DO
20
dplm
=
dpll
21
END IF
22
END IF
23
END
SELECT
24
END FUNCTION
dlegendre
25
END MODULE
spherical_harmonics