跳过导航链接 | |
退出打印视图 | |
![]() |
Oracle Solaris Studio 12.3:性能分析器 MPI 教程 Oracle Solaris Studio 12.3 Information Library (简体中文) |
本附录提供了本教程中使用的样例代码和 Makefile。
示例 B-1 ring_c.c 的样例代码
/* * Copyright (c) 2004-2006 The Trustees of Indiana University and Indiana * University Research and Technology * Corporation. All rights reserved. * Copyright (c) 2006 Cisco Systems, Inc. All rights reserved. * * Simple ring test program */ #include <stdio.h> #include "mpi.h" int main(int argc, char *argv[]) { int rank, size, next, prev, message, tag = 201; /* Start up MPI */ MPI_Init(&argc, &argv); MPI_Comm_rank(MPI_COMM_WORLD, &rank); MPI_Comm_size(MPI_COMM_WORLD, &size); /* Calculate the rank of the next process in the ring. Use the modulus operator so that the last process "wraps around" to rank zero. */ next = (rank + 1) % size; prev = (rank + size - 1) % size; /* If we are the "master" process (i.e., MPI_COMM_WORLD rank 0), put the number of times to go around the ring in the message. */ if (0 == rank) { message = 10; printf("Process 0 sending %d to %d, tag %d (%d processes in ring)\n", message, next, tag, size); MPI_Send(&message, 1, MPI_INT, next, tag, MPI_COMM_WORLD); printf("Process 0 sent to %d\n", next); } /* Pass the message around the ring. The exit mechanism works as follows: the message (a positive integer) is passed around the ring. Each time it passes rank 0, it is decremented. When each processes receives a message containing a 0 value, it passes the message on to the next process and then quits. By passing the 0 message first, every process gets the 0 message and can quit normally. */ while (1) { MPI_Recv(&message, 1, MPI_INT, prev, tag, MPI_COMM_WORLD, MPI_STATUS_IGNORE); if (0 == rank) { --message; printf("Process 0 decremented value: %d\n", message); } MPI_Send(&message, 1, MPI_INT, next, tag, MPI_COMM_WORLD); if (0 == message) { printf("Process %d exiting\n", rank); break; } } /* The last process does one extra send to process 0, which needs to be received before the program can exit */ if (0 == rank) { MPI_Recv(&message, 1, MPI_INT, prev, tag, MPI_COMM_WORLD, MPI_STATUS_IGNORE); } /* All done */ MPI_Finalize(); return 0; }
以下代码显示了可用于生成 ring_c.c 的 Makefile 的内容,该 Makefile 包含在样例代码中。
示例 B-2 用于生成 ring 程序的 Makefile
# # Copyright (c) 2004-2005 The Trustees of Indiana University and Indiana # University Research and Technology # Corporation. All rights reserved. # Copyright (c) 2004-2005 The University of Tennessee and The University # of Tennessee Research Foundation. All rights # reserved. # Copyright (c) 2004-2005 High Performance Computing Center Stuttgart, # University of Stuttgart. All rights reserved. # Copyright (c) 2004-2005 The Regents of the University of California. # All rights reserved. # Copyright (c) 2006-2007 Sun Microsystems, Inc. All rights reserved. # $COPYRIGHT$ # # Additional copyrights may follow # # $HEADER$ # # Use the Open MPI-provided wrapper compilers. Note that gmake # requires the CXX macro, while other versions of make (such as Sun's # make) require the CCC macro. CC = mpicc CXX = mpic++ CCC = mpic++ F77 = mpif77 FC = mpif90 # Using -g is not necessary, but it is helpful for example programs, # especially if users want to examine them with debuggers. Note that # gmake requires the CXXFLAGS macro, while other versions of make # (such as Sun's make) require the CCFLAGS macro. CFLAGS = -g CXXFLAGS = -g CCFLAGS = -g F77FLAGS = -g FCFLAGS = -g # Example programs to build EXAMPLES = hello_c hello_cxx hello_f77 hello_f90 \ ring_c ring_cxx ring_f77 ring_f90 connectivity_c # Default target. Always build the C example. Only build the others # if Open MPI was build with the relevant language bindings. all: hello_c ring_c connectivity_c @ if test "`ompi_info --parsable | grep bindings:cxx:yes`" != ""; then \ $(MAKE) hello_cxx ring_cxx; \ fi @ if test "`ompi_info --parsable | grep bindings:f77:yes`" != ""; then \ $(MAKE) hello_f77 ring_f77; \ fi @ if test "`ompi_info --parsable | grep bindings:f90:yes`" != ""; then \ $(MAKE) hello_f90 ring_f90; \ fi # The usual "clean" target clean: rm -f $(EXAMPLES) *~ *.o # Don't rely on default rules for the fortran examples hello_f77: hello_f77.f $(F77) $(F77FLAGS) $^ -o $@ ring_f77: ring_f77.f $(F77) $(F77FLAGS) $^ -o $@ hello_f90: hello_f90.f90 $(FC) $(FCFLAGS) $^ -o $@ ring_f90: ring_f90.f90 $(FC) $(FCFLAGS) $^ -o $@