bps.usr (288c044332e33f37503f09b6484fec9d0a55fba1) bps.usr (a2fa791084b022ffd3d077ae9be76b4b170bb84e)
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.

--- 804 unchanged lines hidden (view full) ---

813 $ __FILE__
814 $ //':masssetupf',qf_setup,err)
815 call ceedqfunctionaddinput(qf_setup,'x',ncompx,
816 $ ceed_eval_interp,err)
817 call ceedqfunctionaddinput(qf_setup,'dx',ncompx*ldim,
818 $ ceed_eval_grad,err)
819 call ceedqfunctionaddinput(qf_setup,'weight',ncompu,
820 $ ceed_eval_weight,err)
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.

--- 804 unchanged lines hidden (view full) ---

813 $ __FILE__
814 $ //':masssetupf',qf_setup,err)
815 call ceedqfunctionaddinput(qf_setup,'x',ncompx,
816 $ ceed_eval_interp,err)
817 call ceedqfunctionaddinput(qf_setup,'dx',ncompx*ldim,
818 $ ceed_eval_grad,err)
819 call ceedqfunctionaddinput(qf_setup,'weight',ncompu,
820 $ ceed_eval_weight,err)
821 call ceedqfunctionaddoutput(qf_setup,'rho',ncompu,
821 call ceedqfunctionaddoutput(qf_setup,'qdata',ncompu,
822 $ ceed_eval_none,err)
823 call ceedqfunctionaddoutput(qf_setup,'rhs',ncompu,
824 $ ceed_eval_interp,err)
825
826 call ceedqfunctioncreateinterior(ceed,1,massf,
827 $ __FILE__
828 $ //':massf',qf_mass,err)
829 call ceedqfunctionaddinput(qf_mass,'u',ncompu,
830 $ ceed_eval_interp,err)
822 $ ceed_eval_none,err)
823 call ceedqfunctionaddoutput(qf_setup,'rhs',ncompu,
824 $ ceed_eval_interp,err)
825
826 call ceedqfunctioncreateinterior(ceed,1,massf,
827 $ __FILE__
828 $ //':massf',qf_mass,err)
829 call ceedqfunctionaddinput(qf_mass,'u',ncompu,
830 $ ceed_eval_interp,err)
831 call ceedqfunctionaddinput(qf_mass,'rho',ncompu,
831 call ceedqfunctionaddinput(qf_mass,'qdata',ncompu,
832 $ ceed_eval_none,err)
833 call ceedqfunctionaddoutput(qf_mass,'v',ncompu,
834 $ ceed_eval_interp,err)
835
836C Create ceed operators
837 call ceedoperatorcreate(ceed,qf_setup,
838 $ ceed_null,ceed_null,op_setup,err)
839 call ceedoperatorsetfield(op_setup,'x',erstrctx,
840 $ ceed_notranspose,basisx,ceed_vector_active,err)
841 call ceedoperatorsetfield(op_setup,'dx',erstrctx,
842 $ ceed_notranspose,basisx,ceed_vector_active,err)
843 call ceedoperatorsetfield(op_setup,'weight',erstrctx,
844 $ ceed_notranspose,basisx,ceed_vector_none,err)
832 $ ceed_eval_none,err)
833 call ceedqfunctionaddoutput(qf_mass,'v',ncompu,
834 $ ceed_eval_interp,err)
835
836C Create ceed operators
837 call ceedoperatorcreate(ceed,qf_setup,
838 $ ceed_null,ceed_null,op_setup,err)
839 call ceedoperatorsetfield(op_setup,'x',erstrctx,
840 $ ceed_notranspose,basisx,ceed_vector_active,err)
841 call ceedoperatorsetfield(op_setup,'dx',erstrctx,
842 $ ceed_notranspose,basisx,ceed_vector_active,err)
843 call ceedoperatorsetfield(op_setup,'weight',erstrctx,
844 $ ceed_notranspose,basisx,ceed_vector_none,err)
845 call ceedoperatorsetfield(op_setup,'rho',erstrctw,
845 call ceedoperatorsetfield(op_setup,'qdata',erstrctw,
846 $ ceed_notranspose,ceed_basis_collocated,
847 $ ceed_vector_active,err)
848 call ceedoperatorsetfield(op_setup,'rhs',erstrctu,
849 $ ceed_notranspose,basisu,vec_rhs,err)
850
851 call ceedoperatorcreate(ceed,qf_mass,
852 $ ceed_null,ceed_null,op_mass,err)
853 call ceedoperatorsetfield(op_mass,'u',erstrctu,
854 $ ceed_notranspose,basisu,ceed_vector_active,err)
846 $ ceed_notranspose,ceed_basis_collocated,
847 $ ceed_vector_active,err)
848 call ceedoperatorsetfield(op_setup,'rhs',erstrctu,
849 $ ceed_notranspose,basisu,vec_rhs,err)
850
851 call ceedoperatorcreate(ceed,qf_mass,
852 $ ceed_null,ceed_null,op_mass,err)
853 call ceedoperatorsetfield(op_mass,'u',erstrctu,
854 $ ceed_notranspose,basisu,ceed_vector_active,err)
855 call ceedoperatorsetfield(op_mass,'rho',erstrctw,
855 call ceedoperatorsetfield(op_mass,'qdata',erstrctw,
856 $ ceed_notranspose,ceed_basis_collocated,
857 $ vec_qdata,err)
858 call ceedoperatorsetfield(op_mass,'v',erstrctu,
859 $ ceed_notranspose,basisu,ceed_vector_active,err)
860
861C Compute setup data
862 call ceedvectorsetarray(vec_rhs,ceed_mem_host,
863 $ ceed_use_pointer,r1,offset,err)

--- 266 unchanged lines hidden (view full) ---

1130 $ __FILE__
1131 $ //':diffsetupf'//char(0),qf_setup,err)
1132 call ceedqfunctionaddinput(qf_setup,'x',ncompx,
1133 $ ceed_eval_interp,err)
1134 call ceedqfunctionaddinput(qf_setup,'dx',ncompx*ldim,
1135 $ ceed_eval_grad,err)
1136 call ceedqfunctionaddinput(qf_setup,'weight',ncompu,
1137 $ ceed_eval_weight,err)
856 $ ceed_notranspose,ceed_basis_collocated,
857 $ vec_qdata,err)
858 call ceedoperatorsetfield(op_mass,'v',erstrctu,
859 $ ceed_notranspose,basisu,ceed_vector_active,err)
860
861C Compute setup data
862 call ceedvectorsetarray(vec_rhs,ceed_mem_host,
863 $ ceed_use_pointer,r1,offset,err)

--- 266 unchanged lines hidden (view full) ---

1130 $ __FILE__
1131 $ //':diffsetupf'//char(0),qf_setup,err)
1132 call ceedqfunctionaddinput(qf_setup,'x',ncompx,
1133 $ ceed_eval_interp,err)
1134 call ceedqfunctionaddinput(qf_setup,'dx',ncompx*ldim,
1135 $ ceed_eval_grad,err)
1136 call ceedqfunctionaddinput(qf_setup,'weight',ncompu,
1137 $ ceed_eval_weight,err)
1138 call ceedqfunctionaddoutput(qf_setup,'rho',ngeo,
1138 call ceedqfunctionaddoutput(qf_setup,'qdata',ngeo,
1139 $ ceed_eval_none,err)
1140 call ceedqfunctionaddoutput(qf_setup,'rhs',ncompu,
1141 $ ceed_eval_interp,err)
1142
1143 call ceedqfunctioncreateinterior(ceed,1,diffusionf,
1144 $ __FILE__
1145 $ //':diffusionf'//char(0),qf_diffusion,err)
1146 call ceedqfunctionaddinput(qf_diffusion,'u',ncompu*ldim,
1147 $ ceed_eval_grad,err)
1139 $ ceed_eval_none,err)
1140 call ceedqfunctionaddoutput(qf_setup,'rhs',ncompu,
1141 $ ceed_eval_interp,err)
1142
1143 call ceedqfunctioncreateinterior(ceed,1,diffusionf,
1144 $ __FILE__
1145 $ //':diffusionf'//char(0),qf_diffusion,err)
1146 call ceedqfunctionaddinput(qf_diffusion,'u',ncompu*ldim,
1147 $ ceed_eval_grad,err)
1148 call ceedqfunctionaddinput(qf_diffusion,'rho',ngeo,
1148 call ceedqfunctionaddinput(qf_diffusion,'qdata',ngeo,
1149 $ ceed_eval_none,err)
1150 call ceedqfunctionaddoutput(qf_diffusion,'v',ncompu*ldim,
1151 $ ceed_eval_grad,err)
1152
1153C Create ceed operators
1154 call ceedoperatorcreate(ceed,qf_setup,
1155 $ ceed_null,ceed_null,op_setup,err)
1156 call ceedoperatorsetfield(op_setup,'x',erstrctx,
1157 $ ceed_notranspose,basisx,ceed_vector_active,err)
1158 call ceedoperatorsetfield(op_setup,'dx',erstrctx,
1159 $ ceed_notranspose,basisx,ceed_vector_active,err)
1160 call ceedoperatorsetfield(op_setup,'weight',erstrctx,
1161 $ ceed_notranspose,basisx,ceed_vector_none,err)
1149 $ ceed_eval_none,err)
1150 call ceedqfunctionaddoutput(qf_diffusion,'v',ncompu*ldim,
1151 $ ceed_eval_grad,err)
1152
1153C Create ceed operators
1154 call ceedoperatorcreate(ceed,qf_setup,
1155 $ ceed_null,ceed_null,op_setup,err)
1156 call ceedoperatorsetfield(op_setup,'x',erstrctx,
1157 $ ceed_notranspose,basisx,ceed_vector_active,err)
1158 call ceedoperatorsetfield(op_setup,'dx',erstrctx,
1159 $ ceed_notranspose,basisx,ceed_vector_active,err)
1160 call ceedoperatorsetfield(op_setup,'weight',erstrctx,
1161 $ ceed_notranspose,basisx,ceed_vector_none,err)
1162 call ceedoperatorsetfield(op_setup,'rho',erstrctw,
1162 call ceedoperatorsetfield(op_setup,'qdata',erstrctw,
1163 $ ceed_notranspose,ceed_basis_collocated,
1164 $ ceed_vector_active,err)
1165 call ceedoperatorsetfield(op_setup,'rhs',erstrctu,
1166 $ ceed_notranspose,basisu,vec_rhs,err)
1167
1168 call ceedoperatorcreate(ceed,qf_diffusion,
1169 $ ceed_null,ceed_null,op_diffusion,err)
1170 call ceedoperatorsetfield(op_diffusion,'u',erstrctu,
1171 $ ceed_notranspose,basisu,ceed_vector_active,err)
1163 $ ceed_notranspose,ceed_basis_collocated,
1164 $ ceed_vector_active,err)
1165 call ceedoperatorsetfield(op_setup,'rhs',erstrctu,
1166 $ ceed_notranspose,basisu,vec_rhs,err)
1167
1168 call ceedoperatorcreate(ceed,qf_diffusion,
1169 $ ceed_null,ceed_null,op_diffusion,err)
1170 call ceedoperatorsetfield(op_diffusion,'u',erstrctu,
1171 $ ceed_notranspose,basisu,ceed_vector_active,err)
1172 call ceedoperatorsetfield(op_diffusion,'rho',erstrctw,
1172 call ceedoperatorsetfield(op_diffusion,'qdata',erstrctw,
1173 $ ceed_notranspose,ceed_basis_collocated,
1174 $ vec_qdata,err)
1175 call ceedoperatorsetfield(op_diffusion,'v',erstrctu,
1176 $ ceed_notranspose,basisu,ceed_vector_active,err)
1177
1178C Compute setup data
1179 call ceedvectorsetarray(vec_rhs,ceed_mem_host,
1180 $ ceed_use_pointer,r1,offset,err)

--- 513 unchanged lines hidden --- 		1173 $ ceed_notranspose,ceed_basis_collocated,
1174 $ vec_qdata,err)
1175 call ceedoperatorsetfield(op_diffusion,'v',erstrctu,
1176 $ ceed_notranspose,basisu,ceed_vector_active,err)
1177
1178C Compute setup data
1179 call ceedvectorsetarray(vec_rhs,ceed_mem_host,
1180 $ ceed_use_pointer,r1,offset,err)

--- 513 unchanged lines hidden ---