vdop.h

Go to the documentation of this file.

#ifndef VDOP_FILE
#define VDOP_FILE
#include <iostream>
#include <mpi.h>                  // MPI
#include <string>
#include <vector>
 
void vddiv(std::vector<double> &x, std::vector<double> const &y,
           std::vector<double> const &z);
 
void vdaxpy(std::vector<double> &x, std::vector<double> const &y,
            double alpha, std::vector<double> const &z );
 
 
double dscapr(std::vector<double> const &x, std::vector<double> const &y);
 
 
inline
double L2_scapr(std::vector<double> const &x, std::vector<double> const &y)
{
    return dscapr(x, y) / x.size();
}
 
 
double par_scalar(std::vector<double> const &x, std::vector<double> const &y,
                   MPI_Comm const& icomm=MPI_COMM_WORLD);
 
 
 
/*  ReadId : Input and broadcast of an integer */
inline
int ReadIn(std::string const &ss = std::string(), MPI_Comm const &icomm = MPI_COMM_WORLD)
{
    MPI_Barrier(icomm);
    int myrank;                             /*  my rank number */
    MPI_Comm_rank(icomm, &myrank);
    int id;
 
    if (myrank == 0) {
        std::cout << "\n\n  " << ss << " :  Which process do you want to debug ? \n";
        std::cin >> id;
    }
    MPI_Bcast(&id, 1, MPI_INT, 0, icomm);
 
    return id;
}
 
//void DebugVector(std::vector<double> const &v);
template <class T>
void DebugVector(std::vector<T> const &v, std::string const &ss = std::string(), MPI_Comm const &icomm = MPI_COMM_WORLD)
{
    MPI_Barrier(icomm);
    int numprocs;                               /*  # processes    */
    MPI_Comm_size(icomm, &numprocs);
    int myrank;                             /*  my rank number */
    MPI_Comm_rank(icomm, &myrank);
 
    int readid = ReadIn(ss);                            /* Read  readid    */
 
    while ( (0 <= readid) && (readid < numprocs) ) {
        if (myrank == readid) {
            std::cout << "\n\n process " << readid;
            std::cout << "\n ....  " << ss << " (nnode = " << v.size() << ")\n";
            for (size_t j = 0; j < v.size(); ++j) {
                std::cout.setf(std::ios::right, std::ios::adjustfield);
                std::cout << v[j] << "   ";
            }
            std::cout << std::endl;
            fflush(stdout);
        }
 
        readid = ReadIn(ss, icomm);                     /* Read  readid    */
    }
    MPI_Barrier(icomm);
    return;
}
 
bool CompareVectors(std::vector<double> const &x, int n, double const y[], double const eps = 0.0);
 
 
template <class T>
std::ostream &operator<<(std::ostream &s, std::vector<T> const &v)
{
    for (auto vp : v) {
        s << vp << "  ";
    }
    return s;
}
 
void ExchangeAll(std::vector<double> const &xin, std::vector<double> &yout, MPI_Comm const &icomm = MPI_COMM_WORLD);
 
void ExchangeAllInPlace(std::vector<double> &xin, MPI_Comm const &icomm = MPI_COMM_WORLD);
 
 
 
#endif