| bps.usr (86a4271f02baafab377de162285a9b890a25166a) | bps.usr (e2b2c771c7282b3b06c9e412ce543de8dc2e42e1) |
|---|---|
| 1C Copyright (c) 2017, Lawrence Livermore National Security, LLC. Produced at 2C the Lawrence Livermore National Laboratory. LLNL-CODE-734707. All Rights 3C reserved. See files LICENSE and NOTICE for details. 4C 5C This file is part of CEED, a collection of benchmarks, miniapps, software 6C libraries and APIs for efficient high-order finite element and spectral 7C element discretizations for exascale applications. For more information and 8C source code availability see http://github.com/ceed. --- 699 unchanged lines hidden (view full) --- 708 parameter (lt=lx1*ly1*lz1*lelt) 709 parameter (ld=lxd*lyd*lzd*lelt) 710 common /vcrns/ u1(lt),r1(lt),r2(lt),r3(lt) 711 common /vcrny/ e1(lt) 712 common /vcrvh/ h1(lt),h2(lx*lelt),pap(3) 713 real*8 coords(ldim*lx*lelt) 714 715 logical ifh3 | 1C Copyright (c) 2017, Lawrence Livermore National Security, LLC. Produced at 2C the Lawrence Livermore National Laboratory. LLNL-CODE-734707. All Rights 3C reserved. See files LICENSE and NOTICE for details. 4C 5C This file is part of CEED, a collection of benchmarks, miniapps, software 6C libraries and APIs for efficient high-order finite element and spectral 7C element discretizations for exascale applications. For more information and 8C source code availability see http://github.com/ceed. --- 699 unchanged lines hidden (view full) --- 708 parameter (lt=lx1*ly1*lz1*lelt) 709 parameter (ld=lxd*lyd*lzd*lelt) 710 common /vcrns/ u1(lt),r1(lt),r2(lt),r3(lt) 711 common /vcrny/ e1(lt) 712 common /vcrvh/ h1(lt),h2(lx*lelt),pap(3) 713 real*8 coords(ldim*lx*lelt) 714 715 logical ifh3 |
| 716 integer*8 ndof | 716 integer*8 nnode |
| 717 integer ceed,err,test 718 character*64 spec 719 | 717 integer ceed,err,test 718 character*64 spec 719 |
| 720 integer p,q,ncomp,edof,ldof | 720 integer p,q,ncomp,enode,lnode |
| 721 integer vec_p1,vec_ap1,vec_qdata,vec_coords,vec_rhs 722 integer erstrctu,erstrctx,erstrctw 723 integer basisu,basisx 724 integer qf_mass,qf_setup 725 integer op_mass,op_setup 726 real*8 x,y,z 727 integer*8 offset 728 --- 45 unchanged lines hidden (view full) --- 774 q=p+1 775 ncomp=1 776 call ceedbasiscreatetensorh1lagrange(ceed,ndim,ndim,p,q, 777 $ ceed_gauss,basisx,err) 778 call ceedbasiscreatetensorh1lagrange(ceed,ndim,ncomp,p,q, 779 $ ceed_gauss,basisu,err) 780 781C Create ceed element restrictions for mesh and computation | 721 integer vec_p1,vec_ap1,vec_qdata,vec_coords,vec_rhs 722 integer erstrctu,erstrctx,erstrctw 723 integer basisu,basisx 724 integer qf_mass,qf_setup 725 integer op_mass,op_setup 726 real*8 x,y,z 727 integer*8 offset 728 --- 45 unchanged lines hidden (view full) --- 774 q=p+1 775 ncomp=1 776 call ceedbasiscreatetensorh1lagrange(ceed,ndim,ndim,p,q, 777 $ ceed_gauss,basisx,err) 778 call ceedbasiscreatetensorh1lagrange(ceed,ndim,ncomp,p,q, 779 $ ceed_gauss,basisu,err) 780 781C Create ceed element restrictions for mesh and computation |
| 782 edof=p**ldim 783 ldof=edof*nelt*ncomp 784 call ceedelemrestrictioncreateidentity(ceed,nelt,edof,ldof, | 782 enode=p**ldim 783 lnode=enode*nelt*ncomp 784 call ceedelemrestrictioncreateidentity(ceed,nelt,enode,lnode, |
| 785 $ ldim,erstrctx,err) | 785 $ ldim,erstrctx,err) |
| 786 call ceedelemrestrictioncreateidentity(ceed,nelt,edof,ldof, | 786 call ceedelemrestrictioncreateidentity(ceed,nelt,enode,lnode, |
| 787 $ ncomp,erstrctu,err) 788 call ceedelemrestrictioncreateidentity(ceed,nelt,q**ldim, 789 $ nelt*q**ldim,1,erstrctw,err) 790 791C Create ceed vectors | 787 $ ncomp,erstrctu,err) 788 call ceedelemrestrictioncreateidentity(ceed,nelt,q**ldim, 789 $ nelt*q**ldim,1,erstrctw,err) 790 791C Create ceed vectors |
| 792 call ceedvectorcreate(ceed,ldof,vec_p1,err) 793 call ceedvectorcreate(ceed,ldof,vec_ap1,err) 794 call ceedvectorcreate(ceed,ldof,vec_rhs,err) | 792 call ceedvectorcreate(ceed,lnode,vec_p1,err) 793 call ceedvectorcreate(ceed,lnode,vec_ap1,err) 794 call ceedvectorcreate(ceed,lnode,vec_rhs,err) |
| 795 call ceedvectorcreate(ceed,ldim*lx*nelt,vec_coords,err) 796 call ceedvectorcreate(ceed,nelt*q**ldim,vec_qdata,err) 797 798 offset=0 799 call ceedvectorsetarray(vec_coords,ceed_mem_host, 800 $ ceed_use_pointer,coords,offset,err) 801 802C Create ceed qfunctions for masssetupf and massf --- 88 unchanged lines hidden (view full) --- 891 write(6,*) "UNCONVERGED CG" 892 endif 893 if (dabs(er1)>5e-3) then 894 write(6,*) "ERROR IS TOO LARGE" 895 endif 896 endif 897 898 nx = nx1-1 | 795 call ceedvectorcreate(ceed,ldim*lx*nelt,vec_coords,err) 796 call ceedvectorcreate(ceed,nelt*q**ldim,vec_qdata,err) 797 798 offset=0 799 call ceedvectorsetarray(vec_coords,ceed_mem_host, 800 $ ceed_use_pointer,coords,offset,err) 801 802C Create ceed qfunctions for masssetupf and massf --- 88 unchanged lines hidden (view full) --- 891 write(6,*) "UNCONVERGED CG" 892 endif 893 if (dabs(er1)>5e-3) then 894 write(6,*) "ERROR IS TOO LARGE" 895 endif 896 endif 897 898 nx = nx1-1 |
| 899 ndof = nelgt ! ndofs 900 ndof = ndof*(nx**ldim) ! ndofs 901 nppp = ndof/np ! ndofs/proc | 899 nnode = nelgt ! nnodes 900 nnode = nnode*(nx**ldim) ! nnodes 901 nppp = nnode/np ! nnodes/proc |
| 902 | 902 |
| 903 dofpss = ndof/telaps ! DOF/sec - scalar form | 903 nodepss = nnode/telaps ! DOF/sec - scalar form |
| 904 titers = telaps/maxits ! time per iteration 905 tppp_s = titers/nppp ! time per iteraton per local point 906 907 if (nid.eq.0) write(6,1) 'case scalar:' | 904 titers = telaps/maxits ! time per iteration 905 tppp_s = titers/nppp ! time per iteraton per local point 906 907 if (nid.eq.0) write(6,1) 'case scalar:' |
| 908 $ ,np,nx,nelt,nelgt,ndof,nppp,maxits,telaps,dofpss,titers,tppp_s | 908 $ ,np,nx,nelt,nelgt,nnode,nppp,maxits,telaps,nodepss,titers,tppp_s |
| 909 910C Solve libCEED algebraic RHS 911 if (nid.eq.0) write (6,*) "libCEED algebraic RHS" 912 maxit = 100 913 tstart = dnekclock() 914 call cggos(u1,r2,h1,h2,vmult,binvm1,tol,ceed,op_mass, 915 $ vec_p1,vec_ap1,maxit,'bp1') 916 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 926 write(6,*) "UNCONVERGED CG" 927 endif 928 if (dabs(er1)>1e-5) then 929 write(6,*) "ERROR IS TOO LARGE" 930 endif 931 endif 932 933 nx = nx1-1 | 909 910C Solve libCEED algebraic RHS 911 if (nid.eq.0) write (6,*) "libCEED algebraic RHS" 912 maxit = 100 913 tstart = dnekclock() 914 call cggos(u1,r2,h1,h2,vmult,binvm1,tol,ceed,op_mass, 915 $ vec_p1,vec_ap1,maxit,'bp1') 916 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 926 write(6,*) "UNCONVERGED CG" 927 endif 928 if (dabs(er1)>1e-5) then 929 write(6,*) "ERROR IS TOO LARGE" 930 endif 931 endif 932 933 nx = nx1-1 |
| 934 ndof = nelgt ! ndofs 935 ndof = ndof*(nx**ldim) ! ndofs 936 nppp = ndof/np ! ndofs/proc | 934 nnode = nelgt ! nnodes 935 nnode = nnode*(nx**ldim) ! nnodes 936 nppp = nnode/np ! nnodes/proc |
| 937 | 937 |
| 938 dofpss = ndof/telaps ! DOF/sec - scalar form | 938 nodepss = nnode/telaps ! DOF/sec - scalar form |
| 939 titers = telaps/maxits ! time per iteration 940 tppp_s = titers/nppp ! time per iteraton per local point 941 942 if (nid.eq.0) write(6,1) 'case scalar:' | 939 titers = telaps/maxits ! time per iteration 940 tppp_s = titers/nppp ! time per iteraton per local point 941 942 if (nid.eq.0) write(6,1) 'case scalar:' |
| 943 $ ,np,nx,nelt,nelgt,ndof,nppp,maxits,telaps,dofpss,titers,tppp_s | 943 $ ,np,nx,nelt,nelgt,nnode,nppp,maxits,telaps,nodepss,titers,tppp_s |
| 944 945C Solve Nek5000 algebraic RHS 946 if (nid.eq.0) write (6,*) "Nek5000 algebraic RHS" 947 maxit = 100 948 tstart = dnekclock() 949 call cggos(u1,r3,h1,h2,vmult,binvm1,tol,ceed,op_mass, 950 $ vec_p1,vec_ap1,maxit,'bp1') 951 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 961 write(6,*) "UNCONVERGED CG" 962 endif 963 if (dabs(er1)>1e-5) then 964 write(6,*) "ERROR IS TOO LARGE" 965 endif 966 endif 967 968 nx = nx1-1 | 944 945C Solve Nek5000 algebraic RHS 946 if (nid.eq.0) write (6,*) "Nek5000 algebraic RHS" 947 maxit = 100 948 tstart = dnekclock() 949 call cggos(u1,r3,h1,h2,vmult,binvm1,tol,ceed,op_mass, 950 $ vec_p1,vec_ap1,maxit,'bp1') 951 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 961 write(6,*) "UNCONVERGED CG" 962 endif 963 if (dabs(er1)>1e-5) then 964 write(6,*) "ERROR IS TOO LARGE" 965 endif 966 endif 967 968 nx = nx1-1 |
| 969 ndof = nelgt ! ndofs 970 ndof = ndof*(nx**ldim) ! ndofs 971 nppp = ndof/np ! ndofs/proc | 969 nnode = nelgt ! nnodes 970 nnode = nnode*(nx**ldim) ! nnodes 971 nppp = nnode/np ! nnodes/proc |
| 972 | 972 |
| 973 dofpss = ndof/telaps ! DOF/sec - scalar form | 973 nodepss = nnode/telaps ! DOF/sec - scalar form |
| 974 titers = telaps/maxits ! time per iteration 975 tppp_s = titers/nppp ! time per iteraton per local point 976 977 if (nid.eq.0) write(6,1) 'case scalar:' | 974 titers = telaps/maxits ! time per iteration 975 tppp_s = titers/nppp ! time per iteraton per local point 976 977 if (nid.eq.0) write(6,1) 'case scalar:' |
| 978 $ ,np,nx,nelt,nelgt,ndof,nppp,maxits,telaps,dofpss,titers,tppp_s | 978 $ ,np,nx,nelt,nelgt,nnode,nppp,maxits,telaps,nodepss,titers,tppp_s |
| 979 980 1 format(a12,i7,i3,i7,i10,i14,i10,i4,1p4e13.5) 981 3 format(i3,i9,e12.4,1x,a8,i9) 982 983C Destroy ceed handles 984 call ceedvectordestroy(vec_p1,err) 985 call ceedvectordestroy(vec_ap1,err) 986 call ceedvectordestroy(vec_rhs,err) --- 28 unchanged lines hidden (view full) --- 1015 parameter (lt=lx1*ly1*lz1*lelt) 1016 parameter (ld=lxd*lyd*lzd*lelt) 1017 common /vcrns/ u1(lt),r1(lt),r2(lt),r3(lt) 1018 common /vcrny/ e1(lt) 1019 common /vcrvh/ h1(lt),h2(ld),pap(3) 1020 real*8 coords(ldim*lx*lelt) 1021 1022 logical ifh3 | 979 980 1 format(a12,i7,i3,i7,i10,i14,i10,i4,1p4e13.5) 981 3 format(i3,i9,e12.4,1x,a8,i9) 982 983C Destroy ceed handles 984 call ceedvectordestroy(vec_p1,err) 985 call ceedvectordestroy(vec_ap1,err) 986 call ceedvectordestroy(vec_rhs,err) --- 28 unchanged lines hidden (view full) --- 1015 parameter (lt=lx1*ly1*lz1*lelt) 1016 parameter (ld=lxd*lyd*lzd*lelt) 1017 common /vcrns/ u1(lt),r1(lt),r2(lt),r3(lt) 1018 common /vcrny/ e1(lt) 1019 common /vcrvh/ h1(lt),h2(ld),pap(3) 1020 real*8 coords(ldim*lx*lelt) 1021 1022 logical ifh3 |
| 1023 integer*8 ndof | 1023 integer*8 nnode |
| 1024 integer ceed,err,test 1025 character*64 spec 1026 | 1024 integer ceed,err,test 1025 character*64 spec 1026 |
| 1027 integer p,q,ncomp,edof,ldof | 1027 integer p,q,ncomp,enode,lnode |
| 1028 integer vec_p1,vec_ap1,vec_qdata,vec_coords,vec_rhs 1029 integer erstrctu,erstrctx,erstrctw 1030 integer basisu,basisx 1031 integer qf_diffusion,qf_setup 1032 integer op_diffusion,op_setup 1033 integer ii,i,ngeo 1034 real*8 x,y,z 1035 integer*8 offset --- 53 unchanged lines hidden (view full) --- 1089 q=p+1 1090 ncomp=1 1091 call ceedbasiscreatetensorh1lagrange(ceed,ldim,3*ncomp,p,q, 1092 $ ceed_gauss,basisx,err) 1093 call ceedbasiscreatetensorh1lagrange(ceed,ldim,ncomp,p,q, 1094 $ ceed_gauss,basisu,err) 1095 1096C Create ceed element restrictions for mesh and computation | 1028 integer vec_p1,vec_ap1,vec_qdata,vec_coords,vec_rhs 1029 integer erstrctu,erstrctx,erstrctw 1030 integer basisu,basisx 1031 integer qf_diffusion,qf_setup 1032 integer op_diffusion,op_setup 1033 integer ii,i,ngeo 1034 real*8 x,y,z 1035 integer*8 offset --- 53 unchanged lines hidden (view full) --- 1089 q=p+1 1090 ncomp=1 1091 call ceedbasiscreatetensorh1lagrange(ceed,ldim,3*ncomp,p,q, 1092 $ ceed_gauss,basisx,err) 1093 call ceedbasiscreatetensorh1lagrange(ceed,ldim,ncomp,p,q, 1094 $ ceed_gauss,basisu,err) 1095 1096C Create ceed element restrictions for mesh and computation |
| 1097 edof=p**ldim 1098 ldof=edof*nelt*ncomp | 1097 enode=p**ldim 1098 lnode=enode*nelt*ncomp |
| 1099 ngeo=(ldim*(ldim+1))/2 | 1099 ngeo=(ldim*(ldim+1))/2 |
| 1100 call ceedelemrestrictioncreateidentity(ceed,nelt,edof,ldof, | 1100 call ceedelemrestrictioncreateidentity(ceed,nelt,enode,lnode, |
| 1101 $ ldim,erstrctx,err) | 1101 $ ldim,erstrctx,err) |
| 1102 call ceedelemrestrictioncreateidentity(ceed,nelt,edof,ldof, | 1102 call ceedelemrestrictioncreateidentity(ceed,nelt,enode,lnode, |
| 1103 $ ncomp,erstrctu,err) 1104 call ceedelemrestrictioncreateidentity(ceed,nelt,q**ldim, 1105 $ nelt*q**ldim,ngeo,erstrctw,err) 1106 1107C Create ceed vectors | 1103 $ ncomp,erstrctu,err) 1104 call ceedelemrestrictioncreateidentity(ceed,nelt,q**ldim, 1105 $ nelt*q**ldim,ngeo,erstrctw,err) 1106 1107C Create ceed vectors |
| 1108 call ceedvectorcreate(ceed,ldof,vec_p1,err) 1109 call ceedvectorcreate(ceed,ldof,vec_ap1,err) 1110 call ceedvectorcreate(ceed,ldof,vec_rhs,err) | 1108 call ceedvectorcreate(ceed,lnode,vec_p1,err) 1109 call ceedvectorcreate(ceed,lnode,vec_ap1,err) 1110 call ceedvectorcreate(ceed,lnode,vec_rhs,err) |
| 1111 call ceedvectorcreate(ceed,ldim*lx*nelt,vec_coords,err) 1112 call ceedvectorcreate(ceed,nelt*ngeo*q**ldim,vec_qdata,err) 1113 1114 offset=0 1115 call ceedvectorsetarray(vec_coords,ceed_mem_host, 1116 $ ceed_use_pointer,coords,offset,err) 1117 1118C Create ceed qfunctions for diffsetupf and diffusionf --- 91 unchanged lines hidden (view full) --- 1210 write(6,*) "UNCONVERGED CG" 1211 endif 1212 if (dabs(er1)>1e-3) then 1213 write(6,*) "ERROR IS TOO LARGE" 1214 endif 1215 endif 1216 1217 nx = nx1-1 | 1111 call ceedvectorcreate(ceed,ldim*lx*nelt,vec_coords,err) 1112 call ceedvectorcreate(ceed,nelt*ngeo*q**ldim,vec_qdata,err) 1113 1114 offset=0 1115 call ceedvectorsetarray(vec_coords,ceed_mem_host, 1116 $ ceed_use_pointer,coords,offset,err) 1117 1118C Create ceed qfunctions for diffsetupf and diffusionf --- 91 unchanged lines hidden (view full) --- 1210 write(6,*) "UNCONVERGED CG" 1211 endif 1212 if (dabs(er1)>1e-3) then 1213 write(6,*) "ERROR IS TOO LARGE" 1214 endif 1215 endif 1216 1217 nx = nx1-1 |
| 1218 ndof = nelgt ! ndofs 1219 ndof = ndof*(nx**ldim) ! ndofs 1220 nppp = ndof/np ! ndofs/proc | 1218 nnode = nelgt ! nnodes 1219 nnode = nnode*(nx**ldim) ! nnodes 1220 nppp = nnode/np ! nnodes/proc |
| 1221 | 1221 |
| 1222 dofpss = ndof/telaps ! DOF/sec - scalar form | 1222 nodepss = nnode/telaps ! DOF/sec - scalar form |
| 1223 titers = telaps/maxits ! time per iteration 1224 tppp_s = titers/nppp ! time per iteraton per local point 1225 1226 if (nid.eq.0) write(6,1) 'case scalar:' | 1223 titers = telaps/maxits ! time per iteration 1224 tppp_s = titers/nppp ! time per iteraton per local point 1225 1226 if (nid.eq.0) write(6,1) 'case scalar:' |
| 1227 $ ,np,nx,nelt,nelgt,ndof,nppp,maxits,telaps,dofpss,titers,tppp_s | 1227 $ ,np,nx,nelt,nelgt,nnode,nppp,maxits,telaps,nodepss,titers,tppp_s |
| 1228 1229C Solve libCEED algebraic RHS 1230 if (nid.eq.0) write (6,*) "libCEED algebraic RHS" 1231 maxit = 100 1232 tstart = dnekclock() 1233 call cggos(u1,r2,h1,h2,vmult,binvm1,tol,ceed,op_diffusion, 1234 $ vec_p1,vec_ap1,maxit,'bp3') 1235 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 1245 write(6,*) "UNCONVERGED CG" 1246 endif 1247 if (dabs(er1)>1e-9) then 1248 write(6,*) "ERROR IS TOO LARGE" 1249 endif 1250 endif 1251 1252 nx = nx1-1 | 1228 1229C Solve libCEED algebraic RHS 1230 if (nid.eq.0) write (6,*) "libCEED algebraic RHS" 1231 maxit = 100 1232 tstart = dnekclock() 1233 call cggos(u1,r2,h1,h2,vmult,binvm1,tol,ceed,op_diffusion, 1234 $ vec_p1,vec_ap1,maxit,'bp3') 1235 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 1245 write(6,*) "UNCONVERGED CG" 1246 endif 1247 if (dabs(er1)>1e-9) then 1248 write(6,*) "ERROR IS TOO LARGE" 1249 endif 1250 endif 1251 1252 nx = nx1-1 |
| 1253 ndof = nelgt ! ndofs 1254 ndof = ndof*(nx**ldim) ! ndofs 1255 nppp = ndof/np ! ndofs/proc | 1253 nnode = nelgt ! nnodes 1254 nnode = nnode*(nx**ldim) ! nnodes 1255 nppp = nnode/np ! nnodes/proc |
| 1256 | 1256 |
| 1257 dofpss = ndof/telaps ! DOF/sec - scalar form | 1257 nodepss = nnode/telaps ! DOF/sec - scalar form |
| 1258 titers = telaps/maxits ! time per iteration 1259 tppp_s = titers/nppp ! time per iteraton per local point 1260 1261 if (nid.eq.0) write(6,1) 'case scalar:' | 1258 titers = telaps/maxits ! time per iteration 1259 tppp_s = titers/nppp ! time per iteraton per local point 1260 1261 if (nid.eq.0) write(6,1) 'case scalar:' |
| 1262 $ ,np,nx,nelt,nelgt,ndof,nppp,maxits,telaps,dofpss,titers,tppp_s | 1262 $ ,np,nx,nelt,nelgt,nnode,nppp,maxits,telaps,nodepss,titers,tppp_s |
| 1263 1264C Solve Nek5000 algebraic RHS 1265 if (nid.eq.0) write (6,*) "Nek5000 algebraic RHS" 1266 maxit = 100 1267 tstart = dnekclock() 1268 call cggos(u1,r3,h1,h2,vmult,binvm1,tol,ceed,op_diffusion, 1269 $ vec_p1,vec_ap1,maxit,'bp3') 1270 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 1280 write(6,*) "UNCONVERGED CG" 1281 endif 1282 if (dabs(er1)>1e-9) then 1283 write(6,*) "ERROR IS TOO LARGE" 1284 endif 1285 endif 1286 1287 nx = nx1-1 | 1263 1264C Solve Nek5000 algebraic RHS 1265 if (nid.eq.0) write (6,*) "Nek5000 algebraic RHS" 1266 maxit = 100 1267 tstart = dnekclock() 1268 call cggos(u1,r3,h1,h2,vmult,binvm1,tol,ceed,op_diffusion, 1269 $ vec_p1,vec_ap1,maxit,'bp3') 1270 tstop = dnekclock() --- 9 unchanged lines hidden (view full) --- 1280 write(6,*) "UNCONVERGED CG" 1281 endif 1282 if (dabs(er1)>1e-9) then 1283 write(6,*) "ERROR IS TOO LARGE" 1284 endif 1285 endif 1286 1287 nx = nx1-1 |
| 1288 ndof = nelgt ! ndofs 1289 ndof = ndof*(nx**ldim) ! ndofs 1290 nppp = ndof/np ! ndofs/proc | 1288 nnode = nelgt ! nnodes 1289 nnode = nnode*(nx**ldim) ! nnodes 1290 nppp = nnode/np ! nnodes/proc |
| 1291 | 1291 |
| 1292 dofpss = ndof/telaps ! DOF/sec - scalar form | 1292 nodepss = nnode/telaps ! DOF/sec - scalar form |
| 1293 titers = telaps/maxits ! time per iteration 1294 tppp_s = titers/nppp ! time per iteraton per local point 1295 1296 if (nid.eq.0) write(6,1) 'case scalar:' | 1293 titers = telaps/maxits ! time per iteration 1294 tppp_s = titers/nppp ! time per iteraton per local point 1295 1296 if (nid.eq.0) write(6,1) 'case scalar:' |
| 1297 $ ,np,nx,nelt,nelgt,ndof,nppp,maxits,telaps,dofpss,titers,tppp_s | 1297 $ ,np,nx,nelt,nelgt,nnode,nppp,maxits,telaps,nodepss,titers,tppp_s |
| 1298 1299 1 format(a12,i7,i3,i7,i10,i14,i10,i4,1p4e13.5) 1300 3 format(i3,i9,e12.4,1x,a8,i9) 1301 1302C Destroy ceed handles 1303 call ceedvectordestroy(vec_p1,err) 1304 call ceedvectordestroy(vec_ap1,err) 1305 call ceedvectordestroy(vec_rhs,err) --- 378 unchanged lines hidden --- | 1298 1299 1 format(a12,i7,i3,i7,i10,i14,i10,i4,1p4e13.5) 1300 3 format(i3,i9,e12.4,1x,a8,i9) 1301 1302C Destroy ceed handles 1303 call ceedvectordestroy(vec_p1,err) 1304 call ceedvectordestroy(vec_ap1,err) 1305 call ceedvectordestroy(vec_rhs,err) --- 378 unchanged lines hidden --- |