²é¿´: 1610  |  »Ø¸´: 7
±¾Ìû²úÉú 1 ¸ö ³ÌÐòÇ¿Ìû £¬µã»÷ÕâÀï½øÐв鿴
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

kathy2008

ľ³æ (ÕýʽдÊÖ)

[ÇóÖú] ÈçºÎ´Ó¸ß˹Êä³öÎļþ¿ìËÙÌá³ö pai ¹ìµÀÐÅÏ¢¡£

ÈçÌâ¡£´Ó¸ß˹Êä³öÎļþÌá³öÁËeigenvectorÄÇÒ»²¿·Ö³öÀ´£¬¼´¸½¼þ1¡£ÏÖÔÚÐèÒªµÃµ½ pai ¹ìµÀÐÅÏ¢¡£¼´¸½¼þ2¡£ ¸½¼þ2 ¶ÔÓ¦ÓÚ¸½¼þ1µÄ32ºÅ£¬35ºÅ£¬38ºÅ£¬39ºÅ£¬40ºÅ£¬41ºÅ¹ìµÀ£¨Õ¼¾Ý¹ìµÀ£©£¬42ºÅÒ»Ö±µ½47ºÅ£¨·ÇÕ¼¾Ý¹ìµÀ£©µÄ2PxÖµ¡£ÇóһС³ÌÐò¡£ÇëÖ¸µã¡£Ð»Ð»¡£
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

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

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

snoopyzhao

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

¡ï
jjdg(½ð±Ò+1): ¸Ðл²ÎÓë 2011-07-04 12:46:36
ÎÒ²Å×¢Òâµ½£¬Äã°ÑÎÒÔÚ¶þÂ¥¸ø³öµÄ³ÌÐòÖеÄ
CODE:
read(line,fm) px(i,(j+1):(j+n))

¸Ä³ÉÁË
CODE:
read(line,fm) ppx(i,(j+1):(j+n))

º¦µÃÎÒŪÁ˰ëÌì²ÅÖªµÀΪɶ³öÀ´µÄ½á¹û×ÜÊDz»¶Ô¡£
7Â¥2011-07-04 09:48:55
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 8 ¸ö»Ø´ð

snoopyzhao

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

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

¡ï ¡ï
ben_ladeng(½ð±Ò+2): ºÜÏêϸ£¬´ýÂ¥Ö÷ÆÀ¶¨ºó½±Àø³ÌÐòÇ¿Ìû 2011-06-18 17:42:21
kathy2008(½ð±Ò+10): 2011-06-19 13:00:43
΢³¾¡¢ÃÎÏë(³ÌÐòÇ¿Ìû+1): 2011-06-19 17:04:28
´ó¸ÅÕâ¸öÑù×Ó¡£Ö»ÊÇÐèÒªÊÖ¹¤ÊäÈë¹ìµÀºÅ£¨ÕâÑù¿ÉÄÜÁé»îһЩ£©£¬Ã¿´ÎÊäÈëÒ»¸ö¹ìµÀÐòºÅ£¬»Ø³µ£¬ÊäÈë 0 Ôò½áÊøÕû¸ö³ÌÐò¡­¡­
CODE:
program ei
real, dimension(:,:), allocatable :: px,ppx
character(len=256) :: line
character(len=40) :: fm
integer :: nrow, ncol, i, j, k, ios

open(unit=12, file='eigenvector.out', status='old')
open(unit=13, file='2px.out', status='new')

do
   read(12,'(a)', iostat=ios) line
   if (ios /= 0) exit
   if (index(line,'EIGENVALUES') /= 0) then
      nrow=0
      ncol=0
      do
         read(12,'(a)', iostat=ios) line
         if (ios /= 0) exit
         if (line(1:4) == '    ') exit
         ncol=ncol+1
         if (index(line, '2PX') /= 0) nrow=nrow+1
      end do
      exit
   end if
end do

!write (*,*) nrow, ncol
rewind (12)

allocate(px(nrow,ncol),ppx(nrow,ncol))

i=0
j=0
do
   read(12,'(a)', iostat=ios) line
   if (ios /= 0) exit
   if (i == nrow) then
       i=0
       j=j+n
   end if
   if (index(line, '2PX') /= 0) then
       line = line(21:)
!      write(*,*) trim(line)
       i=i+1
       n = len_trim(line)/10
       write(fm,'(a,i0,a)') '(', n, 'f10.5)'
!      write(*,*) j
       read(line,fm) px(i,(j+1):(j+n))
   end if
end do

k=0
do
   write(*,*) 'please input a number between 1 and ', nrow, 'end the program by 0.'
   read(*,*) i
   if(i==0) exit
   k=k+1
   ppx(:,k) = px(:,i)
end do

!write(*,*) k/5, mod(k,5)

if (k>=5) then
   do j=1,k/5
      do i=1,nrow
         write(13,'(5f10.5)') ppx(i,(j-1)*5+1:j*5)
      end do
      write(13,*)
   end do
end if
if (mod(k,5) /=0) then
   write(fm,'(a,i0,a)') '(', mod(k,5), 'f10.5)'
   do i=1,nrow
      write(13, fm) ppx(i,(k/5*5+1):k)
   end do
