return localInnerProduct(s1.elem(), s2.elem());
typename BinaryReturn<ILattice<T1,N>, ILattice<T2,N>, FnLocalInnerProduct>::Type_t  d;  for(int i=0; i < N; ++i) d.elem(i) = localInnerProduct(s1.elem(i), s2.elem(i));  return d;
localInnerProduct(const IScalar<T1>& s1, const IScalar<T2>& s2){  return localInnerProduct(s1.elem(), s2.elem());}