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
|
*DECK KINT05
SUBROUTINE KINT05(NBMIX,NEL,LL4,SGD,VOL,MAT,F2,F3)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Multiplication of a matrix by a vector in mesh centered finite
* difference approximation (hexagonal geometry). Special version for
* Trivac.
*
*Copyright:
* Copyright (C) 2010 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. Hebert
*
*Parameters: input
* NBMIX maximum number of material mixtures.
* NEL total number of finite elements.
* LL4 number of unknowns (order of the system matrices).
* SGD cross section per material mixture.
* VOL volumes.
* MAT index-number of the mixture type assigned to each volume.
* F2 vector to multiply.
*
*Parameters: output
* F3 result of the multiplication.
*
*-----------------------------------------------------------------------
*
*----
* SUBROUTINE ARGUMENTS
*----
INTEGER NBMIX,NEL,LL4,MAT(NEL)
REAL SGD(NBMIX),VOL(NEL),F2(LL4),F3(LL4)
*
KEL=0
DO 10 K=1,NEL
L=MAT(K)
IF(L.EQ.0) GO TO 10
KEL=KEL+1
VOL0=VOL(K)
IF(VOL0.EQ.0.0) GO TO 10
F3(KEL)=F3(KEL)+SGD(L)*VOL0*F2(KEL)
10 CONTINUE
RETURN
END
|
