src/modules/DCRComplexSolve.F90

64   INTEGER, POINTER :: NodeIndexes(:)  local
238         NodeIndexes => CurrentElement % NodeIndexes
240         ElementNodes % x(1:n) = Solver % Mesh % Nodes % x(NodeIndexes)
241         ElementNodes % y(1:n) = Solver % Mesh % Nodes % y(NodeIndexes)
242         ElementNodes % z(1:n) = Solver % Mesh % Nodes % z(NodeIndexes)
256              'Amatrix', Material, n, NodeIndexes )
259              'Avector 1', Material, n, NodeIndexes )
262              'Avector 2', Material, n, NodeIndexes )
265              'Ascalar 1', n, NodeIndexes, GotIt)
268              'Ascalar 2', n, NodeIndexes, GotIt)
279                 'Pressure Source 1', n, NodeIndexes, GotIt )
282                 'Pressure Source 2', n, NodeIndexes, GotIt )
297                   PressurePerm(NodeIndexes) )
326               NodeIndexes => CurrentElement % NodeIndexes
328               IF ( ANY( PressurePerm(NodeIndexes) == 0 ) ) CYCLE
330               ElementNodes % x(1:n) = Solver % Mesh % Nodes % x(NodeIndexes)
331               ElementNodes % y(1:n) = Solver % Mesh % Nodes % y(NodeIndexes)
332               ElementNodes % z(1:n) = Solver % Mesh % Nodes % z(NodeIndexes)
335                    'Wave Flux 1', n, NodeIndexes, GotIt )
338                    'Wave Flux 2', n, NodeIndexes, GotIt )
341                    'Bvector', Model % BCs(i) % Values,  n, NodeIndexes )
344                    'Bscalar 1', n, NodeIndexes, GotIt)
347                    'Bscalar 2', n, NodeIndexes, GotIt)
362                       PressurePerm(NodeIndexes) )
426    SUBROUTINE InputTensor( Tensor, IsScalar, Name, Material, n, NodeIndexes )
429       INTEGER :: n, NodeIndexes(:)  local
446       CALL ListGetRealArray( Material, Name, Hwrk, n, NodeIndexes, stat )
478    SUBROUTINE InputVector( Tensor, IsScalar, Name, Material, n, NodeIndexes )
481       INTEGER :: n, NodeIndexes(:)  local
498       CALL ListGetRealArray( Material, Name, Hwrk, n, NodeIndexes, stat )
838      ELSE IF ( ANY( Perm( Element % NodeIndexes ) <= 0 ) ) THEN
843      IF ( ANY( Perm( Element % NodeIndexes ) <= 0 ) ) RETURN
850      EdgeNodes % x = Mesh % Nodes % x(Edge % NodeIndexes)
851      EdgeNodes % y = Mesh % Nodes % y(Edge % NodeIndexes)
852      EdgeNodes % z = Mesh % Nodes % z(Edge % NodeIndexes)
856      Nodes % x = Mesh % Nodes % x(Element % NodeIndexes)
857      Nodes % y = Mesh % Nodes % y(Element % NodeIndexes)
858      Nodes % z = Mesh % Nodes % z(Element % NodeIndexes)
866           IF ( Edge % NodeIndexes(l) == Element % NodeIndexes(k) ) THEN
890              Model % BCs(j) % Values, Pn, Element % NodeIndexes )
893              'Bscalar 1', Pn, Element % NodeIndexes, GotIt)
896              'Bscalar 2', Pn, Element % NodeIndexes, GotIt)
899           'Wave Flux 1', En, Edge % NodeIndexes, gotIt )
903           'Wave Flux 2', En, Edge % NodeIndexes, gotIt )
913         Pressure(1,1:Pn) = Quant( 2*Perm(Element % NodeIndexes)-1 )
914         Pressure(2,1:Pn) = Quant( 2*Perm(Element % NodeIndexes)-0 )
1067    SUBROUTINE InputVector( Tensor, IsScalar, Name, Material, n, NodeIndexes )  argument
1070       INTEGER :: n, NodeIndexes(:)  local
1087       CALL ListGetRealArray( Material, Name, Hwrk, n, NodeIndexes, stat )
1192      EdgeNodes % x = Mesh % Nodes % x(Edge % NodeIndexes)
1193      EdgeNodes % y = Mesh % Nodes % y(Edge % NodeIndexes)
1194      EdgeNodes % z = Mesh % Nodes % z(Edge % NodeIndexes)
1244            IF ( ANY( Perm( Element % NodeIndexes ) <= 0 ) ) CYCLE
1253                  IF ( Edge % NodeIndexes(j) == Element % NodeIndexes(k) ) THEN
1269            Nodes % x(1:Pn) = Mesh % Nodes % x(Element % NodeIndexes)
1270            Nodes % y(1:Pn) = Mesh % Nodes % y(Element % NodeIndexes)
1271            Nodes % z(1:Pn) = Mesh % Nodes % z(Element % NodeIndexes)
1282                 'Amatrix', Model % Materials(k) % Values, Pn, Element % NodeIndexes )
1285                 'Avector 1', Model % Materials(k) % Values, Pn, Element % NodeIndexes )
1288                 'Avector 2', Model % Materials(k) % Values, Pn, Element % NodeIndexes )
1291                 'Ascalar 1', Pn, Element % NodeIndexes, GotIt)
1294                 'Ascalar 2', Pn, Element % NodeIndexes, GotIt)
1315               Pressure(k,1:Pn) = Quant( 2*Perm(Element % NodeIndexes)-2+k )
1367    SUBROUTINE InputTensor( Tensor, IsScalar, Name, Material, n, NodeIndexes )  argument
1370       INTEGER :: n, NodeIndexes(:)  local
1387       CALL ListGetRealArray( Material, Name, Hwrk, n, NodeIndexes, stat )
1419    SUBROUTINE InputVector( Tensor, IsScalar, Name, Material, n, NodeIndexes )  argument
1422       INTEGER :: n, NodeIndexes(:)  local
1439       CALL ListGetRealArray( Material, Name, Hwrk, n, NodeIndexes, stat )
1512      IF ( ANY( Perm( Element % NodeIndexes ) <= 0 ) ) RETURN
1531      Nodes % x = Mesh % Nodes % x(Element % NodeIndexes)
1532      Nodes % y = Mesh % Nodes % y(Element % NodeIndexes)
1533      Nodes % z = Mesh % Nodes % z(Element % NodeIndexes)
1541      Pressione(1,1:n) = Quant( 2*Perm(Element % NodeIndexes)-1 )
1542      Pressione(2,1:n) = Quant( 2*Perm(Element % NodeIndexes)-0 )
1552           'Amatrix', Material, n, Element % NodeIndexes )
1555           'Avector 1', Material, n, Element % NodeIndexes )
1558           'Avector 2', Material, n, Element % NodeIndexes )
1561           'Ascalar 1', n, Element % NodeIndexes, GotIt)
1564           'Ascalar 2', n, Element % NodeIndexes, GotIt)
1575              'Pressure Source 1', n, Element % NodeIndexes, GotIt )
1578              'Pressure Source 2', n, Element % NodeIndexes, GotIt )
1720    SUBROUTINE InputTensor( Tensor, IsScalar, Name, Material, n, NodeIndexes )  argument
1723       INTEGER :: n, NodeIndexes(:)  local
1740       CALL ListGetRealArray( Material, Name, Hwrk, n, NodeIndexes, stat )
1772    SUBROUTINE InputVector( Tensor, IsScalar, Name, Material, n, NodeIndexes )  argument
1775       INTEGER :: n, NodeIndexes(:)  local
1792       CALL ListGetRealArray( Material, Name, Hwrk, n, NodeIndexes, stat )