biggstab

#include "bicgstab.h"
#include "ui_bicgstab.h"

BiCGStab::BiCGStab(QWidget *parent) :
    QMainWindow(parent),
    ui(new Ui::BiCGStab)
{
    ui->setupUi(this);

    QPixmap qpixmap("imagens/teste2.png");
    ui->label->setPixmap(qpixmap);
}

BiCGStab::~BiCGStab()
{
    delete ui;
}

void BiCGStab::executar(QMainWindow *pai)
{
    ui->lineEdit->clear();
    ui->lineEditErroRel->setText("0.00001");
    ui->lineEdit_2->clear();
    ui->lineEdit_3->setText(QDir::currentPath());
    ui->lineEditNomeArqSai->setText("RESULTADO_BiCGStab.txt");

    this->retorno = pai;
    this->show();
}

void BiCGStab::on_pushButtonSair_clicked()
{
    QString msg = "Tem certeza que deseja sair do programa?";

    if(!QMessageBox::question(this, "Confirmar Saída", msg, "OK", "Cancelar"))
    {
        this->close();
        this->retorno->close();
    }
}

void BiCGStab::on_pushButtonVoltar_clicked()
{
    retorno->setVisible(true);
    this->close();
}

void BiCGStab::on_pushButtonProcurar_clicked()
{
    QString path;

    path = QDir::currentPath();

    QString file = QFileDialog::getOpenFileName(this, tr("Open File"), path,
                    tr("Arquivos Texto (*.txt)"));

    if(!file.isEmpty())
    {
        ui->lineEdit->setText(file);
    }
}

void BiCGStab::on_pushButtonProcurar_2_clicked()
{
    QString path;

    path = QDir::currentPath();

    QString file = QFileDialog::getOpenFileName(this, tr("Open File"), path,
                    tr("Arquivos Texto (*.txt)"));

    if(!file.isEmpty())
    {
        ui->lineEdit_2->setText(file);
    }
}

void BiCGStab::on_pushButtonProcurar_3_clicked()
{
    QString path = QDir::currentPath();

    QString dir = QFileDialog::getExistingDirectory(this, tr("Diretório"), path);

    if (!dir.isEmpty())
    {
        ui->lineEdit_3->setText(dir);
    }
}

void BiCGStab::on_commandLinkButton_clicked()
{
    if(!ui->lineEdit->text().isEmpty() && !ui->lineEdit_2->text().isEmpty() && !ui->lineEditErroRel->text().isEmpty() && !ui->lineEdit_3->text().isEmpty())
    {
        try
        {
            this->realizaOperacoes();

            QMessageBox::warning(this, tr("Sucesso!"), "A operação foi executada com sucesso! Confira os dados no arquivo " + ui->lineEditNomeArqSai->text() + " gerado no diretório escolhido.");

            retorno->setVisible(true);
            this->close();
        }

        catch(string erro)
        {
            QString descricao = QString::fromStdString(erro);

            QMessageBox::warning(this, tr("Erro"), descricao);
        }
    }

    else
    {
        QMessageBox::warning(this, tr("Atenção!"), "Favor preencher todos os campos antes de executar.");
    }
}

void BiCGStab::realizaOperacoes()
{
    try
    {
        string nomeArquivo, linha, dados;
        int cont = 0, qtdlinhas, qtdcolunas;
        float erro = ui->lineEditErroRel->text().toFloat();
        nIteracoes = 4000;
        //float aux;

        nomeArquivo = ui->lineEdit->text().toStdString();

        ifstream arquivo;

        arquivo.open(nomeArquivo.c_str());

        // verificando erro de abertura do arquivo de entrada
        if(!arquivo.is_open())
        {
            string a = "ERRO - Arquivo nao existe ou nao se encontra no mesmo diretorio do executavel!";
            throw a;
        }

        //qtdlinhas = 0;
        //qtdcolunas = 0;

        arquivo>> qtdlinhas;
        arquivo>> qtdcolunas;                      
        //getline(arquivo, linha);

        dados = "";

        /*
        while(!arquivo.eof()) //repeticao para destrinchar os dados da matriz
        {
            //qtdlinhas++;

            dados += linha + " ";

            getline(arquivo, linha);

            while(linha == "" || linha[0] == ' ') //le do arquivo ate encontrar uma linha que nao seja vazia
            {
                getline(arquivo, linha);

                if(arquivo.eof())
                    break;
            }
        }

        if(linha != "") //se a ultima linha lida do arquivo nao for vazia, entao 'dados' recebe linha para ficar completo com a ultima linha da matriz
        {
            //qtdlinhas++;
            dados += linha + " ";
        }
        */

        for(cont = 0; cont < (qtdlinhas * qtdcolunas); cont++)
        {
            arquivo>> linha;
            dados += linha + " ";
        }

        arquivo.close();


        matrizA = new Matriz(qtdlinhas, qtdcolunas, dados);

        if(matrizA->getQuantidadeColunas() == 0)
        {
            string erro = "ESTE ARQUIVO NAO CONTEM DADOS!\n\n\n";
            throw erro;
        }
        /////////////////////////////////////////////////////////////////



        // abrindo um arquivo de entrada - leitura do vetor b
        nomeArquivo = ui->lineEdit_2->text().toStdString();

        ifstream arquivo2;

        arquivo2.open(nomeArquivo.c_str());

        // verificando erro de abertura do arquivo de entrada
        if(!arquivo2.is_open())
        {
            string a = "ERRO - Arquivo nao existe ou nao se encontra no mesmo diretorio do executavel!";
            throw a;
        }

        qtdlinhas = 0;
        qtdcolunas = 1; //o vetor b sempre tera apenas uma coluna

        /*
        getline(arquivo2, linha); //faz a leitura da primeira linha na variavel 'linha'

        while(linha == "" || linha[0] == ' ') //le do arquivo ate encontrar uma linha que nao seja vazia
        {
            getline(arquivo2, linha);

            if(arquivo2.eof())
                break;
        }

        dados = "";

        while(!arquivo2.eof()) //repeticao para destrinchar os dados da matriz
        {
            qtdlinhas++;

            dados += linha + " ";

            getline(arquivo2, linha);

            while(linha == "" || linha[0] == ' ') //le do arquivo ate encontrar uma linha que nao seja vazia
            {
                getline(arquivo2, linha);

                if(arquivo2.eof())
                    break;
            }
        }

        if(linha != "") //se a ultima linha lida do arquivo nao for vazia, entao 'dados' recebe linha ficar completo com a ultima linha da matriz
        {
            qtdlinhas++;
            dados += linha + " ";
        }
        */

        dados = "";
        arquivo2>> linha;

        while(!arquivo2.eof())
        {
            dados += linha + " ";
            qtdlinhas++;
            arquivo2>> linha;
        }

        dados += linha + " ";
        qtdlinhas++;

        arquivo2.close();


        if(matrizA->getQuantidadeLinhas() == qtdlinhas) //a quantidade de linhas da matriz A sempre deve ser igual a quantidade de linhas do vetor b
        {
            matrizB = new Matriz(qtdlinhas, qtdcolunas, dados);
            matrizX = new Matriz(matrizA->getQuantidadeColunas(), 1); 

            //clock_t start, end;

			QString msg = "Se esta máquina possui uma placa gráfica da NVidia, ao clicar em 'Sim', "
				"o método será executado em paralelo pela GPU existente na placa. Dessa forma, poderá haver um grande aumento de desempenho no tempo de execução.\n\n"
				"Deseja ativar o processamento paralelo?";

			if(!QMessageBox::question(this, "Confirmar Saída", msg, "Sim", "Não"))
			{				
				matrizX = bicgstab_kernel(matrizA, matrizB, nIteracoes, erro, tempoExec); //executa no device							
			}

			else
			{
				//start = clock();
				matrizX = matrizA->BiCGStab(matrizB, nIteracoes, erro, tempoExec);
				//matrizX = matrizA->BiCGStab((*matrizB), nIteracoes, erro);
				//end = clock();
				//tempoExec = ((double)(end - start) / ((double)CLOCKS_PER_SEC));
			}



            ofstream arquivoSaida;
            dados = ui->lineEdit_3->text().toStdString();
            dados += "/" + ui->lineEditNomeArqSai->text().toStdString();

            arquivoSaida.open(dados.c_str());

            if(!arquivoSaida.is_open())
            {
                string a = "ERRO NA ABERTURA DO ARQUIVO dados.txt";
                throw a;
            }


            arquivoSaida<< "Quantidade Iterações: " << nIteracoes;
			arquivoSaida<< "\nResíduo: " << erro;
            arquivoSaida<< "\nTempo de execução do método: " << setw(10) << setfill('0') << setprecision(10) << tempoExec << " segundos";

            arquivoSaida<< "\n\nvetor x:\n\n";
            matrizX->imprimeVetor(arquivoSaida); //esse metodo imprime os dados no arquivo passado no parametro

            arquivoSaida<< "\n\nMatriz A:\n\n";
            matrizA->imprimeMatriz(arquivoSaida);
            arquivoSaida<< "\n\nvetor b:\n\n";
            matrizB->imprimeMatriz(arquivoSaida);

            arquivoSaida.close();
        }

        else
        {
            string erro = "IMPOSSIVEL RESOLVER ESTE SISTEMA! A QUANTIDADE DE LINHAS DA MATRIZ 'A' DEVE SER IGUAL A "
                          "QUANTIDADE DE LINHAS DO VETOR 'b'. POR FAVOR, CORRIJA OS DADOS NOS ARQUIVOS!";

            throw erro;
        }

        delete matrizA;
        delete matrizB;
        delete matrizX;
    }

    catch(string erro)
    {
        throw erro;
    }
}