/[MITgcm]/MITgcm/pkg/autodiff/global_max_ad.F
ViewVC logotype

Annotation of /MITgcm/pkg/autodiff/global_max_ad.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.2 - (hide annotations) (download)
Sat Jul 13 03:30:30 2002 UTC (21 years, 10 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint51k_post, checkpoint47e_post, checkpoint57m_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint46l_post, checkpoint57g_pre, checkpoint46g_pre, checkpoint47c_post, checkpoint50c_post, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint46f_post, checkpoint52d_pre, checkpoint57g_post, checkpoint48e_post, checkpoint56b_post, checkpoint50c_pre, checkpoint57y_post, checkpoint46b_post, checkpoint52j_pre, checkpoint51o_pre, checkpoint54d_post, checkpoint54e_post, checkpoint51l_post, checkpoint48i_post, checkpoint57r_post, checkpoint46l_pre, checkpoint57d_post, checkpoint57i_post, checkpoint52l_post, checkpoint52k_post, checkpoint59, checkpoint58, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint58f_post, checkpoint52f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint50b_pre, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint58y_post, checkpoint48b_post, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint51t_post, checkpoint58t_post, c49_autodiff, checkpoint51n_post, checkpoint55i_post, checkpoint58m_post, checkpoint57l_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint47a_post, checkpoint57t_post, checkpoint55c_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint52e_pre, checkpoint57v_post, checkpoint57f_post, checkpoint52e_post, checkpoint51n_pre, checkpoint47d_post, checkpoint53d_post, checkpoint46d_pre, checkpoint60, checkpoint57a_post, checkpoint48d_post, checkpoint57h_pre, checkpoint48f_post, checkpoint52b_pre, checkpoint54b_post, checkpoint46j_pre, checkpoint58w_post, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint57y_pre, checkpoint55g_post, checkpoint48h_post, checkpoint51q_post, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint52b_post, checkpoint52c_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, checkpoint58o_post, checkpoint48a_post, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint58p_post, checkpoint58q_post, checkpoint52f_pre, checkpoint55d_post, checkpoint58e_post, checkpoint47j_post, checkpoint54a_pre, checkpoint53c_post, checkpoint55d_pre, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, branch-exfmods-tag, branchpoint-genmake2, checkpoint54a_post, checkpoint46e_pre, checkpoint55h_post, checkpoint58n_post, checkpoint51r_post, checkpoint48c_post, checkpoint46b_pre, checkpoint51i_post, checkpoint57e_post, checkpoint55b_post, checkpoint51b_post, checkpoint51c_post, checkpoint46c_pre, checkpoint53a_post, checkpoint46, checkpoint47b_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint46h_pre, checkpoint52d_post, checkpoint53g_post, checkpoint46m_post, checkpoint57p_post, checkpint57u_post, checkpoint46a_pre, checkpoint50g_post, checkpoint57q_post, eckpoint57e_pre, checkpoint46g_post, checkpoint58k_post, checkpoint52a_pre, checkpoint58v_post, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint56a_post, checkpoint58l_post, checkpoint53f_post, checkpoint57h_done, checkpoint52j_post, checkpoint47f_post, checkpoint50e_post, checkpoint46i_post, checkpoint57j_post, checkpoint57f_pre, checkpoint46c_post, checkpoint58g_post, branch-netcdf, checkpoint50d_pre, checkpoint58x_post, checkpoint52n_post, checkpoint53b_pre, checkpoint46e_post, checkpoint59j, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint51e_post, checkpoint57a_pre, checkpoint55a_post, checkpoint47, checkpoint48, checkpoint49, checkpoint57o_post, checkpoint46h_post, checkpoint51o_post, checkpoint57k_post, checkpoint51f_pre, checkpoint48g_post, checkpoint53b_post, checkpoint47h_post, checkpoint52a_post, checkpoint57w_post, checkpoint58i_post, checkpoint51g_post, ecco_c52_e35, checkpoint57x_post, checkpoint46d_post, checkpoint50b_post, checkpoint58c_post, checkpoint58u_post, checkpoint51m_post, checkpoint53d_pre, checkpoint58s_post, checkpoint55e_post, checkpoint54c_post, checkpoint51a_post, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-exfmods-curt, branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.1: +191 -0 lines 
Merging from release1_p5
o flow directives:
  - new for S/R do_the_model_io
  - enhanced for S/R checkpoint
  - new gmredi.flow
o Introduced CPP option ALLOW_AUTODIFF_MONITOR to
  disable adcommon.h/g_common.h by default.
o fixed problem for adjoint of global_max
1 heimbach 1.2 C $Header: /u/gcmpack/MITgcm/pkg/autodiff/Attic/global_max_ad.F,v 1.1.2.1 2002/04/17 01:35:52 heimbach Exp $
2     C $Name: release1_p5 $
3    
4     C-- File global_max.F: Routines that perform global max reduction on an array
5     C of thread values.
6     C
7     C Adjoint routines are identical to forward routines,
8     C but parameter list is reversed to be consistent with
9     C call statement generated by TAMC/TAF
10     C P. Heimbach, 16-Apr-2002
11     C
12     C Contents
13     C o global_admax_r4
14     C o global_admax_r8
15     #include "CPP_EEOPTIONS.h"
16    
17     CBOP
18    
19     C !ROUTINE: GLOBAL_ADMAX_R4
20    
21     C !INTERFACE:
22     SUBROUTINE GLOBAL_ADMAX_R4(
23     I myThid,
24     U maxPhi
25     & )
26     IMPLICIT NONE
27     C !DESCRIPTION:
28     C *==========================================================*
29     C | SUBROUTINE GLOBAL_ADMAX_R4
30     C | o Handle max for real*4 data.
31     C *==========================================================*
32     C | Perform max on array of one value per thread and then
33     C | max result of all the processes.
34     C | Notes
35     C | =====
36     C | Within a process only one thread does the max, each
37     C | thread is assumed to have already maxed its local data.
38     C | The same thread also does the inter-process max for
39     C | example with MPI and then writes the result into a shared
40     C | location. All threads wait until the max is avaiailable.
41     C *==========================================================*
42    
43     C !USES:
44     C == Global data ==
45     #include "SIZE.h"
46     #include "EEPARAMS.h"
47     #include "EESUPPORT.h"
48     #include "GLOBAL_MAX.h"
49    
50     C !INPUT/OUTPUT PARAMETERS:
51     C == Routine arguments ==
52     C maxPhi :: Result of max.
53     C myThid :: My thread id.
54     Real*4 maxPhi
55     INTEGER myThid
56    
57     C !LOCAL VARIABLES:
58     C == Local variables ==
59     C I :: Loop counters
60     C mpiRC :: MPI return code
61     INTEGER I
62     Real*4 tmp
63     #ifdef ALLOW_USE_MPI
64     INTEGER mpiRC
65     #endif /* ALLOW_USE_MPI */
66     CEOP
67    
68     CALL BAR2( myThid )
69     C-- write local max into array
70     phiGMRS(1,myThid) = maxPhi
71    
72     C-- Can not start until everyone is ready
73     CALL BAR2( myThid )
74    
75     C-- Max within the process first
76     _BEGIN_MASTER( myThid )
77     tmp = phiGMRS(1,1)
78     DO I=2,nThreads
79     tmp = MAX(tmp,phiGMRS(1,I))
80     ENDDO
81     maxPhi = tmp
82     #ifdef ALLOW_USE_MPI
83     #ifndef ALWAYS_USE_MPI
84     IF ( usingMPI ) THEN
85     #endif
86     CALL MPI_Allreduce(tmp,maxPhi,1,MPI_REAL,MPI_MAX,
87     & MPI_COMM_WORLD,mpiRC)
88     #ifndef ALWAYS_USE_MPI
89     ENDIF
90     #endif
91     #endif /* ALLOW_USE_MPI */
92     phiGMRS(1,1) = maxPhi
93     _END_MASTER( myThid )
94     C--
95     CALL BAR2( myThid )
96    
97     C-- set result for every process
98     maxPhi = phiGMRS(1,1)
99     CALL BAR2( myThid )
100    
101     RETURN
102     END
103    
104     CBOP
105    
106     C !ROUTINE: GLOBAL_ADMAX_R8
107    
108     C !INTERFACE:
109     SUBROUTINE GLOBAL_ADMAX_R8(
110     I myThid,
111     O maxPhi
112     & )
113     IMPLICIT NONE
114     C !DESCRIPTION:
115     C *==========================================================*
116     C | SUBROUTINE GLOBAL_ADMAX_R8
117     C | o Handle max for real*8 data.
118     C *==========================================================*
119     C | Perform max on array of one value per thread and then
120     C | max result of all the processes.
121     C | Notes
122     C | =====
123     C | Within a process only one thread does the max, each
124     C | thread is assumed to have already maxed its local data.
125     C | The same thread also does the inter-process max for
126     C | example with MPI and then writes the result into a shared
127     C | location. All threads wait until the max is avaiailable.
128     C *==========================================================*
129    
130     C !USES:
131     C === Global data ===
132     #include "SIZE.h"
133     #include "EEPARAMS.h"
134     #include "EESUPPORT.h"
135     #include "GLOBAL_MAX.h"
136    
137     C !INPUT/OUTPUT PARAMETERS:
138     C === Routine arguments ===
139     C maxPhi :: Result of max.
140     C myThid :: My thread id.
141     Real*8 maxPhi
142     INTEGER myThid
143    
144     C !LOCAL VARIABLES:
145     C === Local variables ===
146     C I :: Loop counters
147     C mpiRC :: MPI return code
148     INTEGER I
149     Real*8 tmp
150     #ifdef ALLOW_USE_MPI
151     INTEGER mpiRC
152     #endif /* ALLOW_USE_MPI */
153     CEOP
154    
155     CALL BAR2( myThid )
156     C-- write local max into array
157     phiGMRL(1,myThid) = maxPhi
158    
159     C-- Can not start until everyone is ready
160     CALL BAR2( myThid )
161    
162     C-- Max within the process first
163     _BEGIN_MASTER( myThid )
164     tmp = phiGMRL(1,1)
165     DO I=2,nThreads
166     tmp = MAX(tmp,phiGMRL(1,I))
167     ENDDO
168     maxPhi = tmp
169     #ifdef ALLOW_USE_MPI
170     #ifndef ALWAYS_USE_MPI
171     IF ( usingMPI ) THEN
172     #endif
173     CALL MPI_Allreduce(tmp,maxPhi,1,MPI_DOUBLE_PRECISION,MPI_MAX,
174     & MPI_COMM_WORLD,mpiRC)
175     #ifndef ALWAYS_USE_MPI
176     ENDIF
177     #endif
178     #endif /* ALLOW_USE_MPI */
179     C-- Write solution to place where all threads can see it
180     phiGMRL(1,1) = maxPhi
181     _END_MASTER( myThid )
182    
183     C-- Do not leave until we are sure that the max is done
184     CALL BAR2( myThid )
185    
186     C-- set result for every process
187     maxPhi = phiGMRL(1,1)
188     CALL BAR2( myThid )
189    
190     RETURN
191     END
  ViewVC Help
Powered by ViewVC 1.1.22