Skip to content
Snippets Groups Projects
Commit 59c92154 authored by Stephan Hilb's avatar Stephan Hilb
Browse files

fix performance for image projection

parent e398949b
No related branches found
No related tags found
No related merge requests found
......@@ -177,8 +177,9 @@ function project_img(space::FeSpace, img)
for jdim in 1:nrdims, ldofj in 1:nldofs
gdofj = space.dofmap[jdim, ldofj, cell]
phij = SArray{Tuple{space.size...}}(qphi[:, jdim, ldofj, k])
dphij = SArray{Tuple{space.size..., d}}(dqphi[:, :, jdim, ldofj, k] * delmapinv)
phij = SArray{Tuple{space.size...}}(view(qphi, :, jdim, ldofj, k))
dphij = SArray{Tuple{space.size..., d}}(
SArray{Tuple{nrdims, d}}(view(dqphi, :, :, jdim, ldofj, k)) * delmapinv)
lv = qw[k] * l(x, phij, dphij; opvalues...) * intel
b[gdofj] += lv
......@@ -187,8 +188,9 @@ function project_img(space::FeSpace, img)
for idim in 1:nrdims, ldofi in 1:nldofs
gdofi = space.dofmap[idim, ldofi, cell]
phii = SArray{Tuple{space.size...}}(qphi[:, idim, ldofi, k])
dphii = SArray{Tuple{space.size..., d}}(dqphi[:, :, idim, ldofi, k] * delmapinv)
phii = SArray{Tuple{space.size...}}(view(qphi, :, idim, ldofi, k))
dphii = SArray{Tuple{space.size..., d}}(
SArray{Tuple{nrdims, d}}(view(dqphi, :, :, idim, ldofi, k)) * delmapinv)
av = qw[k] * a(x, phii, dphii, phij, dphij; opvalues...) * intel
push!(I, gdofi)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment