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 ); _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
