← Back to branch summary

~ubuntu-branches/ubuntu/precise/code-saturne/precise

« back to all changes in this revision

Viewing changes to tests/fvm_interface_test.c

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2011-11-24 00:00:08 UTC
  • mfrom: (6.1.9 sid)
  • Revision ID: package-import@ubuntu.com-20111124000008-2vo99e38267942q5
Tags: 2.1.0-3
Install a missing file

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/fvm_interface_test.c
 
1
/*============================================================================
 
2
 * Unit test for fvm_interface.c;
 
3
 *============================================================================*/
 
4
 
 
5
/*
 
6
  This file is part of Code_Saturne, a general-purpose CFD tool.
 
7
 
 
8
  Copyright (C) 1998-2011 EDF S.A.
 
9
 
 
10
  This program is free software; you can redistribute it and/or modify it under
 
11
  the terms of the GNU General Public License as published by the Free Software
 
12
  Foundation; either version 2 of the License, or (at your option) any later
 
13
  version.
 
14
 
 
15
  This program is distributed in the hope that it will be useful, but WITHOUT
 
16
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
 
17
  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
 
18
  details.
 
19
 
 
20
  You should have received a copy of the GNU General Public License along with
 
21
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
 
22
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
 
23
*/
 
24
 
 
25
/*----------------------------------------------------------------------------*/
 
26
 
 
27
#if defined(HAVE_CONFIG_H)
 
28
#include "cs_config.h"
 
29
#endif
 
30
 
 
31
#include <assert.h>
 
32
#include <stdarg.h>
 
33
#include <stdio.h>
 
34
#include <stdlib.h>
 
35
#include <string.h>
 
36
 
 
37
#include <bft_error.h>
 
38
#include <bft_mem.h>
 
39
#include <bft_printf.h>
 
40
 
 
41
#include "fvm_interface.h"
 
42
#include "fvm_parall.h"
 
43
#include "fvm_periodicity.h"
 
44
 
 
45
/*---------------------------------------------------------------------------*/
 
46
 
 
47
static fvm_periodicity_t *
 
48
_init_periodicity(void)
 
49
{
 
50
  double t1[3] = {1., 0., 0.};
 
51
  double t2[3] = {0., 1., 0.};
 
52
  double t3[3] = {0., 0., 1.};
 
53
 
 
54
  /* double t4[3] = {1.02, 0., 0.}; */
 
55
  /* double a1[3] = {0., 0., 1.}; */
 
56
  /* double i1[3] = {0., 0., 0.}; */
 
57
 
 
58
  fvm_periodicity_t *p = NULL;
 
59
 
 
60
  p = fvm_periodicity_create(0.1);
 
61
 
 
62
  fvm_periodicity_add_translation(p, 1, t1);
 
63
  fvm_periodicity_add_translation(p, 2, t2);
 
64
  fvm_periodicity_add_translation(p, 3, t3);
 
65
 
 
66
  /* fvm_periodicity_add_rotation(p, 2, 90., a1, i1); */
 
67
  /* fvm_periodicity_add_translation(p, 3, t3); */
 
68
 
 
69
  fvm_periodicity_combine(p, 1);
 
70
 
 
71
  return p;
 
72
}
 
73
 
 
74
static fvm_interface_set_t *
 
75
_periodic_is(int                       comm_size,
 
76
             int                       comm_rank,
 
77
             fvm_gnum_t                n_side,
 
78
             int                       n_periodic_lists,
 
79
             const fvm_periodicity_t  *perio)
 
