²é¿´: 1471  |  »Ø¸´: 4
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

physliebe

гæ (СÓÐÃûÆø)

[ÇóÖú] moduleÎÞ·¨±àÒ룿 ÒÑÓÐ1È˲ÎÓë

Óøß˹ÏûÔª·¨Çó½âÏßÐÔ·½³Ì×é
CODE:
module gauss
    contains
    subroutine solve(A,b,x,N)
        implicit real*8(a-z)
        integer::i,k,N
        real*8::A(N,N),b(N),x(N)
        real*8::Aup(N,N),bup(N)
        real*8::Ab(N,N+1)
        Ab(1:N,1:N)=A
        Ab(:,N+1)=b
        !+++++++++++++++++++++++++++++++++++++++++
        £¡¸ß˹ÏûÈ¥·¨µÄºËÐIJ¿·Ö
        do k=1,N-1
            do i=k+1,N
                temp=Ab(i,k)/Ab(k,k)
                Ab(i,:)=Ab(i,:)-temp*Ab(k,:)
            end do
        end do
        Aup(:,:)=Ab(1:N,1:N)
        bup(:)=Ab(:,N+1)
        !+++++++++++++++++++++++++++++++++++++++++
        £¡Çó½âÉÏÈý½Ç¾ØÕóÐ͵ÄÏÈÐÔ·½³Ì×é
        call uptri(Aup,bup,x,N)   
        end subroutine solve
     subroutine uptri(A,b,x,N)        
        implicit real*8(a-z)         
        integer::i,j,N
        real*8::A(N,N),b(N),x(N)      
        x(N)=b(N)/A(N,N)              
        !»Ø´ú²¿·Ö
        do i=N-1,1,-1               
            x(i)=b(i)
            do j=i+1,N               
                x(i)=x(i)-a(i,j)*x(j)
            end do
            x(i)=x(i)/A(i,i)
        end do
    end subroutine uptri
end module gauss
program main
    use gauss
    implicit real*(a-z)
    integer,parameter::N=4
    integer::i,j
    real*8::A(N,N),b(N),x(N)
    open(unit=11,file='fin.txt')
    open(unit=12,file='fout.txt')
    read(11,*) ((A(i,j),j=1,N),i=1,N),b
    call solve(A,b,x,N)   !xÔÚÕâÀïÃæµÄ×÷ÓÃ
    write(12,101) x
    101 format(T5,'¸ß˹ÏûÈ¥·¨¼ÆËã½á¹û',/,T4,4(f12.8))
end program main

´íÎóÐÅÏ¢£º
end module gauss
                1
Fatal Error: Can't open module file 'D:\fortran_programming\gauss_1\main.os_output_dir/gauss.mod0' for writing at (1): No such file or directory
Process terminated with status 1 (0 minute(s), 0 second(s))
1 error(s), 0 warning(s) (0 minute(s), 0 second(s))

[ Last edited by jjdg on 2014-2-9 at 10:01 ]
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

Whenever you feel like criticizing any one,just remember that all the people in this world haven¡¯t had the advantages that you¡¯ve had.
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

physliebe

гæ (СÓÐÃûÆø)

ÒýÓûØÌû:
3Â¥: Originally posted by snoopyzhao at 2014-02-09 05:22:20
ÖÁÉÙÓà gfortran ±àÒëÄãÌù³öÀ´µÄ´úÂë²»Ö¹ÄãÉÏÃæÌáµ½µÄ´íÎ󣬵«°´ÕÕ³ö´íÐÅϢȥÐÞ¸ÄÓÖ²»»áÓÐʲôÎÊÌâ

1£©ÄãÔÚ³ÌÐòÖÐʹÓÃÖÐÎĵĸÐ̾ºÅÈ¥×¢ÊÍ´úÂëÊDz»ÐеÄ
2£©Ö÷³ÌÐòÖеÄÒþʽÉùÃ÷²»È«£¬ËùÒÔ´úÂëÊÇÓÐÎÊÌâµÄ

Ïû³ýÁË ...

лл£¬Ö÷ÒªÊÇÎÒÓÃÁËcodeblocks13.12×÷Ϊ¼¯³É¿ª·¢»·¾³£¬Ã»ÓÐÅäÖúÃToolchain excutables¡£±ê׼д·¨ÊDz»ÊÇreal,dimension(8)ÄØ£¿
Whenever you feel like criticizing any one,just remember that all the people in this world haven¡¯t had the advantages that you¡¯ve had.
4Â¥2014-02-09 12:28:40
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 5 ¸ö»Ø´ð

physliebe

гæ (СÓÐÃûÆø)

΢Ц±íÇ鴦ΪðºÅ¡®£º¡¯
Whenever you feel like criticizing any one,just remember that all the people in this world haven¡¯t had the advantages that you¡¯ve had.
2Â¥2014-02-08 20:56:04
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

snoopyzhao

ÖÁ×ðľ³æ (Ö°Òµ×÷¼Ò)

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

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
jjdg: ½ð±Ò+1, ¸Ðл²ÎÓë 2014-02-09 10:02:02
physliebe: ½ð±Ò+5 2014-02-09 12:21:04
ÖÁÉÙÓà gfortran ±àÒëÄãÌù³öÀ´µÄ´úÂë²»Ö¹ÄãÉÏÃæÌáµ½µÄ´íÎ󣬵«°´ÕÕ³ö´íÐÅϢȥÐÞ¸ÄÓÖ²»»áÓÐʲôÎÊÌâ

