| /phasta/phSolver/AMG/ |
| H A D | ramg_smoother.f | 12 integer,intent(in),dimension(nshg+1) :: colm 13 integer,intent(in),dimension(nnz_tot) :: rowp 14 real(kind=8),intent(in),dimension(9,nnz_tot) :: lhsK 15 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 16 integer,intent(in),dimension(nlwork) :: ilwork 17 integer,intent(in),dimension(nshg) :: iBC,iper 18 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 20 real(kind=8),intent(in),dimension(nshg) :: v 62 integer,intent(in) :: Anshg,Annz_tot 63 integer,intent(in),dimension(Anshg+1) :: Acolm [all …]
|
| H A D | ramg_driver.f | 25 integer, intent(in) :: clevel 29 integer,intent(in),dimension(nshg+1) :: colm 30 integer,intent(in),dimension(nnz_tot) :: rowp 31 real(kind=8),intent(in),dimension(9,nnz_tot) :: lhsK 32 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 33 integer,intent(in),dimension(nlwork) :: ilwork 34 integer,intent(in),dimension(nshg) :: iBC,iper 35 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 113 integer,intent(in) :: level 114 real(kind=8),intent(in),dimension(amg_nshg(level)) :: rhs [all …]
|
| H A D | ramg_mls.f | 18 integer,intent(in) :: level 19 integer,intent(in) :: flagfb 20 integer,intent(in),dimension(nshg+1) :: colm 21 integer,intent(in),dimension(nnz_tot) :: rowp 23 real(kind=8),intent(in),dimension(9,nnz_tot) :: lhsK 24 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 26 integer,intent(in),dimension(nlwork) :: ilwork 27 integer,intent(in),dimension(nshg) :: iBC,iper 28 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 198 integer,intent(in) :: level [all …]
|
| H A D | ramg_ggb.f | 10 integer,intent(in),dimension(nshg+1) :: colm 11 integer,intent(in),dimension(nnz_tot) :: rowp 13 real(kind=8),intent(in),dimension(9,nnz_tot) :: lhsK 14 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 16 integer,intent(in),dimension(nlwork) :: ilwork 17 integer,intent(in),dimension(nshg) :: iBC,iper 18 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 202 integer,intent(in),dimension(nlwork) :: ilwork 203 integer,intent(in) :: asize 205 integer,intent(in) :: level [all …]
|
| H A D | ramg_cheby.f | 16 integer,intent(in),dimension(nlwork) :: ilwork 17 integer,intent(in),dimension(nshg) :: iBC,iper 18 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 19 integer,intent(in),dimension(nshg+1) :: colm 20 integer,intent(in),dimension(nnz_tot) :: rowp 22 real(kind=8),intent(in),dimension(9,nnz_tot) :: lhsK 23 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 26 real(kind=8),intent(in),dimension(amg_nshg(level)) :: v,r 28 integer,intent(in) :: fwdbck 113 integer,intent(in),dimension(nshg+1) :: colm [all …]
|
| H A D | ramg_tools.f | 28 integer,intent(in) :: level1,level2 30 real(kind=8),intent(in),dimension(amg_nshg(level2)) :: VC 32 integer,intent(in),dimension(nlwork) :: ilwork 33 integer,intent(in),dimension(nshg) :: iBC,iper 34 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 68 integer,intent(in) :: level1,level2 69 real(kind=8),intent(in),dimension(amg_nshg(level1)) :: vf 109 integer,intent(in) :: gnshg,gnnz_tot 110 integer,intent(in),dimension(gnshg+1) :: gcolm 111 integer,intent(in),dimension(gnnz_tot) :: growp [all …]
|
| H A D | ramg_ITAI.f | 8 integer,intent(in) :: m,l,n 9 integer,intent(in) :: a_colm(m+1),b_colm(l+1) 11 integer,intent(in) :: a_rowp(a_colm(m+1)-1) 12 integer,intent(in) :: b_rowp(b_colm(l+1)-1) 38 integer,intent(in) :: m,l,n 39 integer,intent(in) :: a_colm(m+1),b_colm(l+1) 40 integer,intent(in) :: c_colm(m+1) 42 integer,intent(in) :: a_rowp(a_colm(m+1)-1) 43 integer,intent(in) :: b_rowp(b_colm(l+1)-1) 44 logical,intent(in) :: diagflag [all …]
|
| H A D | ramg_control.f | 12 integer,intent(in),dimension(nshg+1) :: colm 13 integer,intent(in),dimension(nnz_tot) :: rowp 14 real(kind=8),intent(in),dimension(9,nnz_tot) :: lhsK 15 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 17 integer, intent(in), dimension(nlwork) :: ilwork 18 integer, intent(in),dimension(nshg) :: iBC,iper 19 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 73 integer,intent(in),dimension(nlwork) :: ilwork 74 integer,intent(in),dimension(nshg) :: iBC,iper 75 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC [all …]
|
| H A D | ramg_paratools.f | 17 integer,intent(in) :: annz_tot 18 integer,intent(in),dimension(nshg+1) :: acolm 19 integer,intent(in),dimension(annz_tot) :: arowp 20 integer,intent(in) :: redun 21 real(kind=8),intent(in),dimension(redun,annz_tot) :: alhsP 27 integer,intent(in),dimension(nlwork) :: ilwork 28 integer,intent(in),dimension(nshg) :: iBC,iper 29 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 493 real(kind=8),intent(in),dimension(nshg) :: p 495 integer,intent(in),dimension(nshg+1) :: colm [all …]
|
| H A D | amgread.txt | 5 Poisson Equation in PHASTA. If you want to start to modify or add 13 you continue. row-wise sparse storage is used everywhere in this 20 *** Files related to AMG in the ../common/ and ../incompressible/ 39 **** in AMG directory ****** 46 every function call (though it's still tedious in some places). 49 dimentions. This is used, as shown in this file later, in submatrices 52 Then we have our multilevel AMG stored in amg_A_* arrays as sparse 54 are same as ilwork and nlwork, but in higher levels with smaller size. 63 CF-spliting is done in a seperate function ramg_CFsplit. Other than the 75 "coarse first, fine second". This order will be used in Gauss-Seidel [all …]
|
| H A D | ramg_coarse.f | 21 integer,intent(in) :: level1,level2 22 real(kind=8),intent(in) :: eps_str 23 integer,intent(in) :: interp_flag 24 real(kind=8),intent(in) :: trunc_tol 26 integer,intent(in),dimension(nlwork) :: ilwork 27 integer,intent(in),dimension(nshg) :: iBC,iper 28 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC 471 real(kind=8),dimension(0:ilen-1),intent(in) :: wght 500 integer,intent(in) :: popid,ilen 501 real(kind=8),dimension(0:ilen-1),intent(in) :: wght [all …]
|
| H A D | ramg_extract.f | 28 integer,intent(in),dimension(nshg+1) :: colm 29 integer,intent(in),dimension(nnz_tot) :: rowp 30 real(kind=8),intent(in),dimension(9,nnz_tot) :: lhsK 31 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 38 integer, intent(in), dimension(nlwork) :: ilwork 39 integer, intent(in),dimension(nshg) :: iBC,iper 40 real(kind=8),intent(in),dimension(nshg,ndofBC) :: BC
|
| /phasta/phSolver/common/ |
| H A D | phIO.f | 8 character(c_char), intent(in) :: fname(*) 25 character(c_char), intent(in) :: phrase(*) 27 integer(c_int), intent(in) :: nvals 28 character(c_char), intent(in) :: datatype(*) 29 character(c_char), intent(in) :: iotype(*) 38 character(c_char), intent(in) :: phrase(*) 40 integer(c_int), intent(in) :: nitems 41 integer(c_int), intent(in) :: ndata 42 character(c_char), intent(in) :: datatype(*) 43 character(c_char), intent(in) :: iotype(*) [all …]
|
| H A D | scatnu.f | 13 do in = 1,nshl 15 xmudmi(:,int) = xmudmi(:,int) + xnut(ien(:,in)) * strl(:,int) 16 & *shp(in,int)
|
| H A D | phasta.h | 8 int phasta(phSolver::Input& ctrl, GRStream* in); 10 int phasta(phSolver::Input& ctrl, GRStream* in, RStream* out);
|
| H A D | phStr.f | 8 integer(c_int), value, intent(in) :: val 17 real(c_double), value, intent(in) :: val
|
| /phasta/phSolver/incompressible/ |
| H A D | bc3lhs.f | 46 in = abs(ien(iel,inod)) 47 if (ibits(iBC(in),3,3) .eq. 0) goto 5000 ! NO velocity BC's 48 if (ibits(iBC(in),3,3) .eq. 7) goto 5000 ! 3 components ok 55 if ( ibits(iBC(in),3,3) .eq. 1) then 85 & - BC(in,4) * xKebe(iel,irem1,i,inod) 87 & - BC(in,4) * xKebe(iel,irem2,i,inod) 89 & - BC(in,4) * xKebe(iel,irem3,i,inod) 105 & - BC(in,5) * xKebe(iel,irem1,i,inod) 107 & - BC(in,5) * xKebe(iel,irem2,i,inod) 109 & - BC(in,5) * xKebe(iel,irem3,i,inod) [all …]
|
| H A D | bc3global.f | 28 do in=1,nshg 29 i0 = (in-1)*4 33 if ( btest(iBC(in),2) ) then 44 if ( ibits(iBC(in),3,3) .eq. 1 ) then 52 if ( ibits(iBC(in),3,3) .eq. 2 ) then 61 if ( ibits(iBC(in),3,3) .eq. 3 ) then 72 if ( ibits(iBC(in),3,3) .eq. 4 ) then 81 if ( ibits(iBC(in),3,3) .eq. 5 ) then 93 if ( ibits(iBC(in),3,3) .eq. 6 ) then 105 if ( ibits(iBC(in),3,3) .eq. 7 ) then
|
| H A D | getdiff.f | 195 double precision, intent(in), dimension(npro,nshl,ndof) :: 197 double precision, intent(in), dimension(npro,nshl) :: 282 double precision, intent(in), dimension(npro,nshl,ndof) :: 284 double precision, intent(in), dimension(npro,nshl) :: 340 double precision, intent(in), dimension(npro,nshl,ndof) :: 342 double precision, intent(in), dimension(npro,nshl) :: 344 double precision, intent(in), dimension(npro) ::
|
| /phasta/phSolver/compressible/ |
| H A D | bc3lhs.f | 44 in = ien(iel,inod) 45 if (iBC(in) .eq. 0) goto 5000 55 if ( btest(iBC(in),2) ) then 73 if ( ibits(iBC(in),3,3) .eq. 1) then 76 & - BC(in,4) * EGmass(iel,i2,i) 78 & - BC(in,5) * EGmass(iel,i2,i) 82 & - BC(in,4) * EGmass(iel,i,i2) 84 & - BC(in,5) * EGmass(iel,i,i2) 95 if ( ibits(iBC(in),3,3) .eq. 2 ) then 98 & - BC(in,4) * EGmass(iel,i3,i) [all …]
|
| /phasta/converterIO/ |
| H A D | README.md | 21 …e single geombc and restart file per part), where # is the number of parts in the partitioned mesh. 25 …the #-procs_case directory, we simply need to run the converterO2N executable in the following way: 29 which will read the PHASTA files in the #-procs_case directory and generate a new directory named #… 39 - the number of part in the mesh (i.e. #), 42 - the number of fields (integer and double) in both the geombc and restart files to convert, and 43 - the list of all the fields in both the geombc and restart files to convert. 50 …files must always preceed the integer field in the list and the fields counter needs also to be in… 52 …n grep, awk, etc and needs to be applied to ALL the phasta files contained in the #-procs_case dir… 53 …there is no garantee that a single geombc file will contain all the topologies present in the mesh. 61 …format to the Posix format, we simply need to run the converterN2O executable in the following way: [all …]
|
| /phasta/phSolver/ |
| H A D | CHANGES | 10 VERSION needs to be specified if the source base is in the dirtory structure 26 …If you are looking to build AMG code you will need to execute the Makefile in that directory betwe… 31 Fixed several places that compressible code was broken in the 3 previous commits. 47 CVS: Committing in . 85 small but essential change in genbkb.f 100 iotype for arrays can now be binary or ascii (set by user in solver.inp). 104 phasta became a function in common that is called by main in either incompressible or compressible … 153 makefiles and as such has a makefile in each of the three 155 or compressible). Many files have changed but almost all in a superficial way 211 3/ mesh dual info now in phasta I/O format [all …]
|
| /phasta/ |
| H A D | License.txt | 4 Redistribution and use in source and binary forms, with or without modification, are permitted prov… 8 2. Redistributions in binary form must reproduce the above copyright notice, this list of condition…
|
| /phasta/svLS/ |
| H A D | svLS_STRUCT.h | 9 ! paragraphs appear in all copies. 57 ! Task in FORTRAN indexing (USE) 74 ! Included in the computations (IN) 96 ! All following are in (USE) 125 ! Number of non-zero in lhs (IN) 171 ! Res. rduction in last itr. (OUT)
|
| /phasta/phSolver/Test_C/ |
| H A D | solver.inp | 9 # To find the Key Phrases, first look in input.config. There you will 12 # in this file. In most cases, the acceptable inputs also appear in 64 Include Viscous Correction in Stabilization: False # if p=1 idiff=1
|