80
{
 
81
  fvm_gnum_t ii, jj, kk, couple_id;
 
82
 
 
83
  fvm_gnum_t n_max = n_side * n_side * n_side;
 
84
  fvm_lnum_t n_elements = ((n_max+1)*4) / (3*comm_size);
 
85
  int periodicity_num[3] = {1, 2, 3};
 
86
 
 
87
  fvm_lnum_t *n_periodic_couples = NULL;
 
88
  fvm_gnum_t **periodic_couples = NULL;
 
89
  fvm_gnum_t *global_number = NULL;
 
90
 
 
91
  fvm_interface_set_t *ifset = NULL;
 
92
 
 
93
  if (comm_size > 1) {
 
94
 
 
95
    BFT_MALLOC(global_number, n_elements, fvm_gnum_t);
 
96
 
 
97
    for (ii = 0; ii < (fvm_gnum_t)n_elements; ii++) {
 
98
      global_number[ii] = (n_elements * 3 / 4) * comm_rank + ii + 1;
 
99
    }
 
100
    if (comm_rank == comm_size -1 && global_number[n_elements - 1] < n_max)
 
101
      global_number[n_elements - 1] = n_max;
 
102
 
 
103
  }
 
104
 
 
105
  BFT_MALLOC(n_periodic_couples, n_periodic_lists, fvm_lnum_t);
 
106
  BFT_MALLOC(periodic_couples, n_periodic_lists, fvm_gnum_t *);
 
107
 
 
108
  for (ii = 0; ii < (fvm_gnum_t)n_periodic_lists; ii++) {
 
109
 
 
110
    n_periodic_couples[ii] = 0;
 
111
    periodic_couples[ii] = NULL;
 
112
 
 
113
    if (comm_rank == 0 || comm_rank == 1) {
 
114
 
 
115
      n_periodic_couples[ii] = n_side*n_side;
 
116
      BFT_MALLOC(periodic_couples[ii], n_periodic_couples[ii]*2, fvm_gnum_t);
 
117
 
 
118
      couple_id = 0;
 
119
 
 
120
      switch(ii) {
 
121
      case 0:
 
122
        for (jj = 0; jj < n_side*n_side; jj++) {
 
123
          periodic_couples[ii][couple_id*2] = (jj*n_side) + 1;
 
124
          periodic_couples[ii][couple_id*2+1]
 
125
            = periodic_couples[ii][couple_id*2] + (n_side-1);
 
126
          couple_id++;
 
127
        }
 
128
        break;
 
129
      case 1:
 
130
        for (jj = 0; jj < n_side; jj++) {
 
131
          for (kk = 0; kk < n_side; kk++) {
 
132
            periodic_couples[ii][couple_id*2] = jj*n_side*n_side + kk + 1;
 
133
            periodic_couples[ii][couple_id*2+1]
 
134
              = periodic_couples[ii][couple_id*2] + (n_side-1)*n_side;
 
135
            couple_id++;
 
136
          }
 
137
        }
 
138
        break;
 
139
      case 2:
 
140
        for (jj = 0; jj < n_side*n_side; jj++) {
 
141
          periodic_couples[ii][couple_id*2] = jj + 1;
 
142
          periodic_couples[ii][couple_id*2+1]
 
143
            = periodic_couples[ii][couple_id*2] + (n_side-1)*n_side*n_side;
 
144
          couple_id++;
 
145
        }
 
146
        break;
 
147
      default:
 
148
        break;
 
149
      }
 
150
 
 
151
      n_periodic_couples[ii] = couple_id;
 
152
      BFT_REALLOC(periodic_couples[ii], n_periodic_couples[ii]*2, fvm_gnum_t);
 
153
 
 
154
    }
 
155
 
 
156
  }
 
157
 
 
158
  ifset = fvm_interface_set_create(n_elements,
 
159
                                   NULL,
 
160
                                   global_number,
 
161
                                   perio,
 
162
                                   n_periodic_lists,
 
163
                                   periodicity_num,
 
164
                                   n_periodic_couples,
 
165
                                   (const fvm_gnum_t **const)periodic_couples);
 
166
 
 
167
  for (ii = 0; ii < (fvm_gnum_t)n_periodic_lists; ii++)
 
168
    BFT_FREE(periodic_couples[ii]);
 
169
  BFT_FREE(periodic_couples);
 
170
  BFT_FREE(n_periodic_couples);
 
171
 
 
172
  if (comm_size > 1)
 
173
    BFT_FREE(global_number);
 
174
 
 
175
  return ifset;
 
176
}
 
177
 
 
178
/*----------------------------------------------------------------------------
 
179
 * Fonction d'impression d'un message sur la sortie standard
 
180
 *----------------------------------------------------------------------------*/
 
181
 
 
182
static int _bft_printf_proxy
 
183
(
 
184
 const char     *const format,
 
185
       va_list         arg_ptr
 
186
)
 
187
{
 
188
  static FILE *f = NULL;
 
189
 
 
190
  if (f == NULL) {
 
191
    char filename[64];
 
192
    int rank = 0;
 
193
#if defined(HAVE_MPI)
 
194
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
 
195
#endif
 
196
    sprintf (filename, "fvm_interface_test_out.%d", rank);
 
197
    f = fopen(filename, "w");
 
198
    assert(f != NULL);
 
199
  }
 
200
 
 
201
  return vfprintf(f, format, arg_ptr);
 
202
}
 
203
 
 
204
/*---------------------------------------------------------------------------*/
 
205
 
 
206
int
 
207
main (int argc, char *argv[])
 
208
{
 
209
  char mem_trace_name[32];
 
210
  int size = 1;
 
211
  int rank = 0;
 
212
  fvm_interface_set_t *ifset = NULL;
 
213
  fvm_periodicity_t *perio = NULL;
 
214
 
 
215
#if defined(HAVE_MPI)
 
216
 
 
217
  int ii;
 
218
 
 
219
  MPI_Status status;
 
220
 
 
221
  int sync = 1;
 
222
 
 
223
  fvm_lnum_t n_elements = 20;
 
224
  fvm_gnum_t *global_number = NULL;
 
225
 
 
226
  /* Initialization */
 
227
 
 
228
  MPI_Init(&argc, &argv);
 
229
 
 
230
  fvm_parall_set_mpi_comm(MPI_COMM_WORLD);
 
231
 
 
232
  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
 
233
  MPI_Comm_size(MPI_COMM_WORLD, &size);
 
234
 
 
235
  if (size > 1)
 
236
    sprintf(mem_trace_name, "fvm_interface_test_mem.%d", rank);
 
237
  else
 
238
    strcpy(mem_trace_name, "fvm_interface_test_mem");
 
239
  bft_mem_init(mem_trace_name);
 
240
 
 
241
  bft_printf_proxy_set(_bft_printf_proxy);
 
242
 
 
243
  /* Build arbitrary interface */
 
244
 
 
245
  BFT_MALLOC(global_number, n_elements, fvm_gnum_t);
 
246
 
 
247
  for (ii = 0; ii < n_elements; ii++) {
 
248
    global_number[ii] = (n_elements * 3 / 4) * rank + ii + 1;
 
249
  }
 
250
 
 
251
  ifset = fvm_interface_set_create(n_elements,
 
252
                                   NULL,
 
253
                                   global_number,
 
254
                                   NULL,
 
255
                                   0,
 
256
                                   NULL,
 
257
                                   NULL,
 
258
                                   NULL);
 
259
 
 
260
  BFT_FREE(global_number);
 
261
 
 
262
  /* Serialize dump of interfaces */
 
263
 
 
264
  if (rank > 0)
 
265
    MPI_Recv(&sync, 1, MPI_INT, rank - 1, 0, MPI_COMM_WORLD, &status);
 
266
 
 
267
  bft_printf("Interface on rank %d:\n\n", rank);
 
268
 
 
269
  fvm_interface_set_dump(ifset);
 
270
 
 
271
  if (rank < size - 1)
 
272
    MPI_Send(&sync, 1, MPI_INT, rank + 1, 0, MPI_COMM_WORLD);
 
273
 
 
274
  /* We are finished for this interface */
 
275
 
 
276
  ifset = fvm_interface_set_destroy(ifset);
 
277
 
 
278
#else
 
279
 
 
280
  strcpy(mem_trace_name, "fvm_interface_test_mem");
 
281
  bft_mem_init(mem_trace_name);
 
282
 
 
283
  bft_printf_proxy_set(_bft_printf_proxy);
 
284
 
 
285
#endif /* (HAVE_MPI) */
 
286
 
 
287
  /* Now build interface with periodicity */
 
288
 
 
289
  perio = _init_periodicity();
 
290
 
 
291
  /* Dump periodicity info for first rank (identical on all ranks) */
 
292
 
 
293
  if (rank == 0)
 
294
    fvm_periodicity_dump(perio);
 
295
 
 
296
  ifset = _periodic_is(size,
 
297
                       rank,
 
298
                       4, /* n vertices/cube side */
 
299
                       3, /* n_periodic_lists */
 
300
                       perio);
 
301
 
 
302
#if defined(HAVE_MPI)
 
303
 
 
304
  /* Serialize dump of interfaces */
 
305
 
 
306
  if (rank > 0)
 
307
    MPI_Recv(&sync, 1, MPI_INT, rank - 1, 0, MPI_COMM_WORLD, &status);
 
308
 
 
309
#endif /* (HAVE_MPI) */
 
310
 
 
311
  bft_printf("Periodic Interface on rank %d:\n\n", rank);
 
312
 
 
313
  fvm_interface_set_dump(ifset);
 
314
 
 
315
#if defined(HAVE_MPI)
 
316
 
 
317
  if (rank < size - 1)
 
318
    MPI_Send(&sync, 1, MPI_INT, rank + 1, 0, MPI_COMM_WORLD);
 
319
 
 
320
#endif /* (HAVE_MPI) */
 
321
 
 
322
  /* We are finished */
 
323
 
 
324
  ifset = fvm_interface_set_destroy(ifset);
 
325
  fvm_periodicity_destroy(perio);
 
326
 
 
327
  bft_mem_end();
 
328
 
 
329
#if defined(HAVE_MPI)
 
330
 
 
331
  MPI_Finalize();
 
332
 
 
333
#endif
 
334
 
 
335
  exit (EXIT_SUCCESS);
 
336
}