1£©ÄãÔÚ³ÌÐòÖÐʹÓÃÖÐÎĵĸÐ̾ºÅÈ¥×¢ÊÍ´úÂëÊDz»ÐеÄ
2£©Ö÷³ÌÐòÖеÄÒþʽÉùÃ÷²»È«£¬ËùÒÔ´úÂëÊÇÓÐÎÊÌâµÄ

Ïû³ýÁËÉÏÃæµÄÁ½¸öÎÊÌâºÃÏó¾Í¿ÉÒÔÁË

²»ÔÞ³ÉʹÓÃÒþʽÉùÃ÷£¬²»ÔÞ³ÉÓà *8 ÕâÑùµÄ·Ç±ê׼д·¨¡­¡­
3Â¥2014-02-09 05:22:20
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

snoopyzhao

ÖÁ×ðľ³æ (Ö°Òµ×÷¼Ò)

ÒýÓûØÌû:
4Â¥: Originally posted by physliebe at 2014-02-09 12:28:40
лл£¬Ö÷ÒªÊÇÎÒÓÃÁËcodeblocks13.12×÷Ϊ¼¯³É¿ª·¢»·¾³£¬Ã»ÓÐÅäÖúÃToolchain excutables¡£±ê׼д·¨ÊDz»ÊÇreal,dimension(8)ÄØ£¿...

real*8 a µÄ±ê׼д·¨ÊÇ real(kind=kind(1.0d0)) :: a»òÕ߸ü±ê×¼µÄд·¨ÊÇ£º
integer, parameter :: dp = SELECTED_REAL_KIND(15,307)
real(kind=dp) :: a
5Â¥2014-02-09 13:14:01
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 085700×ÊÔ´Óë»·¾³308Çóµ÷¼Á +3 īīĮ 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[¿¼ÑÐ] 274Çóµ÷¼Á +5 S.H1 2026-03-18 5/250 2026-03-18 21:27 by guosr9609
[¿¼ÑÐ] ²ÄÁÏרҵÇóµ÷¼Á +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 314Çóµ÷¼Á +8 ÎÞи¿É»÷µÄ¾ÞÈË 2026-03-12 8/400 2026-03-18 14:50 by haxia
[¿¼ÑÐ] 08¹¤¿Æ 320×Ü·Ö Çóµ÷¼Á +5 À滨çóÍí·ç 2026-03-17 5/250 2026-03-18 14:49 by haxia
[¿¼ÑÐ] 297Çóµ÷¼Á +8 Ï·¾«µ¤µ¤µ¤ 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[¿¼ÑÐ] 302Çóµ÷¼Á +10 ºôºôºô¡£¡£¡£¡£ 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[¿¼ÑÐ] ¹¤¿Æ²ÄÁÏ085601 279Çóµ÷¼Á +6 À§ÓÚÐdz¿ 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[¿¼ÑÐ] 268Çóµ÷¼Á +6 ¼òµ¥µã0 2026-03-17 6/300 2026-03-18 09:04 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 334Çóµ÷¼Á +3 Ö¾´æ¸ßÔ¶ÒâÔÚ»úÐ 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +6 Ϊѧ666 2026-03-16 6/300 2026-03-17 20:15 by peike
[¿¼ÑÐ] 326Çóµ÷¼Á +4 ŵ±´¶û»¯Ñ§½±êéê 2026-03-15 7/350 2026-03-16 17:11 by ŵ±´¶û»¯Ñ§½±êéê
[¿¼ÑÐ] 0703Ò»Ö¾Ô¸211 285·ÖÇóµ÷¼Á +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by ŶŶ123
[¿¼ÑÐ] 080500£¬²ÄÁÏѧ˶302·ÖÇóµ÷¼ÁѧУ +4 ³õʶ¿ÉÀÖ 2026-03-14 5/250 2026-03-14 21:08 by peike
[¿¼ÑÐ] ±¾¿ÆÄϾ©´óѧһ־Ը´¨´óҩѧ327 +3 ÂóÌï¸ûÕß 2026-03-14 3/150 2026-03-14 20:04 by ÍâÐÇÎÄÃ÷
[¿¼ÑÐ] Öпƴó²ÄÁÏר˶319Çóµ÷¼Á +3 ÃÏöβÄÁÏ 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[¿¼ÑÐ] 255Çóµ÷¼Á +3 Àî¼Î»Û£¬ 2026-03-12 4/200 2026-03-14 16:58 by ÓÐÖ»ÀêÅ«
[¿¼ÑÐ] ¸´ÊÔµ÷¼Á +3 ºôºô£¿~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[¿¼ÑÐ] 308 085701 ËÄÁù¼¶ÒѹýÇóµ÷¼Á +7 ÎÂÇÇÇÇÇÇÇÇ 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[¿¼ÑÐ] ¹¤¿Æ278·ÖÇóµ÷¼Á +5 ÖÜÂýÈȰ¡ 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û