int main()
{
int matrizA[M][M];
lerMatriz(matrizA);
int resultado = processarMatriz(matrizA);
imprimirResultado(resultado);
getchar();
return 0;
}
int main (int argc, char* argv[]){
// Variáveis auxiliares para leitura do arquivo
int i, j, k; // Vars. controladoras de iterações
int auxInt; // Armazena temporariamente a Ordem da Matriz
short auxShort; // Armazena temporariamente a fila da matriz a ser avaliada
double auxDouble; // Armazena temporariamente o erro máx. permitido e vals. das matrizes
long auxLong; // Armazena temporariamente o núm. max e num real de iterações
short opEscolhida; // Opção escolhida na tela inicial
char* caminhoMat = telaInicial(-1,&opEscolhida); // Armazena o caminho para o arquivo da matriz
double* resultado = NULL; // Armazena o resultado do método de jacobi richardson
double* resIt = NULL; // Armazena o resultado das NUM_EXEC iterações do método jacobi richardson
double* x0 = NULL; // Vetor x0 de chute inicial
FILE* arqMat = NULL; // Ponteiro para o arquivo de matrizes que será aberto
clock_t begin, end; // Clock de inicio e fim da execucao das iteracoes do jacobi richardson
double* time_spent = (double*)calloc(NUM_EXEC,sizeof(double)); // Tempo gasto para rodar as iteracoes (em segundos)
double mediaTimeSpent = 0; // Media dos tempos gastos
double desvioPadrao;
// Cria a estrutura que armazenará as matrizes
MAT_ENTRADA* m = criarMatEntrada();
arqMat = fopen(caminhoMat,"r");
if(arqMat == NULL){
printf("Erro ao abrir arquivo.");
}
else{
// Le as primeiras informações do arquivo
fscanf(arqMat, "%d", &auxInt); // Lê a ordem da matriz
fscanf(arqMat, "%hd", &auxShort); // Lê a fila que será avaliada
fscanf(arqMat, "%lf", &auxDouble); // Lê o erro permitido
fscanf(arqMat, "%ld", &auxLong); // Lê o núm. max de iterações
// Salva os valores na estrutura
inicializarValsMatEntrada(m,auxInt,auxShort,auxDouble,auxLong,0);
// Aloca espaço para as matrizes A e B e diagonalAux, agora que sabemos a dimensão
alocarMatA(m);
alocarMatB(m);
alocarDiagonalAux(m);
// Insere os elementos na Matriz A
for(i=0;i<getOrdem(m);i++){
for(j=0;j<getOrdem(m);j++){
fscanf(arqMat,"%lf",&auxDouble);
inserirElemMatA(m,auxDouble,i,j);
}
}
// Insere os elementos na Matriz B
for(i=0;i<getOrdem(m);i++){
fscanf(arqMat,"%lf",&auxDouble);
inserirElemMatB(m,auxDouble,i);
}
// Fecha o arquivo agora que não precisamos mais dele
fclose(arqMat);
}
resIt = (double*)calloc(getOrdem(m),sizeof(double));
auxLong = 0;
// Prepara as matrizes A e B dividindo os valores pela diagonal da A
prepararMatrizes(m);
// 10 iterações do método de jacobiRichardson
for(i=0;i<NUM_EXEC;i++){
//printf("Iniciando execucao %d do Jacobi Richardson.\n", i+1);
// Aqui geramos o vetor de chute inicial x0
// Cada pos x0[i] recebe um valores double entre 0 e MAX_DOUBLE
x0 = (double*)malloc(getOrdem(m)*sizeof(double));
srand((unsigned int)time(NULL));
for(k=0;k<getOrdem(m);k++){
x0[k] = ((float)rand()/(float)(RAND_MAX))*MAX_DOUBLE;
}
begin = clock(); // Começa a contar o tempo de execucao
resultado = jacobiRichardson(m,x0); // Execucao do metodo
end = clock(); // Termina de contar contar o tempo de execucao
time_spent[i] = (double)(end - begin) / CLOCKS_PER_SEC; // Calcula o valor em segundos
mediaTimeSpent += time_spent[i];
salvarSaidasItermediarias(m,resultado,opEscolhida);
auxLong += getRealIt(m); // Salva o nº de iteracoes da execução "i"
for(j=0;j<getOrdem(m);j++){
resIt[j] += resultado[j]; // Salva o vetor resultado da execução "i"
}
free(x0);
x0 = NULL;
free(resultado);
resultado = NULL;
//printf("Fim da execucao %d.\n\n", i+1);
}
// Pega a media das 10 iterações
mediaTimeSpent = mediaTimeSpent/NUM_EXEC;
desvioPadrao = calcDesvioPadrao(time_spent,mediaTimeSpent,NUM_EXEC);
auxLong = auxLong/NUM_EXEC;
for(i=0;i<getOrdem(m);i++){
resIt[i] = resIt[i]/NUM_EXEC;
}
// Salva e imprime na tela
salvarSaidaFinal(m,resIt,auxLong,opEscolhida);
imprimirResultado(m,resIt,auxLong);
for(i=0;i<NUM_EXEC;i++){
printf("Tempo de execucao[%d]: %lfs.\n", i, time_spent[i]);
}
printf("Media do tempo de execucao: %lfs.\n", mediaTimeSpent);
printf("Desvio padrao: %lfs.\n",desvioPadrao);
// Libera a memória alocada pela estrutura
destruirMatEntrada(m);
free(time_spent);
return 0;
}
