Revision 02fbe2fa4695e141a1d011a515b7566b9f107d03 authored by Leif Strand on 07 November 2006, 01:33:06 UTC, committed by Leif Strand on 07 November 2006, 01:33:06 UTC
approach to build/install. Also, temporarily deleted 'Exchanger' and all references to it, to avoid confusion. We can always go back in time and retrieve it later.
1 parent c20d5e5
array2d.cc
// -*- C++ -*-
//
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
//<LicenseText>
//
// CitcomS.py by Eh Tan, Eun-seo Choi, and Pururav Thoutireddy.
// Copyright (C) 2002-2005, California Institute of Technology.
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
//
//</LicenseText>
//
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
#include <portinfo>
#include <string>
#include <mpi.h>
#include "journal/diagnostics.h"
#include "Array2D.h"
using namespace std;
using namespace Exchanger;
int main(int argc, char** argv)
{
const string name = "Array2DTest";
journal::info_t info(name);
info.activate();
const int dim = 3;
MPI_Init(&argc, &argv);
MPI_Comm comm = MPI_COMM_WORLD;
int rank;
MPI_Comm_rank(comm, &rank);
Array2D<double,dim> d1(2);
Array2D<int,1> i1(5);
if (rank == 0) {
// testing vector-like behaviors
info << "c'tor" << journal::endl;
d1.print(name);
i1.print(name);
info << "copy c'tor" << journal::endl;
d1[0][0] = 1;
Array2D<double,3> d2(d1);
d2.print(name);
info << "op'tor [][]" << journal::endl;
for (int i=0; i<dim; i++)
d2[i][1] = 10*i;
info << d2[1][1] << journal::endl;
d2.print(name);
info << "swap" << journal::endl;
d2.swap(d1);
d1.print(name);
d2.print(name);
info << "resize" << journal::endl;
i1.resize(3);
i1.print(name);
info << "push_back" << journal::endl;
i1.push_back(3);
i1.push_back(4);
i1.print(name);
d1.push_back(vector<double>(dim, 5));
d1.print(name);
info << "reserve shrink capacity" << journal::endl;
i1.reserve(1000);
info << "capacity = " << i1.capacity() << journal::endl;
i1.shrink();
info << "capacity = " << i1.capacity() << journal::endl;
}
// testing send/receive ...
if (rank == 0) {
d1.sendSize(comm, 1);
} else if(rank == 1) {
info << "sendSize/receiveSize" << journal::endl;
Array2D<double,dim> d3;
info << "received size = " << d3.recvSize(comm, 0) << journal::endl;
}
if (rank == 0) {
d1.send(comm, 1);
} else if(rank == 1) {
info << "blocking send -> blocking receive" << journal::endl;
Array2D<double,dim> d3;
d3.recv(comm, 0);
d3.print(name);
}
if (rank == 0) {
MPI_Request request;
MPI_Status status;
d1.send(comm, 1, request);
MPI_Wait(&request, &status);
} else if(rank == 1) {
info << "non-blocking send -> non-blocking receive" << journal::endl;
Array2D<double,dim> d3;
MPI_Request request;
MPI_Status status;
d3.recv(comm, 0, request);
MPI_Wait(&request, &status);
d3.print(name);
}
if (rank == 0) {
Array2D<double,dim> d3(dim*dim);
for(int i=0; i<dim; i++)
for(int j=0; j<dim*dim; j++)
d3[i][j] = (i+1)*j;
//d3.print(name);
vector<MPI_Request> request(dim);
vector<MPI_Status> status(dim);
for(int i=0; i<dim; i++)
d3.send(comm, 1, i*dim, dim, request[i]);
MPI_Waitall(dim, &request[0], &status[0]);
} else if(rank == 1) {
info << "non-blocking partial send -> non-blocking partial receive" << journal::endl;
Array2D<double,dim> d3(dim*dim);
vector<MPI_Request> request(dim);
vector<MPI_Status> status(dim);
for(int i=0; i<dim; i++)
d3.recv(comm, 0, i*dim, dim, request[i]);
MPI_Waitall(dim, &request[0], &status[0]);
d3.print(name);
}
MPI_Finalize();
}
// version
// $Id$
// End of file
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...