Commit a4cc5f19 authored by Anne Reinarz's avatar Anne Reinarz

fix for 2d vectorised pde

parent e1e8d2d4
......@@ -79,8 +79,10 @@ for (int xyz = 0; xyz < {{(nDof**nDim)-((nDof**nDim)%vectSize)}}; xyz+={{vectSiz
{% if useNCP or useFusedSource %}
AoS2SoA<{{nVar}}, {{nVarPad}}>(gradQ+{{nVarPad}}*xyz, gradQt_x_block); //x
AoS2SoA<{{nVar}}, {{nVarPad}}>(gradQ+{{nVarPad}}*xyz+{{1*nVarPad*(nDof**nDim)}}, gradQt_y_block); //y
{% if nDim == 3 %}
AoS2SoA<{{nVar}}, {{nVarPad}}>(gradQ+{{nVarPad}}*xyz+{{2*nVarPad*(nDof**nDim)}}, gradQt_z_block); //z
{% endif %}
{% endif %}
{% if useFusedSource %}
// FusedSource
#ifdef USE_IPO
......@@ -117,8 +119,10 @@ for (int xyz = 0; xyz < {{(nDof**nDim)-((nDof**nDim)%vectSize)}}; xyz+={{vectSiz
{% if useNCP or useFusedSource %}
AoS2SoA_rest<{{nVar}}, {{nVarPad}}, {{rest_size}}>(gradQ+{{nVarPad*rest_start}}, gradQt_x_block); //x
AoS2SoA_rest<{{nVar}}, {{nVarPad}}, {{rest_size}}>(gradQ+{{nVarPad*rest_start+1*nVarPad*(nDof**nDim)}}, gradQt_y_block); //y
{% if nDim == 3 %}
AoS2SoA_rest<{{nVar}}, {{nVarPad}}, {{rest_size}}>(gradQ+{{nVarPad*rest_start+2*nVarPad*(nDof**nDim)}}, gradQt_z_block); //z
{% endif %}
{% endif %}
{% if useFusedSource %}
// FusedSource
#ifdef USE_IPO
......@@ -212,4 +216,4 @@ for(int xyz = 0; xyz < {{nDof**nDim}}; xyz++) {
}
}
{% endif %}{# non vectorized PDE block#}
{% endwith %}
\ No newline at end of file
{% endwith %}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment