1!----------------------------------------------------------------------- 2 subroutine setup(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,& 3& u15,u16,v1,v2,v3,v4,v5,v6,v7,v8,v9,v10,v11,v12,v13,v14,v15,v16,ierr) 4 real*8 ctx 5 real*8 u1(1) 6 real*8 u2(1) 7 real*8 v1(1) 8 integer q,ierr 9 10 do i=1,q 11 v1(i)=u1(i)*u2(i) 12 enddo 13 14 ierr=0 15 end 16!----------------------------------------------------------------------- 17 subroutine mass(ctx,q,u1,u2,u3,u4,u5,u6,u7,u8,u9,u10,u11,u12,u13,u14,& 18& u15,u16,v1,v2,v3,v4,v5,v6,v7,v8,v9,v10,v11,v12,v13,v14,v15,v16,ierr) 19 real*8 ctx 20 real*8 u1(1) 21 real*8 u2(1) 22 real*8 v1(1) 23 integer q,ierr 24 25 do i=1,q 26 v1(i)=u2(i)*u1(i) 27 enddo 28 29 ierr=0 30 end 31!----------------------------------------------------------------------- 32 program test 33 34 include 'ceedf.h' 35 36 integer ceed,err,i,j 37 integer erestrictx,erestrictu,erestrictxi,erestrictui 38 integer bx,bu 39 integer qf_setup,qf_mass 40 integer op_setup,op_mass 41 integer qdata,x,u,v 42 integer nelem,p,q 43 parameter(nelem=15) 44 parameter(p=5) 45 parameter(q=8) 46 integer nx,nu 47 parameter(nx=nelem+1) 48 parameter(nu=nelem*(p-1)+1) 49 integer indx(nelem*2) 50 integer indu(nelem*p) 51 real*8 arrx(nx) 52 integer*8 voffset,xoffset 53 54 real*8 hv(nu) 55 real*8 total 56 57 character arg*32 58 59 external setup,mass 60 61 call getarg(1,arg) 62 call ceedinit(trim(arg)//char(0),ceed,err) 63 64 call ceedvectorcreate(ceed,nu,u,err) 65 call ceedvectorsetvalue(u,1.d0,err) 66 67 call ceedvectorcreate(ceed,nu,v,err) 68 69 do i=0,nx-1 70 arrx(i+1)=i/(nx-1.d0) 71 enddo 72 call ceedvectorcreate(ceed,nx,x,err) 73 xoffset=0 74 call ceedvectorsetarray(x,ceed_mem_host,ceed_use_pointer,arrx,xoffset,err) 75 76 call ceedvectorcreate(ceed,nelem*q,qdata,err) 77 78 do i=0,nelem-1 79 indx(2*i+1)=i 80 indx(2*i+2)=i+1 81 enddo 82 call ceedelemrestrictioncreate(ceed,nelem,2,nx,1,ceed_mem_host,& 83 & ceed_use_pointer,indx,erestrictx,err) 84 call ceedelemrestrictioncreateidentity(ceed,nelem,2,2*nelem,1,& 85 & erestrictxi,err) 86 87 do i=0,nelem-1 88 do j=0,p-1 89 indu(p*i+j+1)=i*(p-1)+j 90 enddo 91 enddo 92 93 call ceedelemrestrictioncreate(ceed,nelem,p,nu,1,ceed_mem_host,& 94 & ceed_use_pointer,indu,erestrictu,err) 95 call ceedelemrestrictioncreateidentity(ceed,nelem,q,q*nelem,1,& 96 & erestrictui,err) 97 98 call ceedbasiscreatetensorh1lagrange(ceed,1,1,2,q,ceed_gauss,bx,err) 99 call ceedbasiscreatetensorh1lagrange(ceed,1,1,p,q,ceed_gauss,bu,err) 100 101 call ceedqfunctioncreateinterior(ceed,1,setup,& 102 &SOURCE_DIR& 103 &//'t500-operator.h:setup'//char(0),qf_setup,err) 104 call ceedqfunctionaddinput(qf_setup,'_weight',1,ceed_eval_weight,err) 105 call ceedqfunctionaddinput(qf_setup,'dx',1,ceed_eval_grad,err) 106 call ceedqfunctionaddoutput(qf_setup,'rho',1,ceed_eval_none,err) 107 108 call ceedqfunctioncreateinterior(ceed,1,mass,& 109 &SOURCE_DIR& 110 &//'t500-operator.h:mass'//char(0),qf_mass,err) 111 call ceedqfunctionaddinput(qf_mass,'rho',1,ceed_eval_none,err) 112 call ceedqfunctionaddinput(qf_mass,'u',1,ceed_eval_interp,err) 113 call ceedqfunctionaddoutput(qf_mass,'v',1,ceed_eval_interp,err) 114 115 call ceedoperatorcreate(ceed,qf_setup,ceed_null,ceed_null,op_setup,err) 116 call ceedoperatorcreate(ceed,qf_mass,ceed_null,ceed_null,op_mass,err) 117 118 call ceedoperatorsetfield(op_setup,'_weight',erestrictxi,& 119 & ceed_notranspose,bx,ceed_vector_none,err) 120 call ceedoperatorsetfield(op_setup,'dx',erestrictx,& 121 & ceed_notranspose,bx,x,err) 122 call ceedoperatorsetfield(op_setup,'rho',erestrictui,& 123 & ceed_notranspose,ceed_basis_collocated,qdata,err) 124 call ceedoperatorsetfield(op_mass,'rho',erestrictui,& 125 & ceed_notranspose,ceed_basis_collocated,qdata,err) 126 call ceedoperatorsetfield(op_mass,'u',erestrictu,& 127 & ceed_notranspose,bu,u,err) 128 call ceedoperatorsetfield(op_mass,'v',erestrictu,& 129 & ceed_notranspose,bu,v,err) 130 131 call ceedoperatorapply(op_setup,ceed_null,ceed_null,& 132 & ceed_request_immediate,err) 133 call ceedoperatorapply(op_mass,ceed_null,ceed_null,& 134 & ceed_request_immediate,err) 135 136 call ceedvectorgetarrayread(v,ceed_mem_host,hv,voffset,err) 137 total=0. 138 do i=1,nu 139 total=total+hv(voffset+i) 140 enddo 141 if (abs(total-1.)>1.0d-10) then 142! LCOV_EXCL_START 143 write(*,*) 'Computed Area: ',total,' != True Area: 1.0' 144! LCOV_EXCL_STOP 145 endif 146 call ceedvectorrestorearrayread(v,hv,voffset,err) 147 148 call ceedvectordestroy(qdata,err) 149 call ceedvectordestroy(x,err) 150 call ceedvectordestroy(u,err) 151 call ceedvectordestroy(v,err) 152 call ceedoperatordestroy(op_mass,err) 153 call ceedoperatordestroy(op_setup,err) 154 call ceedqfunctiondestroy(qf_mass,err) 155 call ceedqfunctiondestroy(qf_setup,err) 156 call ceedbasisdestroy(bu,err) 157 call ceedbasisdestroy(bx,err) 158 call ceedelemrestrictiondestroy(erestrictu,err) 159 call ceedelemrestrictiondestroy(erestrictx,err) 160 call ceedelemrestrictiondestroy(erestrictui,err) 161 call ceedelemrestrictiondestroy(erestrictxi,err) 162 call ceeddestroy(ceed,err) 163 end 164!----------------------------------------------------------------------- 165