Inclusion de matrices

Version imprimable

bonjour,j'ai un problème en fortran,voilà le code:
---------------------------------------------

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
	do I=1,120
	do J=1,120
      do i=1,30
      do j=1,30
	H6(i,j)=(0.d0)
c      write(*,*) 'H6(i,j)=',H6(i,j)        
	if(i.eq.j) then
	call sub1(A,i,y1,y2)
	H1(i,j)=-(y1+y2)
c	write(*,*) 'H1(i,j)=',H1(i,j)
      call sub5(A,i,y7)
      H4(i,j)=(y7)
c      write(*,*) 'H4(i,j)=',H4(i,j)
      call sub3(A,i,y5)
	H3(i,j)=(y5)
c      write(*,*) 'H3(i,j)=',H3(i,j)      
      call sub11(A,i,y3,y4)
	H2(i,j)=-(y3+y4)
c      write(*,*) 'H2(i,j)=',H2(i,j)      		
	else
	call sub2(A,i,y11,y22)
	H1(i,j)=-(y11+y22)
c	write(*,*) 'H1(i,j)=',H1(i,j)
      call sub4(A,i,y6)
	H3(i,j)=(y6)
c      write(*,*) 'H3(i,j)=',H3(i,j)
      call sub6(A,i,y8)
      H4(i,j)=(y8)
c      write(*,*) 'H4(i,j)=',H4(i,j)      
      call sub22(A,i,y33,y44)
	H2(i,j)=-(y33+y44)
c	write(*,*) 'H2(i,j)=',H2(i,j)      	
	endif
	enddo
	enddo
	Hbv(i,j)=
      enddo

---------------------------------------

je veux inclure N matrices de dimensions (30,30) dans une seule matrice de dimension (120,120).
svp je vous demande de l'aide.merci.

ta grosse matrice est-elle définie par blocs de 30x30?
auquel cas c'est facile:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
 
subroutine inclure(Y,H,I,J)
real(kind=8),dimension(120,120),intent(inout) :: H
real(kind=8),dimension(30,30),intent(in) :: Y
integer,intent(in) :: I,J
  do ii=1,30
    do jj=1,30  
      H(I+II-1,J+JJ-1)=Y(II,JJ)
    end do
  end do
return
end subroutine inclure

sinon tu as besoin d'une matrice d'inclusion de 30x30x2 donnant pour chaque élémant sa position dans la grande matrice:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
 
subroutine inclure2(Y,H,INC)
real(kind=8),dimension(120,120),intent(inout) :: H
real(kind=8),dimension(30,30),intent(in) :: Y
integer,dimension(30,30,2),intent(in) :: INC
  do ii=1,30
    do jj=1,30  
      H(INC(II,JJ,1),INC(II,JJ,2))=Y(II,JJ)
    end do
  end do
return
end subroutine inclure2

03/11/2010, 14h22
anahaithem

tu peux m'expliquer en détails ou bien d'une autre manière.merci

edit :

svp donner moi le code en algorithme.merci
04/11/2010, 09h09
genteur slayer
tu as deux solution la première:
Code:

1 2 3 4 5 6 7 *--------------------* | | | a---------* | | | A | | | *---------* | *--------------------*K
A est collée en bloc dans K
il suffit alors de connaitre les indices du premier élément de A (ici en a)
je les ai appelé I et J dans la première subroutine

sinon sdans le cas n°2, chaque élément de A a une position propre dans K auquel cas, on a besoin de connaitre la position de chaque élément le plus simple est de remplir une matrice de même taille que A contenant des couples d'indices.
on utilise alors la seconde subroutine que je te propose
05/11/2010, 12h43
amira ben ahmed

j'ai pas bien compris l'instruction suivante
H(I+II-1,J+JJ-1)=Y(II,JJ)
merci de m'expliquer plus encore car je l'ai utilisé mais pas de résultat ça dépasse la dimension 120*120
08/11/2010, 09h02
genteur slayer

H(I+II-1,J+JJ-1)=Y(II,JJ)

signifie que l'on colle le bloc Y dans H à la position (I,J) I et J doivent donc être plus petit que (dans ton cas) 120-30=90