separation.h

Go to the documentation of this file.

#ifndef CMR_SEPARATION_H
#define CMR_SEPARATION_H
 
#include <cmr/env.h>
#include <cmr/matrix.h>
#include <cmr/element.h>
 
#include <assert.h>
#include <stdint.h>
 
#ifdef __cplusplus
extern "C" {
#endif
 
typedef enum
{
  CMR_SEPA_FIRST = 0,
  CMR_SEPA_SECOND = 1,
  CMR_SEPA_FLAG_RANK1 = 2,
  CMR_SEPA_FLAG_RANK2 = 4,
  CMR_SEPA_MASK_CHILD = 1,
  CMR_SEPA_MASK_EXTRA = CMR_SEPA_FLAG_RANK1 | CMR_SEPA_FLAG_RANK2
} CMR_SEPA_FLAGS;
 
typedef enum
{
  CMR_SEPA_TYPE_TWO = 2,
  CMR_SEPA_TYPE_THREE_DISTRIBUTED_RANKS = 3,
  CMR_SEPA_TYPE_THREE_CONCENTRATED_RANK = 4
} CMR_SEPA_TYPE;
 
typedef struct
{
  size_t numRows;     
  size_t numColumns;  
  int* rowsFlags;     
  int* columnsFlags;  
  CMR_SEPA_TYPE type; 
} CMR_SEPA;
 
CMR_EXPORT
CMR_ERROR CMRsepaCreate(
  CMR* cmr,           
  size_t numRows,     
  size_t numColumns,  
  CMR_SEPA** psepa    
);
 
CMR_EXPORT
CMR_ERROR CMRsepaFree(
  CMR* cmr,         
  CMR_SEPA** psepa  
);
 
CMR_EXPORT
CMR_ERROR CMRsepaTranspose(
  CMR* cmr,               
  CMR_SEPA* sepa,         
  CMR_SEPA** ptransposed  
);
 
CMR_EXPORT
CMR_ERROR CMRsepaComputeSizes(
  CMR_SEPA* sepa,                 
  size_t* pnumRowsTopLeft,        
  size_t* pnumColumnsTopLeft,     
  size_t* pnumRowsBottomRight,    
  size_t* pnumColumnsBottomRight  
);
 
CMR_EXPORT
CMR_ERROR CMRsepaFindBinaryRepresentatives(
  CMR* cmr,               
  CMR_SEPA* sepa,         
  CMR_CHRMAT* matrix,     
  CMR_CHRMAT* transpose,  
  bool* pswapped,         
  CMR_SUBMAT** pviolator  
);
 
CMR_EXPORT
CMR_ERROR CMRsepaFindBinaryRepresentativesSubmatrix(
  CMR* cmr,               
  CMR_SEPA* sepa,         
  CMR_CHRMAT* matrix,     
  CMR_CHRMAT* transpose,  
  CMR_SUBMAT* submatrix,  
  bool* pswapped,         
  CMR_SUBMAT** pviolator  
);
 
CMR_EXPORT
CMR_ERROR CMRsepaGetRepresentatives(
  CMR_SEPA* sepa,           
  size_t reprRows[2][3],    
  size_t reprColumns[2][3]  
);
 
CMR_EXPORT
CMR_ERROR CMRsepaGetProjection(
  CMR_SEPA* sepa,         
  size_t part,            
  size_t* rowsToPart,     
  size_t* columnsToPart,  
  size_t* pnumPartRows,   
  size_t* pnumPartColumns 
);
 
CMR_EXPORT
CMR_ERROR CMRsepaCheckTernary(
  CMR* cmr,               
  CMR_SEPA* sepa,         
  CMR_CHRMAT* matrix,     
  bool* pisTernary,       
  CMR_SUBMAT** pviolator  
);
 
CMR_EXPORT
CMR_ERROR CMRsepaCheckTernarySubmatrix(
  CMR* cmr,               
  CMR_SEPA* sepa,         
  CMR_CHRMAT* matrix,     
  CMR_SUBMAT* submatrix,  
  bool* pisTernary,       
  CMR_SUBMAT** pviolator  
);
 
CMR_EXPORT
CMR_ERROR CMRonesumCompose(
  CMR* cmr,               
  size_t numMatrices,     
  CMR_CHRMAT** matrices,  
  CMR_CHRMAT** presult    
);
 
CMR_EXPORT
CMR_ERROR CMRtwosumCompose(
  CMR* cmr,                     
  CMR_CHRMAT* first,            
  CMR_CHRMAT* second,           
  size_t* firstSpecialRows,     
  size_t* firstSpecialColumns,  
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns, 
  int8_t characteristic,        
  CMR_CHRMAT** presult          
);
 
CMR_EXPORT
CMR_ERROR CMRtwosumDecomposeFirst(
  CMR* cmr,                   
  CMR_CHRMAT* matrix,         
  CMR_SEPA* sepa,             
  CMR_CHRMAT** pfirst,        
  size_t* firstRowsOrigin,    
  size_t* firstColumnsOrigin, 
  size_t* rowsToFirst,        
  size_t* columnsToFirst,     
  size_t* firstSpecialRows,   
  size_t* firstSpecialColumns 
);
 
CMR_EXPORT
CMR_ERROR CMRtwosumDecomposeSecond(
  CMR* cmr,                     
  CMR_CHRMAT* matrix,           
  CMR_SEPA* sepa,               
  CMR_CHRMAT** psecond,         
  size_t* secondRowsOrigin,     
  size_t* secondColumnsOrigin,  
  size_t* rowsToSecond,         
  size_t* columnsToSecond,      
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns  
);
 
CMR_EXPORT
CMR_ERROR CMRdeltasumCompose(
  CMR* cmr,                     
  CMR_CHRMAT* first,            
  CMR_CHRMAT* second,           
  size_t* firstSpecialRows,     
  size_t* firstSpecialColumns,  
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns, 
  int8_t characteristic,      
  CMR_CHRMAT** presult        
);
 
CMR_EXPORT
CMR_ERROR CMRdeltasumDecomposeEpsilon(
  CMR* cmr,               
  CMR_CHRMAT* matrix,     
  CMR_CHRMAT* transpose,  
  CMR_SEPA* sepa,         
  char* pepsilon          
);
 
CMR_EXPORT
CMR_ERROR CMRdeltasumDecomposeFirst(
  CMR* cmr,                   
  CMR_CHRMAT* matrix,         
  CMR_SEPA* sepa,             
  char epsilon,               
  CMR_CHRMAT** pfirst,        
  size_t* firstRowsOrigin,    
  size_t* firstColumnsOrigin, 
  size_t* rowsToFirst,        
  size_t* columnsToFirst,     
  size_t* firstSpecialRows,   
  size_t* firstSpecialColumns 
);
 
CMR_EXPORT
CMR_ERROR CMRdeltasumDecomposeSecond(
  CMR* cmr,                     
  CMR_CHRMAT* matrix,           
  CMR_SEPA* sepa,               
  char epsilon,                 
  CMR_CHRMAT** psecond,         
  size_t* secondRowsOrigin,     
  size_t* secondColumnsOrigin,  
  size_t* rowsToSecond,         
  size_t* columnsToSecond,      
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns  
);
 
CMR_EXPORT
CMR_ERROR CMRysumCompose(
  CMR* cmr,                     
  CMR_CHRMAT* first,            
  CMR_CHRMAT* second,           
  size_t* firstSpecialRows,     
  size_t* firstSpecialColumns,  
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns, 
  int8_t characteristic,      
  CMR_CHRMAT** presult        
);
 
CMR_EXPORT
CMR_ERROR CMRysumDecomposeEpsilon(
  CMR* cmr,               
  CMR_CHRMAT* matrix,     
  CMR_CHRMAT* transpose,  
  CMR_SEPA* sepa,         
  char* pepsilon          
);
 
 
CMR_EXPORT
CMR_ERROR CMRysumDecomposeFirst(
  CMR* cmr,                   
  CMR_CHRMAT* matrix,         
  CMR_SEPA* sepa,             
  char epsilon,               
  CMR_CHRMAT** pfirst,        
  size_t* firstRowsOrigin,    
  size_t* firstColumnsOrigin, 
  size_t* rowsToFirst,        
  size_t* columnsToFirst,     
  size_t* firstSpecialRows,   
  size_t* firstSpecialColumns 
);
 
CMR_EXPORT
CMR_ERROR CMRysumDecomposeSecond(
  CMR* cmr,                     
  CMR_CHRMAT* matrix,           
  CMR_SEPA* sepa,               
  char epsilon,                 
  CMR_CHRMAT** psecond,         
  size_t* secondRowsOrigin,     
  size_t* secondColumnsOrigin,  
  size_t* rowsToSecond,         
  size_t* columnsToSecond,      
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns  
);
 
CMR_EXPORT
CMR_ERROR CMRthreesumCompose(
  CMR* cmr,                     
  CMR_CHRMAT* first,            
  CMR_CHRMAT* second,           
  size_t* firstSpecialRows,     
  size_t* firstSpecialColumns,  
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns, 
  int8_t characteristic,        
  CMR_CHRMAT** presult          
);
 
 
CMR_EXPORT
CMR_ERROR CMRthreesumDecomposeSearchConnecting(
  CMR* cmr,               
  CMR_CHRMAT* matrix,     
  CMR_CHRMAT* transpose,  
  CMR_SEPA* sepa,         
  size_t* specialRows,    
  size_t* specialColumns, 
  char* pgamma,           
  char* pbeta             
);
 
CMR_EXPORT
CMR_ERROR CMRthreesumDecomposeSignConnecting(
  CMR* cmr,               
  CMR_CHRMAT* matrix,     
  CMR_CHRMAT* transpose,  
  CMR_SEPA* sepa,         
  size_t* specialRows,    
  size_t* specialColumns, 
  char* pgamma,           
  char* pbeta             
);
 
 
CMR_EXPORT
CMR_ERROR CMRthreesumDecomposeFirst(
  CMR* cmr,                   
  CMR_CHRMAT* matrix,         
  CMR_SEPA* sepa,             
  size_t* specialRows,        
  size_t* specialColumns,     
  char beta,                  
  CMR_CHRMAT** pfirst,        
  size_t* firstRowsOrigin,    
  size_t* firstColumnsOrigin, 
  size_t* rowsToFirst,        
  size_t* columnsToFirst,     
  size_t* firstSpecialRows,   
  size_t* firstSpecialColumns 
);
 
CMR_EXPORT
CMR_ERROR CMRthreesumDecomposeSecond(
  CMR* cmr,                     
  CMR_CHRMAT* matrix,           
  CMR_SEPA* sepa,               
  size_t* specialRows,          
  size_t* specialColumns,       
  char gamma,                   
  CMR_CHRMAT** psecond,         
  size_t* secondRowsOrigin,     
  size_t* secondColumnsOrigin,  
  size_t* rowsToSecond,         
  size_t* columnsToSecond,      
  size_t* secondSpecialRows,    
  size_t* secondSpecialColumns  
);
 
#ifdef __cplusplus
}
#endif
 
#endif /* CMR_SEPARATION_H */