linear_algebra_internal.h

Go to the documentation of this file.

#ifndef CMR_LINALG_INTERNAL_H
#define CMR_LINALG_INTERNAL_H
 
#include <cmr/linear_algebra.h>
 
#ifdef __cplusplus
extern "C" {
#endif
 
static inline
int moduloNonnegative(int p, int q)
{
  if (q == 0)
    return p;
  if (q < 0)
    q = -q;
  int r = p % q;
  if (r < 0)
    r += q;
  return r;
}
 
static inline
int moduloTernary(int p, int q)
{
  if (q == 0)
    return p;
  if (q < 0)
    q = -q;
  int r = p % q;
  if (r < 0)
    r += q;
  return (r == 2 && q == 3) ? -1 : r;
}
 
CMR_ERROR CMRintmatComputeUpperDiagonal(
  CMR* cmr,                   
  CMR_INTMAT* matrix,         
  bool invert,                
  size_t* prank,              
  CMR_SUBMAT** ppermutations, 
  CMR_INTMAT** presult,       
  CMR_INTMAT** ptranspose     
);
 
#ifdef __cplusplus
}
#endif
 
#endif /* CMR_LINALG_INTERNAL_H */