yade-dev team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

This was not needed; makeTensorProduct is in the part of
lib/base/Math.hpp that will be removed, as in Eigen you just write
stressTensor=f*branch.transpose() for outer product. If you put back
makeTensorProduct, I will (once we drop wm3) grep the code for such
functions and update it.

> -		// tensorial product f*branch
> -		stressTensor+=makeTensorProduct( f, branch );
> +		// tensorial product f*branch (hand-write the tensor product to prevent matrix instanciation inside the loop by makeTensorProduct)
> +		stressTensor(0,0)+=f(0)*branch(0); stressTensor(1,0)+=f(1)*branch(0); stressTensor(2,0)+=f(2)*branch(0);
> +		stressTensor(0,1)+=f(0)*branch(1); stressTensor(1,1)+=f(1)*branch(1); stressTensor(2,1)+=f(2)*branch(1);
> +		stressTensor(0,2)+=f(0)*branch(2); stressTensor(1,2)+=f(1)*branch(2); stressTensor(2,2)+=f(2)*branch(2);
> +		//stressTensor+=makeTensorProduct( f, branch );