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
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
*DECK TRIMTD
SUBROUTINE TRIMTD(ISPLH,MAXMIX,NEL,LL4,VOL,MAT,SGD,KN,IPW,VEC)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Assembly of a diagonal system matrix for a mesh centered finite
* difference discretization in hexagonal geometry (triangular submeshs).
* Note: system matrix should be initialized by the calling program.
*
*Copyright:
* Copyright (C) 2002 Ecole Polytechnique de Montreal
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version
*
*Author(s): A. Benaboud
*
*Parameters: input
* ISPLH related to the triangular submesh. The number of triangles is
* 6*(ISPLH-1)**2.
* MAXMIX size of array SGD.
* NEL total number of finite elements.
* LL4 order of the system matrices.
* VOL volume of each element.
* MAT mixture index assigned to each hexagon.
* SGD nuclear properties per material mixtures.
* KN element-ordered unknown list.
* IPW permutation matrices.
*
*Parameters: output
* VEC diagonal system matrix.
*
*-----------------------------------------------------------------------
*
*----
* SUBROUTINE ARGUMENTS
*----
INTEGER ISPLH,MAXMIX,NEL,LL4,MAT(NEL),KN(NEL*(18*(ISPLH-1)**2+8)),
1 IPW(LL4)
REAL VOL(NEL),SGD(MAXMIX),VEC(LL4)
*----
* ASSEMBLY OF DIAGONAL MATRIX VEC
*----
NUM1 = 0
NTPH = 6 * (ISPLH-1)**2
NTPL = 1 + 2 * (ISPLH-1)
NVT1 = NTPL + 2 * (ISPLH-2) + NTPH / 2
NVT2 = NTPH - NTPL - (ISPLH-4) * (NTPL+2)
NVT3 = NTPH - (ISPLH-4) * NTPL
IVAL = 3*NTPH+8
IF(ISPLH.EQ.3) NVT2 = NTPH
IF(ISPLH.LE.3) ISAU = 2*(ISPLH-2)
IF(ISPLH.GE.4) ISAU = 6*(ISPLH-3)
ICR = ISAU*(1+2*(ISPLH-2))
DO 40 K=1,NEL
L = MAT(K)
IF(L.EQ.0) GO TO 40
VOL0 = VOL(K)/NTPH
IF(VOL0.EQ.0.0) GO TO 30
DO 20 I = 1,NTPH
*
CALL TRINEI (3,1,1,ISPLH,ICR,I,KK1,KK2,KK3,KEL,IQF,NUM1,
> NTPH,NTPL,NVT1,NVT2,NVT3,IVAL,KN)
*
IND1=IPW(KEL)
VEC(IND1)=VEC(IND1)+SGD(L)*VOL0
20 CONTINUE
30 NUM1=NUM1+IVAL
40 CONTINUE
RETURN
END
|
