anchor.diff:              else failwithf "Cannot add Tensors with shapes %A, %A" a.Shape b.Shape                
           // TODO: implement general broadcasting additions
           else failwithf "Cannot add Tensors with shapes %A, %A" a.Shape b.Shape
  -    static member (+) (a:Tensor, b) = a + a.Create(b)
  -    static member (+) (a, b:Tensor) = b.Create(a) + b
  +    static member (+) (a:Tensor, b) = a + a.CreateLike(b)
  +    static member (+) (a, b:Tensor) = b.CreateLike(a) + b
       member t1.Add(t2:Tensor) = t1 + t2
  -    member t1.Add(t2) = t1 + t1.Create(t2)
  +    member t1.Add(t2) = t1 + t1.CreateLike(t2)
   
       static member (-) (a:Tensor, b:Tensor) =
           if a.Shape = b.Shape then
  
positive.code:             let inline dfTensorRevCT(a,b) = SubTConstT0(b)
              Tensor.OpBinary(a, b, fRaw, fTensor, dfTensorFwdTT, dfTensorFwdTC, dfTensorFwdCT, dfTensorRevTT, dfTensorRevTC, dfTensorRevCT)
          else failwithf "Cannot subtract Tensors with shapes %A, %A" a.Shape b.Shape
      static member (-) (a:Tensor, b) = a - a.Create(b)
      static member (-) (a, b:Tensor) = b.Create(a) - b
      member t1.Sub(t2:Tensor) = t1 - t2
      member t1.Sub(t2) = t1 - t1.Create(t2)
  
      static member (*) (a:Tensor, b:Tensor) =
          if a.Shape = b.Shape then
  
negative.code:             Tensor.OpBinary(a, b, fRaw, fTensor, dfTensorFwdTT, dfTensorFwdTC, dfTensorFwdCT, dfTensorRevTT, dfTensorRevTC, dfTensorRevCT)
          // TODO: implement general broadcasting?
          else failwithf "Cannot add Tensors with shapes %A, %A" a.Shape b.Shape
      static member (*) (a:Tensor, b) = a * a.CreateLike(b)
      static member (*) (a, b:Tensor) = b.CreateLike(a) * b
      member t1.Mul(t2:Tensor) = t1 * t2
      member t1.Mul(t2) = t1 * t1.CreateLike(t2)
  
      static member (/) (a:Tensor, b:Tensor) =
          if a.Shape = b.Shape then
  