a*s
2 楼
【 以下文字转载自 History 讨论区 】
发信人: arabianights (kd), 信区: History
标 题: 虚云 大師活了119歲,真的假的?
发信站: BBS 未名空间站 (Wed Jan 4 12:41:38 2012, 美东)
http://zh.wikipedia.org/wiki/%E9%87%8A%E8%99%9A%E4%BA%91
如果是真的,應該是古今中外最長壽之人了吧
发信人: arabianights (kd), 信区: History
标 题: 虚云 大師活了119歲,真的假的?
发信站: BBS 未名空间站 (Wed Jan 4 12:41:38 2012, 美东)
http://zh.wikipedia.org/wiki/%E9%87%8A%E8%99%9A%E4%BA%91
如果是真的,應該是古今中外最長壽之人了吧
p*u
3 楼
大家说说可能是啥?
amd 4000+ 系统刚装的win 7
内存双通道 ddr2 2*2gb + 2*1gb
显卡 是mc 买的 最便宜的pcie显卡。
硬盘旧的 300gb sata 2 硬盘。 刚刚磁盘整理过的。
大家给点意见吧。谢谢。
amd 4000+ 系统刚装的win 7
内存双通道 ddr2 2*2gb + 2*1gb
显卡 是mc 买的 最便宜的pcie显卡。
硬盘旧的 300gb sata 2 硬盘。 刚刚磁盘整理过的。
大家给点意见吧。谢谢。
s*n
4 楼
搞得多少人浪费多少时间了
z*7
5 楼
为了实现并行计算,怎么在原来的Fortran 95串行代码中插入OpenMP?
3X!
3X!
q*n
7 楼
多修神足的人, 欲住世一劫,也是可以的.
【在 a**********s 的大作中提到】
: 【 以下文字转载自 History 讨论区 】
: 发信人: arabianights (kd), 信区: History
: 标 题: 虚云 大師活了119歲,真的假的?
: 发信站: BBS 未名空间站 (Wed Jan 4 12:41:38 2012, 美东)
: http://zh.wikipedia.org/wiki/%E9%87%8A%E8%99%9A%E4%BA%91
: 如果是真的,應該是古今中外最長壽之人了吧
【在 a**********s 的大作中提到】
: 【 以下文字转载自 History 讨论区 】
: 发信人: arabianights (kd), 信区: History
: 标 题: 虚云 大師活了119歲,真的假的?
: 发信站: BBS 未名空间站 (Wed Jan 4 12:41:38 2012, 美东)
: http://zh.wikipedia.org/wiki/%E9%87%8A%E8%99%9A%E4%BA%91
: 如果是真的,應該是古今中外最長壽之人了吧
m*n
10 楼
看编译器吧。for gnu fortran: http://gcc.gnu.org/onlinedocs/gfortran/OpenMP.html
l*h
11 楼
牛牛应该不会喜欢最贵的牌子的钻戒吧?比如一个很好牌子的小颗钻戒和一个没什么牌
子的大颗钻戒,牛牛选哪个啊?
子的大颗钻戒,牛牛选哪个啊?
z*7
14 楼
这段代码如何用openmp并行运算?(附件是全模块) 并保证不溢出16gb内存
全套程序下载地址如下:
http://www.atomic-theory.uni-jen ... ange/ratip-2012.tgz
subroutine auger_calculate_amplitudes()
!--------------------------------------------------------------------
! calculates for all transitions in turn the required continuum
! spinors and auger amplitudes.
!
! calls: add_csf_to_basis(), anco_calculate_csf_matrix(),
! auger_channel_amplitude(), auger_pure_matrix(),
! auger_transition_properties() cowf_iterate_csp(),
! cowf_set_xk_coefficients(), cowf_set_yk_coefficients(),
! print_configuration_scheme()
! set_configuration_scheme().
!--------------------------------------------------------------------
!
integer :: i, j, n, nw, nocsf
real(kind=dp) :: energy
type(nkappa) :: subshell
integer, dimension(, allocatable :: ndx
!
n = asf_final%csf_set%nocsf + asf_initial%csf_set%nocsf
!
! allocate for a "first time"; it is first dellocated before any usage
allocate( auger_csp%p(1:n_grasp92), auger_csp%q(1:n_grasp92) )
allocate( cowf_csp%p(1:10), cowf_csp%q(1:10) )
allocate( ndx(1:n) )
do i = 1,number_of_transitions
if (transition(i)%energy < zero) then
transition(i)%probability = zero
transition(i)%alpha_2 = zero
transition(i)%alpha_4 = zero
transition(i)%eta_2 = zero
transition(i)%eta_4 = zero
cycle
end if
!
do j = 1,transition(i)%number_of_channels
energy = transition(i)%energy
!!x print *, "auger_calculate_amplitudes - a"
call set_configuration_scheme(asf_final%csf_set,asf_cont%csf_set
,&
-1,transition(i)%channel(j)%kappa,
&
transition(i)%totalj_f,transition(i)%parity_f,
&
transition(i)%totalj_i,transition(i)%parity_i,
&
append=.false.,index=ndx)
!
auger%no_f = asf_cont%csf_set%nocsf
allocate( auger%ndx_f(auger%no_f) )
auger%ndx_f(1:auger%no_f) = ndx(1:auger%no_f)
!
nw = asf_cont%csf_set%nwshells
if (rabs_use_stop .and. nw /= asf_final%csf_set%nwshells + 1)
then
stop "auger_calculate_amplitudes(): program stop a."
end if
! calculate the mcp coefficients for the current coupling scheme
! as well as the d_rs, y_k(ab), and x_k(abcd) coefficients
nocsf = asf_cont%csf_set%nocsf
call anco_calculate_csf_matrix(asf_cont%csf_set,1,nocsf,1,nocsf)
call cowf_set_drs_coefficients(transition(i)%asff, &
asf_cont%csf_set,ndx)
subshell = nkappa(-1,transition(i)%channel(j)%kappa)
call cowf_set_yk_coefficients(subshell,asf_cont%csf_set)
!!x print *, "auger_calculate_amplitudes - f"
call cowf_set_xk_coefficients(subshell,asf_cont%csf_set)
!!x print *, "auger_calculate_amplitudes - g"
!
! now iterate the continuum spinors for this channel
! cowf_solve_homogeneous_eqn = .true.
cowf_start_homogeneous = .true.
cowf_phaseshift_wkb = .true.
cowf_phaseshift_zero_potential = .false.
cowf_phaseshift_coulomb = .false.
!
!! cowf_norm_nonrel = .true.
cowf_norm_wkb = .true.
call cowf_iterate_csp(energy,subshell)
!
auger_csp = cowf_csp
transition(i)%channel(j)%phase = auger_csp%phase
!
! define the 'extended' configuration scheme for calculating
! the auger matrix and allocate memory
call add_csf_to_basis(asf_initial%csf_set,asf_cont%csf_set,
&
transition(i)%totalj_i,transition(i)%parity_i,index=ndx)
if (auger_print_csf_scheme) then
call print_configuration_scheme(6,asf_cont%csf_set)
end if
!
auger%no_i = asf_cont%csf_set%nocsf - auger%no_f
allocate( auger%ndx_i(auger%no_i) )
auger%ndx_i(1:auger%no_i) = ndx(1+auger%no_f:asf_cont%csf_set%
nocsf)
allocate( auger%matrix(1:auger%no_f,1:auger%no_i) )
!
! calculate the 'pure' auger matrix in the given csf scheme
! (not including mixing coefficients)
call auger_pure_matrix(asf_cont%csf_set,i)
!
call auger_channel_amplitude(i,j)
!
deallocate( auger%ndx_f, auger%ndx_i, auger%matrix )
call deallocate_csf_basis(asf_cont%csf_set)
end do
!
! calculates all selected properties for the selected transition
call auger_transition_properties(transition(i))
end do
deallocate( ndx, auger_csp%p, auger_csp%q)
!
end subroutine auger_calculate_amplitudes
!
!
subroutine auger_channel_amplitude(i,j)
!--------------------------------------------------------------------
! calculates the auger amplitude of channel j of transition i
! by summing over the 'pure' auger matrix using the proper weights of
! transition i.
!
! calls:
!--------------------------------------------------------------------
!
integer, intent(in) :: i, j
!
integer :: asfi, asff, l, r, rr, s, ss
real(kind=dp) :: phase, value
!
if (auger_print_main_csf_me) then
print *, " "
print *, "main contribution from initial- and final-state csf "// &
"(abs(c_i*c_f) > 0.01)"
print *, "----------------------------------------------------"// &
"---------------------"
print *, " "
print *, " i-csf f-csf kappa c_i c_f c_i*c_f"/
/&
" c_i*c_f*a_if "
print *, "-------------------------------------------------------"/
/&
"-----------------"
end if
!
asfi = transition(i)%asfi; asff = transition(i)%asff
value = zero
do r = 1,auger%no_f
rr = auger%ndx_f(r)
do s = 1,auger%no_i
ss = auger%ndx_i(s)
value = value + asf_final%asf(asff)%eigenvector(rr) * &
auger%matrix(r,s) * asf_initial%asf(asfi)%eigenvector(ss)
!
if (auger_print_main_csf_me) then
if (abs(asf_final%asf(asff)%eigenvector(rr)*
&
asf_initial%asf(asfi)%eigenvector(ss)) > 0.01_dp .
and.&
abs(asf_final%asf(asff)%eigenvector(rr) *
&
auger%matrix(r,s) *
&
asf_initial%asf(asfi)%eigenvector(ss)) > 0.000001_dp)
then
!
! determine first the radial integrals from the occupation
! of the csf
write(*,1) ss,rr, &
orbital_symmetry(transition(i)%channel(j)%kappa), &
asf_initial%asf(asfi)%eigenvector(ss), &
asf_final%asf(asff)%eigenvector(rr), &
asf_final%asf(asff)%eigenvector(rr)* &
asf_initial%asf(asfi)%eigenvector(ss), &
asf_final%asf(asff)%eigenvector(rr)* &
asf_initial%asf(asfi)%eigenvector(ss)* &
auger%matrix(r,s)
1 format(1x,i7,i10,6x,a2,3x,f6.3,2x,f6.3,5x,f8.5,4x,f9.6)
end if
end if
!
end do
end do
!
if (auger_print_main_csf_me) then
print *, "-------------------------------------------------------"/
/&
"-----------------"
end if
!
l = angular_momentum_l(transition(i)%channel(j)%kappa)
phase = transition(i)%channel(j)%phase
!
transition(i)%channel(j)%amplitude_re = value
transition(i)%channel(j)%amplitude = cmplx(zero,one)**l *
&
exp( -cmplx(zero,one)*phase) * cmplx(value,
zero)
!
print *, "i,j,transition(i)%channel(j)%amplitude = ", &
i,j,transition(i)%channel(j)%amplitude
!
end subroutine auger_channel_amplitude
!
全套程序下载地址如下:
http://www.atomic-theory.uni-jen ... ange/ratip-2012.tgz
subroutine auger_calculate_amplitudes()
!--------------------------------------------------------------------
! calculates for all transitions in turn the required continuum
! spinors and auger amplitudes.
!
! calls: add_csf_to_basis(), anco_calculate_csf_matrix(),
! auger_channel_amplitude(), auger_pure_matrix(),
! auger_transition_properties() cowf_iterate_csp(),
! cowf_set_xk_coefficients(), cowf_set_yk_coefficients(),
! print_configuration_scheme()
! set_configuration_scheme().
!--------------------------------------------------------------------
!
integer :: i, j, n, nw, nocsf
real(kind=dp) :: energy
type(nkappa) :: subshell
integer, dimension(, allocatable :: ndx
!
n = asf_final%csf_set%nocsf + asf_initial%csf_set%nocsf
!
! allocate for a "first time"; it is first dellocated before any usage
allocate( auger_csp%p(1:n_grasp92), auger_csp%q(1:n_grasp92) )
allocate( cowf_csp%p(1:10), cowf_csp%q(1:10) )
allocate( ndx(1:n) )
do i = 1,number_of_transitions
if (transition(i)%energy < zero) then
transition(i)%probability = zero
transition(i)%alpha_2 = zero
transition(i)%alpha_4 = zero
transition(i)%eta_2 = zero
transition(i)%eta_4 = zero
cycle
end if
!
do j = 1,transition(i)%number_of_channels
energy = transition(i)%energy
!!x print *, "auger_calculate_amplitudes - a"
call set_configuration_scheme(asf_final%csf_set,asf_cont%csf_set
,&
-1,transition(i)%channel(j)%kappa,
&
transition(i)%totalj_f,transition(i)%parity_f,
&
transition(i)%totalj_i,transition(i)%parity_i,
&
append=.false.,index=ndx)
!
auger%no_f = asf_cont%csf_set%nocsf
allocate( auger%ndx_f(auger%no_f) )
auger%ndx_f(1:auger%no_f) = ndx(1:auger%no_f)
!
nw = asf_cont%csf_set%nwshells
if (rabs_use_stop .and. nw /= asf_final%csf_set%nwshells + 1)
then
stop "auger_calculate_amplitudes(): program stop a."
end if
! calculate the mcp coefficients for the current coupling scheme
! as well as the d_rs, y_k(ab), and x_k(abcd) coefficients
nocsf = asf_cont%csf_set%nocsf
call anco_calculate_csf_matrix(asf_cont%csf_set,1,nocsf,1,nocsf)
call cowf_set_drs_coefficients(transition(i)%asff, &
asf_cont%csf_set,ndx)
subshell = nkappa(-1,transition(i)%channel(j)%kappa)
call cowf_set_yk_coefficients(subshell,asf_cont%csf_set)
!!x print *, "auger_calculate_amplitudes - f"
call cowf_set_xk_coefficients(subshell,asf_cont%csf_set)
!!x print *, "auger_calculate_amplitudes - g"
!
! now iterate the continuum spinors for this channel
! cowf_solve_homogeneous_eqn = .true.
cowf_start_homogeneous = .true.
cowf_phaseshift_wkb = .true.
cowf_phaseshift_zero_potential = .false.
cowf_phaseshift_coulomb = .false.
!
!! cowf_norm_nonrel = .true.
cowf_norm_wkb = .true.
call cowf_iterate_csp(energy,subshell)
!
auger_csp = cowf_csp
transition(i)%channel(j)%phase = auger_csp%phase
!
! define the 'extended' configuration scheme for calculating
! the auger matrix and allocate memory
call add_csf_to_basis(asf_initial%csf_set,asf_cont%csf_set,
&
transition(i)%totalj_i,transition(i)%parity_i,index=ndx)
if (auger_print_csf_scheme) then
call print_configuration_scheme(6,asf_cont%csf_set)
end if
!
auger%no_i = asf_cont%csf_set%nocsf - auger%no_f
allocate( auger%ndx_i(auger%no_i) )
auger%ndx_i(1:auger%no_i) = ndx(1+auger%no_f:asf_cont%csf_set%
nocsf)
allocate( auger%matrix(1:auger%no_f,1:auger%no_i) )
!
! calculate the 'pure' auger matrix in the given csf scheme
! (not including mixing coefficients)
call auger_pure_matrix(asf_cont%csf_set,i)
!
call auger_channel_amplitude(i,j)
!
deallocate( auger%ndx_f, auger%ndx_i, auger%matrix )
call deallocate_csf_basis(asf_cont%csf_set)
end do
!
! calculates all selected properties for the selected transition
call auger_transition_properties(transition(i))
end do
deallocate( ndx, auger_csp%p, auger_csp%q)
!
end subroutine auger_calculate_amplitudes
!
!
subroutine auger_channel_amplitude(i,j)
!--------------------------------------------------------------------
! calculates the auger amplitude of channel j of transition i
! by summing over the 'pure' auger matrix using the proper weights of
! transition i.
!
! calls:
!--------------------------------------------------------------------
!
integer, intent(in) :: i, j
!
integer :: asfi, asff, l, r, rr, s, ss
real(kind=dp) :: phase, value
!
if (auger_print_main_csf_me) then
print *, " "
print *, "main contribution from initial- and final-state csf "// &
"(abs(c_i*c_f) > 0.01)"
print *, "----------------------------------------------------"// &
"---------------------"
print *, " "
print *, " i-csf f-csf kappa c_i c_f c_i*c_f"/
/&
" c_i*c_f*a_if "
print *, "-------------------------------------------------------"/
/&
"-----------------"
end if
!
asfi = transition(i)%asfi; asff = transition(i)%asff
value = zero
do r = 1,auger%no_f
rr = auger%ndx_f(r)
do s = 1,auger%no_i
ss = auger%ndx_i(s)
value = value + asf_final%asf(asff)%eigenvector(rr) * &
auger%matrix(r,s) * asf_initial%asf(asfi)%eigenvector(ss)
!
if (auger_print_main_csf_me) then
if (abs(asf_final%asf(asff)%eigenvector(rr)*
&
asf_initial%asf(asfi)%eigenvector(ss)) > 0.01_dp .
and.&
abs(asf_final%asf(asff)%eigenvector(rr) *
&
auger%matrix(r,s) *
&
asf_initial%asf(asfi)%eigenvector(ss)) > 0.000001_dp)
then
!
! determine first the radial integrals from the occupation
! of the csf
write(*,1) ss,rr, &
orbital_symmetry(transition(i)%channel(j)%kappa), &
asf_initial%asf(asfi)%eigenvector(ss), &
asf_final%asf(asff)%eigenvector(rr), &
asf_final%asf(asff)%eigenvector(rr)* &
asf_initial%asf(asfi)%eigenvector(ss), &
asf_final%asf(asff)%eigenvector(rr)* &
asf_initial%asf(asfi)%eigenvector(ss)* &
auger%matrix(r,s)
1 format(1x,i7,i10,6x,a2,3x,f6.3,2x,f6.3,5x,f8.5,4x,f9.6)
end if
end if
!
end do
end do
!
if (auger_print_main_csf_me) then
print *, "-------------------------------------------------------"/
/&
"-----------------"
end if
!
l = angular_momentum_l(transition(i)%channel(j)%kappa)
phase = transition(i)%channel(j)%phase
!
transition(i)%channel(j)%amplitude_re = value
transition(i)%channel(j)%amplitude = cmplx(zero,one)**l *
&
exp( -cmplx(zero,one)*phase) * cmplx(value,
zero)
!
print *, "i,j,transition(i)%channel(j)%amplitude = ", &
i,j,transition(i)%channel(j)%amplitude
!
end subroutine auger_channel_amplitude
!
S*P
17 楼
牌子绝对加分
l*a
19 楼
同意简单经典。
我喜欢金的,有美美和我一样的没?
我喜欢金的,有美美和我一样的没?
c*g
21 楼
harry winston
适合金牛⋯⋯哈哈哈哈
适合金牛⋯⋯哈哈哈哈
i*n
23 楼
金牛座mm喜欢什么款式的男人啊?
h*2
26 楼
戒指送出去了,mm很高兴。但是同时也很谨慎……说因为我说是birthday gift才收下
的。她还说我们以后要多一起出去一起见面一起玩~
嗯,基本就这些了~
大家帮忙分析分析~
的。她还说我们以后要多一起出去一起见面一起玩~
嗯,基本就这些了~
大家帮忙分析分析~
相关阅读
C++能否重定义函数?刚刚在茶水间看电视expedia暴涨Java 这是啥语法?为什么你会讨厌百度公司?请教:如何使用浏览器的cookies给PowerQuery无标题现在大规模ML的训练网络还是瓶颈了吗NLP sentiment analysis 有推荐吗[bssd] 我经常过几年就会被很多人围着一起反对抖音,再见!国家正式发声!python 没有替代品bokeh vs Dash, 哪个好?小米上市就是雷军人设崩塌的开始。。。教小孩学Java用什么?寻合伙人/Chief Scientist (computer vision方向)我遇到的最不舒服的语言是Stata从刘强东的“不卖假货”想到了中国的互联网行业中国政府派黑客混入开源代码项目偷偷部署安全漏洞1D to 2D conversion problem本科刷acm好,还是刷文章好