Angular Momentum Coupling Coefficients
Edit Page
Report
Scan day: 03 March 2014 UTC
49
Virus safety - good
Description: Fortran 90 code by Paul Stevenson.
! anglib.f90: angular momentum coupling coefficients in Fortran 90 ! Copyright (C) 1998 Paul Stevenson ! ! This library is free software; you can redistribute it and/or ! modify it under the terms of the GNU Lesser General Public ! License as published by the Free Software Foundation; either ! version 2.1 of the License, or (at your option) any later version. ! ! This library is distributed in the hope that it will be useful, ! but WITHOUT ANY WARRANTY; without even the implied warranty of ! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU ! Lesser General Public License for more details. ! ! You should have received a copy of the GNU Lesser General Public ! License along with this library; if not, write to the Free Software ! Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA ! 02110-1301,USA module anglib ! Library of angular momentum coupling coefficient routines in fortran 90 ! Paul Stevenson, [email protected] integer, parameter :: rk = selected_real_kind(p=15) contains function cleb(j1,m1,j2,m2,j,m) implicit none ! calculate a clebsch-gordan coefficient ! arguments are integer and twice the true value. real(rk) :: cleb,factor,sum integer :: j1,m1,j2,m2,j,m,par,z,zmin,zmax ! some checks for validity (let's just return zero for bogus arguments) if (2*(j1/2)-int(2*(j1/2.0)) /= 2*(abs(m1)/2)-int(2*(abs(m1)/2.0)) .or. & 2*(j2/2)-int(2*(j2/2.0)) /= 2*(abs(m2)/2)-int(2*(abs(m2)/2.0)) .or. & 2*(j/2)-int(2*(j/2.0)) /= 2*(abs(m)/2)-int(2*(abs(m)/2.0)) .or. & j1j1 .or. abs(m2)>j2 .or.& abs(m)>j .or. j1+j2
Size: 1610 chars
Contact Information
Email: —
Phone&Fax: —
Address: —
Extended: —
WEBSITE Info
| Page title: | |
| Keywords: | |
| Description: | |
| IP-address: | 131.227.115.98 |