Problemy z dodawaniem rozszerzenia do typu DenseVector

Question

Próbuję dodać metodę do klasy DenseVector, aby móc odróżnić wektor n razy. Poniższa nie wydają się działać jako rodzaj wnioskowania skarży się, że typ Vector nie jest zgodny z typem DenseVector:

open System 
open System.IO 
open Deedle 
open MathNet.Numerics 
open MathNet.Numerics.LinearAlgebra 
open MathNet.Numerics.LinearAlgebra.Double 
open MathNet.Numerics.Distributions 

[<Extension>] 
type DenseVector with 
    member this.diffVector (v : DenseVector) (n : int) = 
     let rec run (v : DenseVector) (n : int) = 
      match n with 
      | 0 -> v 
      | _ -> run (v.[ 1 .. v.Count-1 ] - v.[ 0 .. (v.Count-1)-1 ]) (n - 1) 
     run v n 

v.[ 0 .. (v.Count-1)-1 ] w powyższym jest przyczyną problemów. Dlaczego jest on wywiedziony jako Vector, a nie DenseVector, co jest przekazywane do funkcji? Czy mój sposób dodania metody rozszerzenia jest poprawny?

Answer 1

Przyczyną problemu jest fakt, że extension method który definiuje metodę GetSlice dla wszystkich Vector<'T> zwraca nie Vector<'T>DenseVector. Więc jeśli używasz cięcia, które działa przez wywołanie GetSlice, otrzymasz Vector, a powyższy kod nie będzie działał zgodnie z oczekiwaniami.

nie jestem pewien, w jaki sposób wewnętrzne pracy MathNet, ale może po prostu być w stanie dokonać rozszerzenia metody pracy dla wszystkich Vector<'T> sz coś takiego:

type Vector with 
    member this.diffVector (v : Vector<'T>) (n : int) = 
     let rec run (v : Vector<'T>) (n : int) = 
      match n with 
      | 0 -> v 
      | _ -> run (v.[ 1 .. v.Count-1 ] - v.[ 0 .. (v.Count-1)-1 ]) (n - 1) 
     run v n 

Alternatywnie, może się okazać, że możesz po prostu spuścić z Vector do DenseVector bezpiecznie przez coś takiego. To może, ale nie musi działać:

[<Extension>] 
type DenseVector with 
    member this.diffVector (v : DenseVector) (n : int) = 
     let rec run (v : DenseVector) (n : int) = 
      match n with 
      | 0 -> v 
      | _ -> let v1 = v.[ 1 .. v.Count-1 ] :?> DenseVector 
        let v2 = v.[ 0 .. (v.Count-1)-1 ] :?> DenseVector 
        run (v1 - v2) (n - 1) 
     run v n