Télécharger addite.eso

Retour à la liste

Numérotation des lignes :

  1. C ADDITE SOURCE PV 20/04/01 21:15:08 10569
  2.  
  3. C Sous-programme elementaire pour effectuer la translation, la rotation
  4. C ou l'affinite d'un objet
  5. C 10/2003 : Prise en compte cas IDIM=1.
  6.  
  7. SUBROUTINE ADDITE(X,IPT1,IPT2,ICPR,IARG)
  8.  
  9. IMPLICIT INTEGER(I-N)
  10. IMPLICIT REAL*8 (A-H,O-Z)
  11.  
  12.  
  13. -INC PPARAM
  14. -INC CCOPTIO
  15. -INC SMELEME
  16. -INC SMCOORD
  17.  
  18. DIMENSION X(*)
  19. SEGMENT ICPR(nbpts)
  20.  
  21. COMMON / CTOURN / XPT1,YPT1,ZPT1,XV1,YV1,ZV1,XV2,YV2,ZV2,
  22. . XVEC,YVEC,ZVEC,ANGLE,ICLE
  23.  
  24. idimp1=IDIM+1
  25. XVECS=XVEC
  26. YVECS=YVEC
  27. ZVECS=ZVEC
  28. COSA=COS(ANGLE)
  29. SINA=SIN(ANGLE)
  30. ANG2=ANGLE*ANGLE
  31.  
  32. C Initialisation du maillage IPT2 transforme de IPT1
  33. NBSOUS=0
  34. IF (IARG.EQ.0) THEN
  35. NBREF=0
  36. ELSE
  37. NBREF=IPT1.LISREF(/1)
  38. ENDIF
  39. NBNN=IPT1.NUM(/1)
  40. NBELEM=IPT1.NUM(/2)
  41. SEGINI IPT2
  42. IPT2.ITYPEL=IPT1.ITYPEL
  43. DO i=1,NBELEM
  44. IPT2.ICOLOR(i)=IPT1.ICOLOR(i)
  45. ENDDO
  46.  
  47. C ICLE / KCLE : indique la transformation a realiser
  48. C KCLE = 1 : translation (operateurs PLUS,MOIN,DEDU 'TRAN')
  49. C KCLE = 2 : rotation (operateurs TOUR,DEDU 'ROTA')
  50. C KCLE = 3 : homothetie (operateur HOMO)
  51. C KCLE = 4 : affinite (operateur AFFI)
  52. C KCLE = 5 : symetrie point (operateur SYME 'POINT')
  53. C KCLE = 6 : symetrie droite (operateur SYME 'DROIT')
  54. C KCLE = 7 : symetrie plan (operateur SYME 'PLAN')
  55. C KCLE = 8 : projection plan (operateur PROJ 'PLAN')
  56. C KCLE = 9 : projection sphere (operateur PROJ 'SPHE')
  57. C KCLE = 10 : projection cylindre (operateur PROJ 'CYLI')
  58. C KCLE = 11 : projection conique (operateur PROJ 'CONI')
  59. C KCLE = 12 : projection torique (operateur PROJ 'TORI')
  60. C KCLE = 13 : projection droite (operateur PROJ 'DROI')
  61. C KCLE = 14 : projection cercle (operateur PROJ 'CERC')
  62. KCLE=ICLE
  63. IF (ICLE.GE.10000) KCLE=ICLE-10000
  64.  
  65. C Reservation de la place dans XCOOR
  66. segdes mcoord
  67. SEGACT MCOORD*mod
  68. NPACT=NBPTS
  69. NXCOUR=NPACT*idimp1
  70. NBPTS=NPACT+NBELEM*NBNN
  71. SEGADJ MCOORD
  72.  
  73. C Boucle sur les noeuds du maillage
  74. C Creation du noeud de IPT2 image de IPT1.NUM(i,j) et stocke dans ICPR
  75. DO 10 j=1,NBELEM
  76. DO 10 i=1,NBNN
  77. C L'image de IPT1.NUM(i,j) a-t-elle deja ete creee ?
  78. IF (ICPR(IPT1.NUM(i,j)).EQ.0) THEN
  79. IREF=(IPT1.NUM(i,j)-1)*idimp1
  80. GOTO (51,52,53,54,55,56,57,58,59,60,61,62,63,64),KCLE
  81. C Translation de vecteur X :
  82. 51 DO k=1,IDIM
  83. XCOOR(NXCOUR+k)=X(k)+XCOOR(IREF+k)
  84. ENDDO
  85. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  86. GOTO 70
  87. C Rotation (2D/3D) :
  88. 52 XD=XCOOR(IREF+1)-XPT1
  89. YD=XCOOR(IREF+2)-YPT1
  90. ZD=0.
  91. IF (IDIM.GE.3) ZD=XCOOR(IREF+3)-ZPT1
  92. XE=XD*XV1+YD*YV1+ZD*ZV1
  93. YE=XD*XV2+YD*YV2+ZD*ZV2
  94. ZE=XD*XVEC+YD*YVEC+ZD*ZVEC
  95. XD=XE*COSA-YE*SINA
  96. YD=XE*SINA+YE*COSA
  97. ZD=ZE
  98. XCOOR(NXCOUR+1)=XD*XV1+YD*XV2+ZD*XVEC+XPT1
  99. XCOOR(NXCOUR+2)=XD*YV1+YD*YV2+ZD*YVEC+YPT1
  100. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=XD*ZV1+YD*ZV2+ZD*ZVEC+ZPT1
  101. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  102. GOTO 70
  103. C Homothetie :
  104. 53 XCOOR(NXCOUR+1)=(XCOOR(IREF+1)-XPT1)*ANGLE+XPT1
  105. IF (IDIM.GE.2) XCOOR(NXCOUR+2)=(XCOOR(IREF+2)-YPT1)*ANGLE+YPT1
  106. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=(XCOOR(IREF+3)-ZPT1)*ANGLE+ZPT1
  107. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)*ANGLE
  108. GOTO 70
  109. C Affinite (2D/3D) :
  110. 54 XD=XCOOR(IREF+1)-XPT1
  111. YD=XCOOR(IREF+2)-YPT1
  112. ZD=0.
  113. IF (IDIM.GE.3) ZD=XCOOR(IREF+3)-ZPT1
  114. XE=XD*XV1+YD*YV1+ZD*ZV1
  115. YE=XD*XV2+YD*YV2+ZD*ZV2
  116. ZE=(XD*XVEC+YD*YVEC+ZD*ZVEC)*ANGLE
  117. XCOOR(NXCOUR+1)=XE*XV1+YE*XV2+ZE*XVEC+XPT1
  118. XCOOR(NXCOUR+2)=XE*YV1+YE*YV2+ZE*YVEC+YPT1
  119. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=XE*ZV1+YE*ZV2+ZE*ZVEC+ZPT1
  120. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  121. GOTO 70
  122. C Symetrie par rapport a un point :
  123. 55 XCOOR(NXCOUR+1)=2*XPT1-XCOOR(IREF+1)
  124. IF (IDIM.GE.2) XCOOR(NXCOUR+2)=2*YPT1-XCOOR(IREF+2)
  125. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=2*ZPT1-XCOOR(IREF+3)
  126. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  127. GOTO 70
  128. C Symetrie par rapport a une droite (2D/3D) :
  129. 56 XD=XCOOR(IREF+1)-XPT1
  130. YD=XCOOR(IREF+2)-YPT1
  131. ZD=0.
  132. IF (IDIM.GE.3) ZD=XCOOR(IREF+3)-ZPT1
  133. PVEC=2*(XD*XVEC+YD*YVEC+ZD*ZVEC)
  134. XCOOR(NXCOUR+1)=XPT1+XVEC*PVEC-XD
  135. XCOOR(NXCOUR+2)=YPT1+YVEC*PVEC-YD
  136. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=ZPT1+ZVEC*PVEC-ZD
  137. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  138. GOTO 70
  139. C Symetrie par rapport a un plan (3D) :
  140. 57 XD=XCOOR(IREF+1)-XPT1
  141. YD=XCOOR(IREF+2)-YPT1
  142. ZD=XCOOR(IREF+3)-ZPT1
  143. PVEC=2*(XD*XVEC+YD*YVEC+ZD*ZVEC)
  144. XCOOR(NXCOUR+1)=XPT1+XD-XVEC*PVEC
  145. XCOOR(NXCOUR+2)=YPT1+YD-YVEC*PVEC
  146. XCOOR(NXCOUR+3)=ZPT1+ZD-ZVEC*PVEC
  147. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  148. GOTO 70
  149. C Projection sur un plan :
  150. 58 XPOIN=XCOOR(IREF+1)
  151. YPOIN=XCOOR(IREF+2)
  152. ZPOIN=XCOOR(IREF+3)
  153. IF (ICLE.GE.10000) THEN
  154. XVECS=XVEC-XPOIN
  155. YVECS=YVEC-YPOIN
  156. ZVECS=ZVEC-ZPOIN
  157. ENDIF
  158. SVECS=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  159. IF (SVECS.EQ.0.) CALL ERREUR(21)
  160. IF (IERR.NE.0) RETURN
  161. SVECS=SQRT(SVECS)
  162. XVECS=XVECS/SVECS
  163. YVECS=YVECS/SVECS
  164. ZVECS=ZVECS/SVECS
  165. DENUM=(XPOIN-XPT1)*XV1+(YPOIN-YPT1)*YV1+(ZPOIN-ZPT1)*ZV1
  166. DENOM=XVECS*XV1+YVECS*YV1+ZVECS*ZV1
  167. IF (DENOM.EQ.0.) CALL ERREUR(21)
  168. IF (IERR.NE.0) RETURN
  169. RAP=-DENUM/DENOM
  170. XCOOR(NXCOUR+1)=XPOIN+RAP*XVECS
  171. XCOOR(NXCOUR+2)=YPOIN+RAP*YVECS
  172. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=ZPOIN+RAP*ZVECS
  173. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  174. GOTO 70
  175. C Projection sur une sphere (3D) :
  176. 59 XPOIN=XCOOR(IREF+1)
  177. YPOIN=XCOOR(IREF+2)
  178. ZPOIN=XCOOR(IREF+3)
  179. IF (ICLE.GE.10000) THEN
  180. XVECS=XVEC-XPOIN
  181. YVECS=YVEC-YPOIN
  182. ZVECS=ZVEC-ZPOIN
  183. ENDIF
  184. SVECS=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  185. IF (SVECS.EQ.0.) CALL ERREUR(21)
  186. IF (IERR.NE.0) RETURN
  187. SVECS=SQRT(SVECS)
  188. XVECS=XVECS/SVECS
  189. YVECS=YVECS/SVECS
  190. ZVECS=ZVECS/SVECS
  191. SCA=XVECS*(XPT1-XPOIN)+YVECS*(YPT1-YPOIN)+ZVECS*(ZPT1-ZPOIN)
  192. XV=XVECS*SCA
  193. YV=YVECS*SCA
  194. ZV=ZVECS*SCA
  195. S2=(XPOIN+XV-XPT1)**2+(YPOIN+YV-YPT1)**2+(ZPOIN+ZV-ZPT1)**2
  196. IF (S2.GT.ANG2) CALL ERREUR(21)
  197. IF (IERR.NE.0) RETURN
  198. C=SQRT(ANG2-S2)
  199. IF (SCA.LT.0.) C=-C
  200. XCOOR(NXCOUR+1)=XPOIN+XV-C*XVECS
  201. XCOOR(NXCOUR+2)=YPOIN+YV-C*YVECS
  202. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=ZPOIN+ZV-C*ZVECS
  203. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  204. GOTO 70
  205. C Projection sur un cylindre (3D) :
  206. 60 XPOIN=XCOOR(IREF+1)
  207. YPOIN=XCOOR(IREF+2)
  208. ZPOIN=XCOOR(IREF+3)
  209. IF (ICLE.GE.10000) THEN
  210. XVECS=XVEC-XPOIN
  211. YVECS=YVEC-YPOIN
  212. ZVECS=ZVEC-ZPOIN
  213. ENDIF
  214. SVECS=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  215. IF (SVECS.EQ.0.) CALL ERREUR(21)
  216. IF (IERR.NE.0) RETURN
  217. SVECS=SQRT(SVECS)
  218. XVECS=XVECS/SVECS
  219. YVECS=YVECS/SVECS
  220. ZVECS=ZVECS/SVECS
  221. XV2=YVECS*ZV1-ZVECS*YV1
  222. YV2=ZVECS*XV1-XVECS*ZV1
  223. ZV2=XVECS*YV1-YVECS*XV1
  224. S2V2=XV2*XV2+YV2*YV2+ZV2*ZV2
  225. IF (S2V2.EQ.0.) CALL ERREUR(21)
  226. IF (IERR.NE.0) RETURN
  227. C2=(XVECS*XV1+YVECS*YV1+ZVECS*ZV1)**2
  228. IF (C2.EQ.1.) CALL ERREUR(21)
  229. IF (IERR.NE.0) RETURN
  230. XV3=XPT1-XPOIN
  231. YV3=YPT1-YPOIN
  232. ZV3=ZPT1-ZPOIN
  233. XV4=YV3*ZV1-ZV3*YV1
  234. YV4=ZV3*XV1-XV3*ZV1
  235. ZV4=XV3*YV1-YV3*XV1
  236. DNUM=(XV4*XVECS+YV4*YVECS+ZV4*ZVECS)**2
  237. DIS2=DNUM/S2V2
  238. IF (IIMPI.EQ.1) WRITE (IOIMP,1000) DIS2,ANGLE
  239. 1000 FORMAT (' DISTANCE AU CARRE DES AXES ',G12.5,
  240. . 'RAYON DU CYLINDRE ',G12.5)
  241. IF (DIS2.GT.ANG2) CALL ERREUR(21)
  242. IF (IERR.NE.0) RETURN
  243. DMU=SQRT((ANG2-DIS2)/(1.-C2))
  244. DNUM=XV2*XV4+YV2*YV4+ZV2*ZV4
  245. DLA=DNUM/S2V2
  246. DMU=SIGN(DMU,DLA)
  247. IF (IIMPI.NE.0) WRITE (IOIMP,1001) DLA,DMU
  248. 1001 FORMAT(' LAMBDA,MU ',2G13.5)
  249. XCOOR(NXCOUR+1)=XPOIN+XVECS*(DLA-DMU)
  250. XCOOR(NXCOUR+2)=YPOIN+YVECS*(DLA-DMU)
  251. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=ZPOIN+ZVECS*(DLA-DMU)
  252. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  253. GOTO 70
  254. C Projection sur un cone (3D) :
  255. 61 XPOIN=XCOOR(IREF+1)
  256. YPOIN=XCOOR(IREF+2)
  257. ZPOIN=XCOOR(IREF+3)
  258. IF (ICLE.GE.10000) THEN
  259. XVECS=XVEC-XPOIN
  260. YVECS=YVEC-YPOIN
  261. ZVECS=ZVEC-ZPOIN
  262. ENDIF
  263. SVECS=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  264. IF (SVECS.EQ.0.) CALL ERREUR(21)
  265. IF (IERR.NE.0) RETURN
  266. SVECS=SQRT(SVECS)
  267. XVECS=XVECS/SVECS
  268. YVECS=YVECS/SVECS
  269. ZVECS=ZVECS/SVECS
  270. XV2=XPOIN-XPT1
  271. YV2=YPOIN-YPT1
  272. ZV2=ZPOIN-ZPT1
  273. VECV1=XVECS*XV1+YVECS*YV1+ZVECS*ZV1
  274. VEC2=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  275. V2V1=XV2*XV1+YV2*YV1+ZV2*ZV1
  276. VECV2=XVECS*XV2+YVECS*YV2+ZVECS*ZV2
  277. V22=XV2*XV2+YV2*YV2+ZV2*ZV2
  278. A=VECV1*VECV1-ANGLE*VEC2
  279. B=2*(V2V1*VECV1-ANGLE*VECV2)
  280. C=V2V1*V2V1-ANGLE*V22
  281. DELTA=B*B-4*A*C
  282. IF (DELTA.LT.0.) CALL ERREUR(21)
  283. IF (IERR.NE.0) RETURN
  284. DEL=SQRT(DELTA)
  285. X1=(-B+DEL)/(2*A)
  286. XX=(-B-DEL)/(2*A)
  287. IF (ABS(X1).LT.ABS(XX)) XX=X1
  288. XCOOR(NXCOUR+1)=XPOIN+XX*XVECS
  289. XCOOR(NXCOUR+2)=YPOIN+XX*YVECS
  290. XCOOR(NXCOUR+3)=ZPOIN+XX*ZVECS
  291. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  292. GOTO 70
  293. C Projection sur un tore (3D) :
  294. 62 XPOIN=XCOOR(IREF+1)
  295. YPOIN=XCOOR(IREF+2)
  296. ZPOIN=XCOOR(IREF+3)
  297. IF (ICLE.GE.10000) THEN
  298. XVECS=XVEC-XPOIN
  299. YVECS=YVEC-YPOIN
  300. ZVECS=ZVEC-ZPOIN
  301. ENDIF
  302. SVECS=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  303. IF (SVECS.EQ.0.) CALL ERREUR(21)
  304. IF (IERR.NE.0) RETURN
  305. SVECS=SQRT(SVECS)
  306. XVECS=XVECS/SVECS
  307. YVECS=YVECS/SVECS
  308. ZVECS=ZVECS/SVECS
  309. PR2=XV2
  310. GR2=YV2
  311. XOP=XPOIN-XPT1
  312. YOP=YPOIN-YPT1
  313. ZOP=ZPOIN-ZPT1
  314. OPV=XOP*XVECS+YOP*YVECS+ZOP*ZVECS
  315. R2=XOP*XOP+YOP*YOP+ZOP*ZOP-GR2-PR2
  316. VA=XVECS*XV1+YVECS*YV1+ZVECS*ZV1
  317. QR2VA2=4*GR2*VA*VA
  318. OPA=XOP*XV1+YOP*YV1+ZOP*ZV1
  319. HR2PV=8*GR2*OPA*VA
  320. R=4*GR2*OPA*OPA-4*PR2*GR2
  321. RLMD=0
  322. C Resolution iterative de l'equation du 4eme degre
  323. DO ITER=1,25
  324. EXP1=RLMD*(RLMD+2*OPV)+R2
  325. FDLM=EXP1*EXP1+QR2VA2*RLMD*RLMD+HR2PV*RLMD+R
  326. FPDLM=4*EXP1*(RLMD+OPV)+QR2VA2*2*RLMD+HR2PV
  327. IF (FPDLM.EQ.0.) CALL ERREUR(40)
  328. IF (IERR.NE.0) RETURN
  329. CORR=FDLM/FPDLM
  330. IF (IIMPI.NE.0) WRITE(IOIMP,1016) ITER,RLMD,CORR
  331. 1016 FORMAT(' ITER ',I6,' LAMBDA ',G12.5,' CORR ',G12.5)
  332. RLMD=RLMD-CORR
  333. IF (RLMD.NE.0.) THEN
  334. IF (ABS(CORR/RLMD).LT.1E-5) GOTO 66
  335. ENDIF
  336. ENDDO
  337. CALL ERREUR(40)
  338. RETURN
  339. 66 XCOOR(NXCOUR+1)=XPOIN+RLMD*XVECS
  340. XCOOR(NXCOUR+2)=YPOIN+RLMD*YVECS
  341. XCOOR(NXCOUR+3)=ZPOIN+RLMD*ZVECS
  342. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  343. GOTO 70
  344. C Projection sur une droite (2D/3D) :
  345. 63 XPOIN=XCOOR(IREF+1)
  346. YPOIN=XCOOR(IREF+2)
  347. ZPOIN=0.
  348. IF (ICLE.GE.10000) THEN
  349. XVECS=XVEC-XPOIN
  350. YVECS=YVEC-YPOIN
  351. ZVECS=0.
  352. ENDIF
  353. SVECS=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  354. IF (SVECS.EQ.0.) CALL ERREUR(21)
  355. IF (IERR.NE.0) RETURN
  356. SVECS=SQRT(SVECS)
  357. XVECS=XVECS/SVECS
  358. YVECS=YVECS/SVECS
  359. ZVECS=ZVECS/SVECS
  360. DENUM=(XPOIN-XPT1)*XV1+(YPOIN-YPT1)*YV1
  361. DENOM=XVECS*XV1+YVECS*YV1+ZVECS*ZV1
  362. IF (DENOM.EQ.0.) CALL ERREUR(21)
  363. IF (IERR.NE.0) RETURN
  364. RAP=-DENUM/DENOM
  365. XCOOR(NXCOUR+1)=XPOIN+RAP*XVECS
  366. XCOOR(NXCOUR+2)=YPOIN+RAP*YVECS
  367. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=ZPOIN+RAP*ZVECS
  368. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  369. GOTO 70
  370. C Projection sur un cercle (2D/3D) :
  371. 64 XPOIN=XCOOR(IREF+1)
  372. YPOIN=XCOOR(IREF+2)
  373. ZPOIN=0.
  374. IF (ICLE.GE.10000) THEN
  375. XVECS=XVEC-XPOIN
  376. YVECS=YVEC-YPOIN
  377. ZVECS=0.
  378. ENDIF
  379. SVECS=XVECS*XVECS+YVECS*YVECS+ZVECS*ZVECS
  380. IF (SVECS.EQ.0.) CALL ERREUR(21)
  381. IF (IERR.NE.0) RETURN
  382. SVECS=SQRT(SVECS)
  383. XVECS=XVECS/SVECS
  384. YVECS=YVECS/SVECS
  385. ZVECS=ZVECS/SVECS
  386. SCA=XVECS*(XPT1-XPOIN)+YVECS*(YPT1-YPOIN)
  387. XV=XVECS*SCA
  388. YV=YVECS*SCA
  389. ZV=ZVECS*SCA
  390. S2=(XPOIN+XV-XPT1)**2+(YPOIN+YV-YPT1)**2+(ZPOIN+ZV-ZPT1)**2
  391. IF (S2.GT.ANG2) CALL ERREUR(21)
  392. IF (IERR.NE.0) RETURN
  393. C=SQRT(ANG2-S2)
  394. IF (SCA.LT.0.) C=-C
  395. XCOOR(NXCOUR+1)=XPOIN+XV-C*XVECS
  396. XCOOR(NXCOUR+2)=YPOIN+YV-C*YVECS
  397. IF (IDIM.GE.3) XCOOR(NXCOUR+3)=ZPOIN+ZV-C*ZVECS
  398. XCOOR(NXCOUR+idimp1)=XCOOR(IREF+idimp1)
  399. GOTO 70
  400. C Stockage du point image cree dans XCOOR et dans ICPR
  401. C On remplit le maillage IPT2
  402. 70 NXCOUR=NXCOUR+idimp1
  403. IPT2.NUM(i,j)=NXCOUR/idimp1
  404. ICPR(IPT1.NUM(i,j))=IPT2.NUM(i,j)
  405. ELSE
  406. C On recupere le noeud image dans ICPR et on remplit IPT2
  407. IPT2.NUM(i,j)=ICPR(IPT1.NUM(i,j))
  408. ENDIF
  409. 10 CONTINUE
  410. C Fin de la boucle sur les noeuds de IPT1
  411. C On a cree IPT2 maillage transforme de IPT2 (elements,noeuds)
  412.  
  413. NBPTS=NXCOUR/idimp1
  414. SEGADJ MCOORD
  415.  
  416. RETURN
  417. END
  418.  
  419.  
  420.  
  421.  
  422.  
  423.  
  424.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales