Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021 #include "ColPackHeaders.h"
00022
00023 using namespace std;
00024
00025 namespace ColPack
00026 {
00027
00028
00029 BipartiteGraphBicoloringInterface::~BipartiteGraphBicoloringInterface()
00030 {
00031 BipartiteGraphBicoloring::Clear();
00032
00033 Seed_reset();
00034 }
00035
00036
00037
00038 void BipartiteGraphBicoloringInterface::Clear()
00039 {
00040 BipartiteGraphBicoloring::Clear();
00041
00042 return;
00043 }
00044
00045
00046
00047 void BipartiteGraphBicoloringInterface::Reset()
00048 {
00049 BipartiteGraphBicoloring::Reset();
00050
00051 return;
00052 }
00053
00054
00055
00056 void BipartiteGraphBicoloringInterface::GenerateSeedJacobian(double*** dp3_LeftSeed, int *ip1_LeftSeedRowCount, int *ip1_LeftSeedColumnCount, double*** dp3_RightSeed, int *ip1_RightSeedRowCount, int *ip1_RightSeedColumnCount, string s_OrderingVariant, string s_BicoloringVariant) {
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066 Bicoloring(s_OrderingVariant, s_BicoloringVariant);
00067
00068
00069 *dp3_LeftSeed = GetLeftSeedMatrix(ip1_LeftSeedRowCount, ip1_LeftSeedColumnCount);
00070 *dp3_RightSeed = GetRightSeedMatrix(ip1_RightSeedRowCount, ip1_RightSeedColumnCount);
00071
00072 }
00073
00074 void BipartiteGraphBicoloringInterface::GenerateSeedJacobian_unmanaged(double*** dp3_LeftSeed, int *ip1_LeftSeedRowCount, int *ip1_LeftSeedColumnCount, double*** dp3_RightSeed, int *ip1_RightSeedRowCount, int *ip1_RightSeedColumnCount, string s_OrderingVariant, string s_BicoloringVariant) {
00075
00076
00077 Bicoloring(s_OrderingVariant, s_BicoloringVariant);
00078
00079
00080 *dp3_LeftSeed = GetLeftSeedMatrix_unmanaged(ip1_LeftSeedRowCount, ip1_LeftSeedColumnCount);
00081 *dp3_RightSeed = GetRightSeedMatrix_unmanaged(ip1_RightSeedRowCount, ip1_RightSeedColumnCount);
00082 }
00083
00084 int BipartiteGraphBicoloringInterface::Bicoloring(string s_OrderingVariant, string s_BicoloringVariant) {
00085 m_T_Timer.Start();
00086 int i_OrderingStatus = OrderVertices(s_OrderingVariant);
00087 m_T_Timer.Stop();
00088 m_d_OrderingTime = m_T_Timer.GetWallTime();
00089
00090 if(i_OrderingStatus != _TRUE)
00091 {
00092 cerr<<endl;
00093 cerr<<s_OrderingVariant<<" Ordering Failed";
00094 cerr<<endl;
00095
00096 return(1);
00097 }
00098
00099 s_BicoloringVariant = toUpper(s_BicoloringVariant);
00100 m_T_Timer.Start();
00101
00102 int i_ColoringStatus;
00103 if(s_BicoloringVariant == "IMPLICIT_COVERING__STAR_BICOLORING") {
00104 i_ColoringStatus = ImplicitCoveringStarBicoloring();
00105 } else if (s_BicoloringVariant == "EXPLICIT_COVERING__STAR_BICOLORING") {
00106 i_ColoringStatus = ExplicitCoveringStarBicoloring();
00107 } else if (s_BicoloringVariant == "EXPLICIT_COVERING__MODIFIED_STAR_BICOLORING") {
00108 i_ColoringStatus = ExplicitCoveringModifiedStarBicoloring();
00109 } else if (s_BicoloringVariant == "IMPLICIT_COVERING__GREEDY_STAR_BICOLORING") {
00110 i_ColoringStatus = ImplicitCoveringGreedyStarBicoloring();
00111 } else {
00112 cout<<" Unknown Bicoloring Method "<<s_BicoloringVariant<<". Please use a legal Method."<<endl;
00113 m_T_Timer.Stop();
00114 m_d_ColoringTime = m_T_Timer.GetWallTime();
00115 return (_FALSE);
00116 }
00117
00118 m_T_Timer.Stop();
00119 m_d_ColoringTime = m_T_Timer.GetWallTime();
00120 return(i_ColoringStatus);
00121 }
00122
00123 BipartiteGraphBicoloringInterface::BipartiteGraphBicoloringInterface(int i_type, ...) {
00124
00125 Clear();
00126
00127 if (i_type == SRC_WAIT) return;
00128
00129
00130 va_list ap;
00131 va_start(ap,i_type);
00132
00133 if (i_type == SRC_MEM_ADOLC) {
00134
00135 unsigned int ** uip2_JacobianSparsityPattern = va_arg(ap,unsigned int **);
00136 int i_RowCount = va_arg(ap,int);
00137 int i_ColumnCount = va_arg(ap,int);
00138
00139 BuildBPGraphFromRowCompressedFormat(uip2_JacobianSparsityPattern, i_RowCount, i_ColumnCount);
00140 }
00141 else if (i_type == SRC_MEM_ADIC) {
00143 cerr<<"ERR: GraphColoringInterface(): s_inputSource \"ADIC\" is not supported yet"<<endl;
00144
00145 va_end(ap);
00146 return;
00147 }
00148 else if (i_type == SRC_FILE) {
00149
00150 string s_InputFile ( va_arg(ap,char *) );
00151 string s_fileFormat ( va_arg(ap,char *) );
00152
00153 ReadBipartiteGraph(s_InputFile, s_fileFormat);
00154 }
00155 else {
00156 cerr<<"ERR: BipartiteGraphBicoloringInterface(): i_type =\""<< i_type <<"\" unknown or unspecified"<<endl;
00157
00158 va_end(ap);
00159 return;
00160 }
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220
00221
00222 va_end(ap);
00223 return;
00224 }
00225
00226 double** BipartiteGraphBicoloringInterface::GetLeftSeedMatrix(int* ip1_LeftSeedRowCount, int* ip1_LeftSeedColumnCount) {
00227 return BipartiteGraphBicoloring::GetLeftSeedMatrix(ip1_LeftSeedRowCount, ip1_LeftSeedColumnCount);
00228 }
00229
00230 double** BipartiteGraphBicoloringInterface::GetRightSeedMatrix(int* ip1_RightSeedRowCount, int* ip1_RightSeedColumnCount) {
00231 return BipartiteGraphBicoloring::GetRightSeedMatrix(ip1_RightSeedRowCount, ip1_RightSeedColumnCount);
00232 }
00233
00234 void BipartiteGraphBicoloringInterface::GetOrderedVertices(vector<int> &output) {
00235 BipartiteGraphOrdering::GetOrderedVertices(output);
00236 }
00237 }