if(epsilon2d.rank() != denom.rank() && denom.length() == epsilon2d.length()) {
                    epsilon2d = epsilon2d.reshape(denom.shape());
                }