end if

end program ei

2Â¥2011-06-18 16:41:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

kathy2008

ľ³æ (ÕýʽдÊÖ)

¡ï
dubo(½ð±Ò+1): »¶Ó­³£À´³ÌÐòÓïÑÔ°æÌÖÂÛ 2011-07-31 13:36:05
ÒýÓûØÌû:
Originally posted by snoopyzhao at 2011-06-18 16:41:08:
´ó¸ÅÕâ¸öÑù×Ó¡£Ö»ÊÇÐèÒªÊÖ¹¤ÊäÈë¹ìµÀºÅ£¨ÕâÑù¿ÉÄÜÁé»îһЩ£©£¬Ã¿´ÎÊäÈëÒ»¸ö¹ìµÀÐòºÅ£¬»Ø³µ£¬ÊäÈë 0 Ôò½áÊøÕû¸ö³ÌÐò¡­¡­

[code]
program ei
real, dimension(:,, allocatable :: px,ppx
character(len=256 ...

ÀûÓøóÌÐòÌáÈ¡¦Ð¹ìµÀÐÅÏ¢£¬±¨´í¡£ÐÅÏ¢ÈçÏÂ
At line 48 of file eigen-pai-nc3h7-r2.f
Fortran runtime error: Bad value during floating point read

googleҲûÓÐÕÒµ½Ó¦¶ÔÖ®²ß£¬Çë¸ßÊÖÖ¸µã¡£
3Â¥2011-07-03 15:06:57
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

snoopyzhao

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

¡ï
dubo(½ð±Ò+1): »¶Ó­³£À´³ÌÐòÓïÑÔ°æÌÖÂÛ 2011-07-31 13:36:26
ÄãµÄ .out ÎļþÊÇÕ¦Éú³ÉµÄ£¬Õâ´ÎµÄÕâ¸öÎļþ±ÈÉϴεÄÎļþÿһÐÐÇ°Ãæ¶àÁËÒ»¸ö¿Õ¸ñ¡­¡­

ËùÒÔ£¬Äã°Ñ³ÌÐòÖУº
CODE:
line = line(21:)

¸Ä³É
CODE:
line = line(22:)

¾Í¿ÉÒÔÁË¡­¡­
4Â¥2011-07-03 22:28:17
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 302Çóµ÷¼Á +11 ºôºôºô¡£¡£¡£¡£ 2026-03-17 11/550 2026-03-21 08:29 by JourneyLucky
[¿¼ÑÐ] Äϲý´óѧ²ÄÁÏר˶311·ÖÇóµ÷¼Á +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[¿¼ÑÐ] 316Çóµ÷¼Á +6 ÁºÜçö© 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666£¡
[¿¼ÑÐ] 332Çóµ÷¼Á +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÄÏʦ´ó 070300£¨»¯Ñ§£©304·ÖÇóµ÷¼Á +3 0703ÎäÜ·»ÛÑ©304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[¿¼ÑÐ] 304Çóµ÷¼Á +6 ÂüÊâ2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏר˶ӢһÊý¶þ306 +7 z1z2z3879 2026-03-18 7/350 2026-03-20 23:48 by JourneyLucky
[¿¼ÑÐ] 321Çóµ÷¼Á +9 ºÎÈó²É123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[¿¼ÑÐ] 330Çóµ÷¼Á +4 С²Ä»¯±¾¿Æ 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[¿¼ÑÐ] ±¾ÈË¿¼085602 »¯Ñ§¹¤³Ì ר˶ +19 ²»ÖªµÀ½Ðʲô£¡ 2026-03-15 21/1050 2026-03-20 20:48 by zhukairuo
[¿¼ÑÐ] 295²ÄÁÏÇóµ÷¼Á£¬Ò»Ö¾Ô¸Î人Àí¹¤085601ר˶ +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[¿¼ÑÐ] ²ÄÁϹ¤³Ìר˶µ÷¼Á +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[¿¼ÑÐ] 085601ר˶£¬×Ü·Ö342Çóµ÷¼Á£¬µØÇø²»ÏÞ +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[¿¼ÑÐ] ÉúÎïѧ071000 329·ÖÇóµ÷¼Á +3 ÎÒ°®ÉúÎïÉúÎﰮΠ2026-03-17 3/150 2026-03-18 10:12 by macy2011
[¿¼ÑÐ] 334Çóµ÷¼Á +3 Ö¾´æ¸ßÔ¶ÒâÔÚ»úÐ 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[¿¼ÑÐ] 293Çóµ÷¼Á +11 zjlµÄºÅ 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[¿¼ÑÐ] 283Çóµ÷¼Á +3 Ìý·ç¾ÍÊÇÓꣻ 2026-03-16 3/150 2026-03-17 07:41 by ÈÈÇéɳĮ
[¿¼ÑÐ] 070305Çóµ÷¼Á +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[¿¼ÑÐ] ±¾¿ÆÄϾ©´óѧһ־Ը´¨´óҩѧ327 +3 ÂóÌï¸ûÕß 2026-03-14 3/150 2026-03-14 20:04 by ÍâÐÇÎÄÃ÷
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û