diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Trivac/src/KINB04.f | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Trivac/src/KINB04.f')
| -rwxr-xr-x | Trivac/src/KINB04.f | 66 |
1 files changed, 66 insertions, 0 deletions
diff --git a/Trivac/src/KINB04.f b/Trivac/src/KINB04.f new file mode 100755 index 0000000..0517873 --- /dev/null +++ b/Trivac/src/KINB04.f @@ -0,0 +1,66 @@ +*DECK KINB04 + SUBROUTINE KINB04(MAXKN,MAXQF,SGD,NREG,LL4,ISPLH,NBMIX,MAT,KN, + 1 QFR,VOL,F2,F3) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Multiplication of a matrix by a vector in mesh-centered finite- +* difference diffusion approximation (hexagonal geometry). Special +* version for Bivac. +* +*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 +* MAXKN dimension of array KN. +* MAXQF dimension of array QFR. +* SGD mixture-ordered cross sections. +* NREG number of hexagons in Bivac. +* LL4 number of unknowns per group in Bivac. Equal to the number +* of finite elements (hexagons or triangles) excluding the +* virtual elements. +* ISPLH type of hexagonal mesh-splitting: +* =1: hexagonal elements; >1: triangular elements. +* NBMIX number of macro-mixtures. +* MAT mixture index per hexagon. +* KN element-ordered unknown list. +* QFR element-ordered information. +* VOL volume of hexagons. +* F2 vector to multiply. +* +*Parameters: output +* F3 result of the multiplication. +* +*----------------------------------------------------------------------- +* +*---- +* SUBROUTINE ARGUMENTS +*---- + INTEGER MAXKN,MAXQF,NREG,LL4,ISPLH,NBMIX,MAT(NREG),KN(MAXKN) + REAL SGD(NBMIX),QFR(MAXQF),VOL(NREG),F2(LL4),F3(LL4) +* + IF(ISPLH.EQ.1) THEN + NSURF=6 + ELSE + NSURF=3 + ENDIF +*---- +* MULTIPLICATION. +*---- + NUM1=0 + DO 20 IND1=1,LL4 + KHEX=KN(NUM1+NSURF+1) + IF(VOL(KHEX).EQ.0.0) GO TO 10 + L=MAT(KHEX) + F3(IND1)=F3(IND1)+SGD(L)*QFR(NUM1+NSURF+1)*F2(IND1) + 10 NUM1=NUM1+NSURF+1 + 20 CONTINUE + RETURN + END |