[deal.II] Bug in Tensor::operator*()

[Alexander Grayver]

Hello,

There is typo in this routine:

template <int dim, typename Number>
inline
Number
Tensor<1,dim,Number>::operator * (const Tensor<1,dim,Number> &p) const
{
...
    default:
      Number q = values[0] * values[0];
      for (unsigned int i=1; i<dim; ++i)
        q += values[i] * p.values[i];
      return q;
    }
}

I guess it should be:

Number q = values[0] * p.values[0];

--
Regards,
Alexander

_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii