Epetra Package Browser (Single Doxygen Collection) Development
Loading...
Searching...
No Matches
threaded_main.cpp
Go to the documentation of this file.
1//@HEADER
2// ************************************************************************
3//
4// Epetra: Linear Algebra Services Package
5// Copyright 2011 Sandia Corporation
6//
7// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
8// the U.S. Government retains certain rights in this software.
9//
10// Redistribution and use in source and binary forms, with or without
11// modification, are permitted provided that the following conditions are
12// met:
13//
14// 1. Redistributions of source code must retain the above copyright
15// notice, this list of conditions and the following disclaimer.
16//
17// 2. Redistributions in binary form must reproduce the above copyright
18// notice, this list of conditions and the following disclaimer in the
19// documentation and/or other materials provided with the distribution.
20//
21// 3. Neither the name of the Corporation nor the names of the
22// contributors may be used to endorse or promote products derived from
23// this software without specific prior written permission.
24//
25// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36//
37// Questions? Contact Michael A. Heroux (maherou@sandia.gov)
38//
39// ************************************************************************
40//@HEADER
41
42
43// Petra_Comm Test routine
44
45#include <math.h>
46#include <unistd.h>
47#include <string.h>
48#include <ctype.h>
49#include <stdlib.h>
50#ifdef PETRA_MPI
51#include <mpi.h>
52#endif
53#include "Petra_Comm.h"
54#include "Petra_Time.h"
55
56int main(int argc, char *argv[]) {
57
58 // Initialize MPI
59
60 MPI_Init(&argc,&argv);
61 int size, rank; // Number of MPI processes, My process ID
62
63 MPI_Comm_size(MPI_COMM_WORLD, &size);
64 MPI_Comm_rank(MPI_COMM_WORLD, &rank);
65 // I'm alive !!!
66
67
68
69 Petra_Comm & petracomm = *new Petra_Comm( MPI_COMM_WORLD );
70 int MyPID = petracomm.MyPID();
71 int NumProc = petracomm.NumProc();
72 cout << "Processor "<<MyPID<<" of " << NumProc << " is alive."<<endl;
73
74 if (NumProc!=2) {
75 cout << " This special routine only works for 2 processors " << endl;
76 abort();
77 }
78
79
80 Petra_Comm * other_comm;
81 MPI_Status status;
82
83 unsigned int icomm = &petracomm;
84
85 if (MyPID==1) cout << "Address of Petra_Comm object on PE 1 = " << &petracomm << endl;
86
87 if (MyPID==0) MPI_Recv((void *) &other_comm, 1, MPI_UNSIGNED, 1, 99, MPI_COMM_WORLD, &status);
88 else MPI_Send ( (void *) &icomm, 1, MPI_UNSIGNED, 0, 99, MPI_COMM_WORLD);
89
90 if (MyPID==0) cout << "Address of other Petra_Comm object on PE 0 = " << other_comm << endl;
91
92 int otherPID = other_comm->MyPID();
93
94 if (MyPID==0) cout << "Processor "<<MyPID<<" of " << NumProc
95 << " has a neighbor processor with ID "
96 << otherPID << " of " << other_comm->NumProc() <<endl;
97
98 delete &petracomm;
99 MPI_Finalize();
100
101 return 0;
102}
103
int main(int argc, char *argv[])