internal override void MapIndexedToUnchecked <TU>(VectorStorage <TU> target, Func <int, T, TU> f, Zeros zeros, ExistingData existingData)
{
if (target is DenseVectorStorage <TU> denseTarget)
{
CommonParallel.For(0, Data.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
denseTarget.Data[i] = f(i, Data[i]);
}
});
return;
}
// FALL BACK
for (int i = 0; i < Length; i++)
{
target.At(i, f(i, Data[i]));
}
}
// VECTOR COPY
internal override void CopyToUnchecked(VectorStorage <T> target, ExistingData existingData)
{
if (target is DenseVectorStorage <T> denseTarget)
{
if (!ReferenceEquals(this, denseTarget))
{
Array.Copy(Data, 0, denseTarget.Data, 0, Data.Length);
}
return;
}
if (target is SparseVectorStorage <T> sparseTarget)
{
var indices = new List <int>();
var values = new List <T>();
for (int i = 0; i < Data.Length; i++)
{
var item = Data[i];
if (!Zero.Equals(item))
{
values.Add(item);
indices.Add(i);
}
}
sparseTarget.Indices = indices.ToArray();
sparseTarget.Values = values.ToArray();
sparseTarget.ValueCount = values.Count;
return;
}
// FALL BACK
for (int i = 0; i < Data.Length; i++)
{
target.At(i, Data[i]);
}
}
internal virtual void CopySubVectorToUnchecked(VectorStorage <T> target,
int sourceIndex, int targetIndex, int count, ExistingData existingData)
{
if (ReferenceEquals(this, target))
{
var tmp = new T[count];
for (int i = 0; i < tmp.Length; i++)
{
tmp[i] = At(i + sourceIndex);
}
for (int i = 0; i < tmp.Length; i++)
{
At(i + targetIndex, tmp[i]);
}
return;
}
for (int i = sourceIndex, ii = targetIndex; i < sourceIndex + count; i++, ii++)
{
target.At(ii, At(i));
}
}
// VECTOR COPY
internal override void CopyToUnchecked(VectorStorage <T> target, ExistingData existingData)
{
if (target is SparseVectorStorage <T> sparseTarget)
{
CopyToUnchecked(sparseTarget);
return;
}
// FALL BACK
if (existingData == ExistingData.Clear)
{
target.Clear();
}
if (ValueCount != 0)
{
for (int i = 0; i < ValueCount; i++)
{
target.At(Indices[i], Values[i]);
}
}
}