我有两个向量,如下所示:

vdA = { 8.0, 7.0, 6.0 }
vdB = { 0.0, 1.0, 2.0, 3.0 }


我基本上想要一个向量vdX,其结果是将vdA的所有元素与vdB的所有值相加。

vdX = {
        8.0, 9.0, 10.0 11.0,
        7.0, 8.0, 9.0, 10.0,
        6.0, 7.0, 8.0, 9.0
      }


使用MathNet.Numerics,我找不到执行此操作的函数。

在C#中,我编写此代码来执行此操作

Vector<double> vdA = new DenseVector(new[] { 8.0, 7.0, 6.0 });
Vector<double> vdB = new DenseVector(new[] { 0.0, 1.0, 2.0, 3.0 });

List<double> resultSumVector = new List<double>();
foreach (double vectorValueA in vdA.Enumerate())
   foreach (double vectorValueB in vdB.Enumerate())
      resultSumVector.Add(vectorValueA + vectorValueB);
Vector<double> vdX = new DenseVector(resultSumVector.ToArray());


还有其他选择可以使用c#中的Math.Net Numerics更快地完成此任务吗?

最佳答案

您基本上需要一个cross join in Linq。您可以编写一个扩展方法,这样看起来就像是Math.Net方法:

namespace MathNet.Numerics
{
    public static class DenseVectorExtensions
    {
        public static DenseVector AddAlls(this DenseVector vdA, DenseVector vdB)
        {
           return DenseVector.OfEnumerable(
                     vdA.SelectMany(x => vdB, (y, z) => { return y + z; })
                  );
        }
    }
}


用法:

var vdA = new DenseVector(new[] { 8.0, 7.0, 6.0 });
var vdB = new DenseVector(new[] { 0.0, 1.0, 2.0, 3.0 });
var vdX = vdA.AddAlls(vdB);


这并不是特别快。

10-08 05:59