²é¿´: 1691  |  »Ø¸´: 4

superrice

ľ³æ (ÕýʽдÊÖ)

[ÇóÖú] Fortran90 ±àÒ벻ͨ¹ý ÒÑÓÐ2È˲ÎÓë

ÇëÎÊ´óÉñÊÇʲôÇé¿ö£¬»¹ÓÐÔõôÄܰÑÄǸöÌÖÑáµÄĬÈϱíÇé¿ì½Ý¼ü¸øÈ¥µô°¡

subroutine laplacian(f, lapl)

  implicit none

  real(8), intent(in)  :: f(N, N)    ! function whose laplacian is to be taken
  real(8), intent(out) :: lapl(N, N) ! the laplacian


  lapl = 0.0_8


  integer :: ix, iy, k

  integer, parameter :: order = 4
  real(8), allocatable :: c(


  allocate(c(-orderrder))
  c(-orderrder) = (/ &  ! The coefficients of the laplacian...
      -1.785714d-3, 2.539683d-2, -0.2d0, 1.6d0,      &
      -2.847222d0,                                   &
      1.6d0, -0.2d0, 2.539683d-2, -1.785714d-3 /)


  do ix = 1, N
    do iy = 1, N
      lapl(ix, iy) = (0.0_8, 0.0_8)
      do k = -order, order
        if(iy+k>=1 .and. iy+k<=N) lapl(ix, iy) = lapl(ix, iy) + c(k)*f(ix, iy+k)
        if(ix+k>=1 .and. ix+k<=N) lapl(ix, iy) = lapl(ix, iy) + c(k)*f(ix+k, iy)
      end do
    end do
  end do

  lapl = lapl/delta**2


´íÎó£º
integer :: ix, iy, k
                      1
Error: Unexpected data declaration statement at (1)
mesh.f90:129.33:

  integer, parameter :: order = 4
                                 1
Error: Unexpected data declaration statement at (1)
mesh.f90:130.30:

  real(8), allocatable :: c(
                              1
Error: Unexpected data declaration statement at (1)
mesh.f90:133.12:

  allocate(c(-orderrder))
            1
Error: Syntax error in ALLOCATE statement at (1)
mesh.f90:134.2:

  c(-orderrder) = (/ &  ! The coefficients of the laplacian...
  1
Error: Unclassifiable statement at (1)
mesh.f90:181.22:

  integer :: ix, iy, k
                      1
Error: Unexpected data declaration statement at (1)
mesh.f90:183.33:

  integer, parameter :: order = 4
                                 1
Error: Unexpected data declaration statement at (1)
mesh.f90:184.30:

  real(8), allocatable :: c(
                              1
Error: Unexpected data declaration statement at (1)
mesh.f90:186.12:

  allocate(c(-orderrder))
            1
Error: Syntax error in ALLOCATE statement at (1)
mesh.f90:187.2:

  c(-orderrder) = (/ &  ! The coefficients of the laplacian...
  1
Error: Unclassifiable statement at (1)
mesh.f90:192.7:

  do ix = 1, N
       1
Error: Symbol 'ix' at (1) has no IMPLICIT type
mesh.f90:193.9:

    do iy = 1, N
         1
Error: Symbol 'iy' at (1) has no IMPLICIT type
mesh.f90:195.10:

      do k = -order, order
          1
Error: Symbol 'k' at (1) has no IMPLICIT type
mesh.f90:195.19:

      do k = -order, order
                   1
Error: Symbol 'order' at (1) has no IMPLICIT type
mesh.f90:140.7:

  do ix = 1, N
       1
Error: Symbol 'ix' at (1) has no IMPLICIT type
mesh.f90:141.9:

    do iy = 1, N
         1
Error: Symbol 'iy' at (1) has no IMPLICIT type
mesh.f90:143.10:

      do k = -order, order
          1
Error: Symbol 'k' at (1) has no IMPLICIT type
mesh.f90:143.19:

      do k = -order, order
                   1
Error: Symbol 'order' at (1) has no IMPLICIT type
make[2]: *** [mesh.o] Error 1
make[2]: Leaving directory `/public/home/qishengzhang/chaowang/script/qd-0.1.0/src'
make[1]: *** [all-recursive] Error 1
make[1]: Leaving directory `/public/home/qishengzhang/chaowang/script/qd-0.1.0'
make: *** [all] Error 2
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

δ°¶ÍöÁé

Ìú¸Ëľ³æ (ÖøÃûдÊÖ)

ÔÝʱÀ뿪Сľ³æ


¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
Äã¿ÉÒÔÉÏ´«¸½¼þ£¬°ïÄã¿´¿´¡£Õâ¸ö±íÇéÒ²Õæ¶º
2Â¥2016-09-26 15:22:18
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
CODE:
...

3Â¥2016-09-26 16:00:26
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

Â©ÍøÖ®Óæ

Ìú³æ (³õÈëÎÄ̳)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

N,deltaû¶¨Ò壬Óï¾ä˳ÐòÎÊÌ⣺

subroutine laplacian(f,lapl,N,delta)

   implicit none
   integer :: N,ix,iy,k
   real(8):: delta
   real(8),intent(in)  :: f(N,N)    ! function whose laplacian is to be taken
   real(8),intent(out) :: lapl(N,N) ! the laplacian

   integer,parameter :: order = 4
   real(8),allocatable:: C(

   lapl = 0.0_8

   allocate(c(order+order+1))

   c(-orderrder) = (/ &  ! The coefficients of the laplacian...
       -1.785714d-3,2.539683d-2,-0.2d0,1.6d0,     &
       -2.847222d0,                                  &
       1.6d0,-0.2d0,2.539683d-2,-1.785714d-3 /)


   do ix = 1,N
     do iy = 1,N
       lapl(ix,iy) = (0.0_8,0.0_8)
       do k = -order,order
         if(iy+k>=1 .and. iy+k<=N) lapl(ix,iy) = lapl(ix,iy) + c(k)*f(ix,iy+k)
         if(ix+k>=1 .and. ix+k<=N) lapl(ix,iy) = lapl(ix,iy) + c(k)*f(ix+k,iy)
       end do
     end do
   end do

   lapl = lapl/delta**2

end
4Â¥2016-10-02 11:25:05
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

Â©ÍøÖ®Óæ

Ìú³æ (³õÈëÎÄ̳)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

ϵͳ°Ñ"£º"ºÍ""£¬ÒÔ¼°"£º"Óë"o" ºÏ²¢³É±íÇé·ûºÅÁË£¬ÔÚðºÅºóÃæ¼ÓÉϿոñ¾Í¿É±ÜÃâ¡£
subroutine laplacian(f,lapl,N,delta)

   implicit none
   integer :: N,ix,iy,k
   real(8):: delta
   real(8),intent(in)  :: f(N,N)    ! function whose laplacian is to be taken
   real(8),intent(out) :: lapl(N,N) ! the laplacian

   integer,parameter :: order = 4
   real(8),allocatable:: C(: )

   lapl = 0.0_8

   allocate(c(order+order+1))

   c(-order: order) = (/ &  ! The coefficients of the laplacian...
       -1.785714d-3,2.539683d-2,-0.2d0,1.6d0,     &
       -2.847222d0,                                  &
       1.6d0,-0.2d0,2.539683d-2,-1.785714d-3 /)


   do ix = 1,N
     do iy = 1,N
       lapl(ix,iy) = (0.0_8,0.0_8)
       do k = -order,order
         if(iy+k>=1 .and. iy+k<=N) lapl(ix,iy) = lapl(ix,iy) + c(k)*f(ix,iy+k)
         if(ix+k>=1 .and. ix+k<=N) lapl(ix,iy) = lapl(ix,iy) + c(k)*f(ix+k,iy)
       end do
     end do
   end do

   lapl = lapl/delta**2

end
5Â¥2016-10-02 11:39:58
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ superrice µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[»ù½ðÉêÇë] Сľ³æÃ»ÂäÁË£¬³ýÁËÆíµ»Ìû×Ó£¬¼¸ºõ¿´²»µ½ÓмÛÖµµÄÌû×Ó +14 botar 2026-07-14 15/750 2026-07-19 00:15 by beefly
[»ù½ðÉêÇë] ÃæÉϵÄË®ºÜÉî +13 ²¼²¼ºÍÒ»¶þ 2026-07-18 18/900 2026-07-18 20:37 by ÈÕÔÂÌìµØ
[»ù½ðÉêÇë] ¹úÉç¿ÆÍ¨Ñ¶ÆÀÉó +4 DZˮ¹Ú¾ü123 2026-07-18 5/250 2026-07-18 19:58 by chenmodeyu
[»ù½ðÉêÇë] 5¸ö%2F£¬·Ç¹Ì¶¨¶Î22¸ö×Öĸ +8 jaasxulu 2026-07-15 9/450 2026-07-18 17:40 by dxcgbm
[»ù½ðÉêÇë] ÉúÃü¿ÚC13ÃæÉÏ»áÆÀʱ¼ä +7 luckinging 2026-07-14 11/550 2026-07-17 20:43 by ³æ²·1988
[½Ìʦ֮¼Ò] ÄÚÐÄØÑ·¦ +16 Ë®±ùÔÂÔÂÒ°Íà 2026-07-12 16/800 2026-07-17 19:25 by otani
[»ù½ðÉêÇë] ÉúÃü¿ÚÃæÉÏ»áÆÀ½áÊøÁËÂð´ó¼ÒÓÐÖªµÀ½á¹ûÂð£¿ÏÖÔÚ»¹Ã»Óнá¹ûÕý³£²»£¿ +6 ±ÌË®00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[ÂÛÎÄͶ¸å] ͶӢÎÄÆÚ¿¯ÈçºÎ²éÖØ 10+5 ²»ÖªµÀȡʲôÃû× 2026-07-15 6/300 2026-07-17 09:07 by bobvan
[»ù½ðÉêÇë] ²»ÒªÔÙÊý¹ú×ÔÈ»ÉêÇëÊéµÄ filecode µÄ·Ö¸ô·û¸öÊýÁË +14 sadpencil 2026-07-12 21/1050 2026-07-17 00:15 by wangze12014
[Óлú½»Á÷] ²®°·Éϼ׻ùÇóÖú 100+4 002tmumu 2026-07-14 7/350 2026-07-16 19:47 by 002tmumu
[ÂÛÎÄͶ¸å] ¿ç³ö°æÉçÉÌͶ¸å 20+5 ÄߺÃ520 2026-07-13 5/250 2026-07-16 17:24 by ÎÞÊÓ¾üÍÅ
[»ù½ðÉêÇë] ÉúÉúÃüѧ²¿»áÆÀ +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by À²À²Ð¡ÄÝ×Ó
[»ù½ðÉêÇë] H¿Ú£¬ÃæÉÏ£¬Ê±¼ä´Á7ÔÂ7ÈÕ¡£ +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[»ù½ðÉêÇë] »áÆÀ +5 ÎÞÃûÕßµÇɽ 2026-07-15 5/250 2026-07-15 18:15 by ¸Ü¾ÍÊÇÄã¶Ô
[»ù½ðÉêÇë] E¿Ú»áÆÀ +11 ²¼²¼ºÍÒ»¶þ 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[»ù½ðÉêÇë] E¿ÚÇàC»áÆÀ½áÊøÁËÂ𣿠40+3 ССÊ鳿c 2026-07-14 9/450 2026-07-15 16:04 by °²¾²µÄÖªÁË
[ÂÛÎÄͶ¸å] ũҵ»úеѧ±¨Í¶¸åÖÜÆÚ +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[»ù½ðÉêÇë] ¹ú×ÔÈ»ÃæÉÏD¿ÚÆíµ» +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by °×Ë®Öó
[»ù½ðÉêÇë] ûÊÕµ½ÏûÏ¢£¬ÔÙ´ÎÅãÅÜ +3 Ìì¿ÕÐdz½fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[»ù½ðÉêÇë] Ã÷ÌìE¿ÚÃæÉÏ»áÆÀ +8 ²¼²¼ºÍÒ»¶þ 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û