Привет,
моя программа завершает работу с кодом 3. Нет сообщений об ошибках и исключений, и выход не инициируется моим кодом.
Проблема возникает, когда я пытаюсь прочитать очень длинные целочисленные значения из текстового файла (текстовый файл присутствует и правильно открыт, с успешным предварительным чтением).
Я использую очень большие объемы памяти (на самом деле, я думаю, что это может быть причиной, так как я почти уверен, что превышаю лимит памяти в 2 ГБ на процесс). Я также использую библиотеку GMP (или, скорее, MPIR) для умножения бигнумов. Я вполне уверен, что это не проблема с файловым вводом / выводом, поскольку я получил тот же код ошибки в предыдущей версии программы, которая была полностью в памяти.
Система:
MS Visual Studio 2008
MS Windows Vista Home Premium x86
MPIR 2.1.0 rc2
4 ГБ ОЗУ
Откуда может исходить этот код ошибки?
РЕДАКТИРОВАТЬ: это процедура, которая завершается с кодом
void condenseBinSplitFile(const char *sourceFilename, int partCount){
//condense results file into final P and Q
std::string tempFilename;
std::string inputFilename(sourceFilename);
std::string outputFilename(BIN_SPLIT_FILENAME_DATA2);
mpz_class *P = new mpz_class(0);
mpz_class *Q = new mpz_class(0);
mpz_class *PP = new mpz_class(0);
mpz_class *QQ = new mpz_class(0);
FILE *sourceFile;
FILE *resultFile;
fpos_t oldPos;
int swapCount = 0;
while (partCount > 1){
std::cout << partCount << std::endl;
sourceFile = fopen(inputFilename.c_str(), "r");
resultFile = fopen(outputFilename.c_str(), "w");
for (int i=0; i<partCount/2; i++){
//Multiplication order:
//Get Q, skip P
//Get QQ, mul Q and QQ, print Q, delete Q
//Jump back to P, get P
//Mul P and QQ, delete QQ
//Skip QQ, get PP
//Mul P and PP, delete P and PP
//Get Q, skip P
mpz_inp_str(Q->get_mpz_t(), sourceFile, CALC_BASE);
fgetpos(sourceFile, &oldPos);
skipLine(sourceFile);
skipLine(sourceFile);
//Get QQ, mul Q and QQ, print Q, delete Q
mpz_inp_str(QQ->get_mpz_t(), sourceFile, CALC_BASE);
(*Q) *= (*QQ);
mpz_out_str(resultFile, CALC_BASE, Q->get_mpz_t());
fputc('\n', resultFile);
(*Q) = 0;
//Jump back to P, get P
fsetpos(sourceFile, &oldPos);
mpz_inp_str(P->get_mpz_t(), sourceFile, CALC_BASE);
//Mul P and QQ, delete QQ
(*P) *= (*QQ);
(*QQ) = 0;
//Skip QQ, get PP
skipLine(sourceFile);
skipLine(sourceFile);
mpz_inp_str(PP->get_mpz_t(), sourceFile, CALC_BASE);
//Mul P and PP, delete PP, print P, delete P
(*P) += (*PP);
(*PP) = 0;
mpz_out_str(resultFile, CALC_BASE, P->get_mpz_t());
fputc('\n', resultFile);
(*P) = 0;
}
partCount /= 2;
fclose(sourceFile);
fclose(resultFile);
//swap filenames
tempFilename = inputFilename;
inputFilename = outputFilename;
outputFilename = tempFilename;
swapCount++;
}
delete P;
delete Q;
delete PP;
delete QQ;
remove(BIN_SPLIT_FILENAME_RESULTS);
if (swapCount%2 == 0)
rename(sourceFilename, BIN_SPLIT_FILENAME_RESULTS);
else
rename(BIN_SPLIT_FILENAME_DATA2, BIN_SPLIT_FILENAME_RESULTS);
}
EDIT2: версия полностью в памяти, которая также выходит с 3
void binarySplitE(const ULONG first, const ULONG last, mpz_class *P, mpz_class *Q){
//P(first, last) = P(first, mid)*Q(mid, last) + P(mid, last)
//Q(first, last) = Q(first, mid)*Q(mid, last)
if (last - first == 1){
calcP(P, first, last);
calcQ(Q, first, last);
return;
}
ULONG mid = (first+last)/2;
mpz_class *PP = new mpz_class(*P);
mpz_class *QQ = new mpz_class(*Q);
//Calculate P(first, mid) and Q(first, mid)
binarySplitE(first, mid, P, Q);
//Calculate P(mid, last) and Q(mid, last)
binarySplitE(mid, last, PP, QQ);
//P(first, last) == P(first, mid)
*P *= (*QQ);
//P(first, last) == P(first, mid)*Q(mid, last)
*P += (*PP);
//P(first, last) == P(first, mid)*Q(mid, last) + P(mid, last);
//Q(first, last) == Q(first, mid)
*Q *= (*QQ);
//Q(first, last) == Q(first, mid)*Q(mid, last)
delete PP;
delete QQ;
}