1 : C ASSMAT SOURCE BP208322 12/10/03 21:15:00 7517 2 : SUBROUTINE ASSMAT(IPT1,X1,RI1,X2,RI2,X3,RI3,X4,RI4,XMATRI) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : * 6 : ********************************************************************************* 7 : * CONSTITUTION D'UNE MATRICE 2Nx2N POUR LE QZ * 8 : * A PARTIR DE MATRICES PROJETEES NxN * 9 : * _____________________________________________________________________________ * 10 : * * 11 : * DATE : le 24 Fevrier 1995 * 12 : * AUTEUR : Nicolas BENECH * 13 : * _____________________________________________________________________________ * 14 : * * 15 : * MODULE(S) APPELANT(S) : QZCONS * 16 : * * 17 : * MODULE(S) APPELE (S): COPMAT * 18 : * _____________________________________________________________________________ * 19 : * * 20 : * EN ENTREE : * 21 : * IPT1 : maillage permettant d'ordonner les valeurs * 22 : * (Xi, Ri) : coefficient et matrice a stoker a l'angle i * 23 : * (i=1, angle superieur gauche, i=2 angle superieur droit) * 24 : * (i=3, angle inferieur gauche, i=4 angle inferieur droit) * 25 : * * 26 : * EN SORTIE : * 27 : * XMATRI : la matrice constituee * 28 : * (actif et modifiable en entree/sortie) * 29 : * * 30 : ********************************************************************************* 31 : * 32 : -INC CCOPTIO 33 : -INC SMRIGID 34 : * 35 : LOGICAL AFFICH 36 : REAL*8 X1,X2,X3,X4 37 : INTEGER N1, N2 38 : * 39 : * Affichage des messages pour verification 40 : * 41 : AFFICH = .FALSE. 42 : * 43 : c SEGACT, XMATRI 44 : N1= INT (XMATRI.RE(/1)/2) 45 : N2= INT (XMATRI.RE(/2)/2) 46 : c SEGDES, XMATRI 47 : * 48 : IF (AFFICH) WRITE (*,*) 'Coin superieur gauche ...' 48 : CALL COPMAT(IPT1,X1,RI1,0,0,XMATRI) 49 : IF (AFFICH) WRITE (*,*) 'Coin superieur droit ...' 49 : CALL COPMAT(IPT1,X2,RI2,0,N2,XMATRI) 50 : IF (AFFICH) WRITE (*,*) 'Coin inferieur gauche ...' 50 : CALL COPMAT(IPT1,X3,RI3,N1,0,XMATRI) 51 : IF (AFFICH) WRITE (*,*) 'Coin inferieur droit ...' 51 : CALL COPMAT(IPT1,X4,RI4,N1,N2,XMATRI) 52 : * 53 : RETURN 54 : END 55 : 56 : 57 : 58 : 59 : 60 : 61 : 62 : 63 : 64 :
© Cast3M 2003 - All rights reserved.
Disclaimer