routines/jacobi/libjacobi.h

#pragma once
#ifndef LIBJACOBI_H_
#define LIBJACOBI_H_

#include <mpi.h>

/*
 * External variables
 * \param verbose Print out extra debugging info
 * \param n Size of matrix
 * \param max_itt Maximum number of iterations
 * \param rank Rank of current processor
 * \param nprocs Number of nodes in computation
 */
//extern int verbose, max_itt, rank, nprocs;

double* jacobi_compute(double* a, double* b, double* x, int D, int d, int offset, double* diff);

double* jacobi_redistribute(MPI_Comm comm, double* a, int D, int* old_d, int* new_d); 

int jacobi_fill_matrix(double** a, double** b, double** x, int D, int d, int offset, int seed, int rank);

void jacobi_final_test(MPI_Comm comm, double* a, double* b, double* x, int D, int* d, int* offset);

void jacobi_set_D(long long int D);

double jacobi_complexity(long long int d);

#endif /* LIBJACOBI_H_ */