Example #1
0
bool CLinkMatrix::rightMultiply(const C_FLOAT64 & alpha,
                                const CMatrix< C_FLOAT64> & m,
                                CMatrix< C_FLOAT64> & p) const
{
  bool success = true;

  if (m.numCols() != mRowPivots.size())
    {
      return false;
    }

  // p := alpha * (m1, m2) * (I, L0')' = alpha * m1 + alpha * m2 * L0

  p.resize(m.numRows(), getNumIndependent());

  char T = 'N';
  C_INT M = (C_INT) p.numCols();
  C_INT N = (C_INT) p.numRows();
  C_INT K = (C_INT) numRows();

  C_INT LDA = (C_INT) numCols();
  C_INT LDB = (C_INT) m.numCols();
  C_INT LDC = (C_INT) p.numCols();

  // p := m1

  C_FLOAT64 *pRowP = p.array();
  const C_FLOAT64 *pRowM = m.array();
  C_FLOAT64 *pEnd = p.array() + p.size();

  for (; pRowP < pEnd; pRowP += LDC, pRowM += LDB)
    {
      memcpy(pRowP, pRowM, sizeof(C_FLOAT64) * LDC);
    }

  // DGEMM (TRANSA, TRANSB, M, N, K, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
  // C := alpha A B + beta C

  dgemm_(&T, &T, &M, &N, &K,
         const_cast< C_FLOAT64 * >(&alpha),
         const_cast< C_FLOAT64 * >(array()), &LDA,
         const_cast< C_FLOAT64 * >(m.array()) + LDA, &LDB,
         const_cast< C_FLOAT64 * >(&alpha), p.array(), &LDC);

  return success;
}
Example #2
0
C_FLOAT64 * CState::resize(const unsigned C_INT32 & size)
{
  if (mSize != size)
    {
      C_FLOAT64 * pTmp = mpValues;

      if (size > 0)
        {
          mpValues = new C_FLOAT64[size];
          memcpy(mpValues, pTmp, sizeof(C_FLOAT64) * std::min(mSize, size));
        }
      else
        mpValues = NULL;

      pdeletev(pTmp);
      mSize = size;
    }

  updateIterator(getNumIndependent(), getNumDependent(), getNumFixed());

  return mpValues;
}