Method CreateSparseMatrix

linalg/special_matrix.cpp:112–134

  }

  shared_ptr<BaseSparseMatrix> Embedding :: CreateSparseMatrix() const
  {
    Array<int> ai(range.Size()), aj(range.Size());

    for (int i = 0; i < range.Size(); i++)
      {
        ai[i] = range.First()+i;
        aj[i] = i;
      }

    if (IsComplex())
      {
        Array<Complex> vals(range.Size());
        vals = 1.;
        return SparseMatrix<Complex>::CreateFromCOO (ai, aj, vals, Height(), Width());
      }
    else
      {
        Array<double> vals(range.Size());
        vals = 1.;
        return SparseMatrix<double>::CreateFromCOO (ai, aj, vals, Height(), Width());
      }
  }