source: soft/giet_vm/applications/ocean/subblock.C @ 588

Last change on this file since 588 was 581, checked in by laurent, 9 years ago

Adding ocean application, some mathematics functions and distributed locks

File size: 3.9 KB
RevLine 
[581]1/*************************************************************************/
2/*                                                                       */
3/*  Copyright (c) 1994 Stanford University                               */
4/*                                                                       */
5/*  All rights reserved.                                                 */
6/*                                                                       */
7/*  Permission is given to use, copy, and modify this software for any   */
8/*  non-commercial purpose as long as this copyright notice is not       */
9/*  removed.  All other uses, including redistribution in whole or in    */
10/*  part, are forbidden without prior written permission.                */
11/*                                                                       */
12/*  This software is provided with absolutely no warranty and no         */
13/*  support.                                                             */
14/*                                                                       */
15/*************************************************************************/
16
17EXTERN_ENV
18
19#include <stdio.h>
20#include <math.h>
21
22#include "decs.h"
23
24
25void subblock()
26{
27    long i;
28    long j;
29    long k;
30    long xportion;
31    long yportion;
32    long my_num;
33
34/* Determine starting coord and number of points to process in     */
35/* each direction                                                  */
36
37    for (i = 0; i < numlev; i++) {
38        xportion = (jmx[i] - 2) / xprocs;
39        //xextra = (jmx[i] - 2) % xprocs;
40        for (j = 0; j < xprocs; j++) {
41            for (k = 0; k < yprocs; k++) {
42                gp[k * xprocs + j].rel_num_x[i] = xportion;
43            }
44        }
45        yportion = (imx[i] - 2) / yprocs;
46        //yextra = (imx[i] - 2) % yprocs;
47        for (j = 0; j < yprocs; j++) {
48            for (k = 0; k < xprocs; k++) {
49                gp[j * xprocs + k].rel_num_y[i] = yportion;
50            }
51        }
52    }
53
54    for (my_num = 0; my_num < nprocs; my_num++) {
55        for (i = 0; i < numlev; i++) {
56            gp[my_num].rlist[i] = 1;
57            gp[my_num].rljst[i] = 1;
58            gp[my_num].rlien[i] = gp[my_num].rlist[i] + gp[my_num].rel_num_y[i];
59            gp[my_num].rljen[i] = gp[my_num].rljst[i] + gp[my_num].rel_num_x[i];
60            gp[my_num].eist[i] = gp[my_num].rlist[i] + 1;
61            gp[my_num].oist[i] = gp[my_num].rlist[i];
62            gp[my_num].ejst[i] = gp[my_num].rljst[i] + 1;
63            gp[my_num].ojst[i] = gp[my_num].rljst[i];
64        }
65    }
66   
67    for (i = 0; i < nprocs; i++) {
68        gp[i].neighbors[LEFT] = -1;
69        gp[i].neighbors[RIGHT] = -1;
70        gp[i].neighbors[UP] = -1;
71        gp[i].neighbors[DOWN] = -1;
72        gp[i].neighbors[UPLEFT] = -1;
73        gp[i].neighbors[UPRIGHT] = -1;
74        gp[i].neighbors[DOWNLEFT] = -1;
75        gp[i].neighbors[DOWNRIGHT] = -1;
76       
77        if (i >= xprocs) {
78            gp[i].neighbors[UP] = i - xprocs;
79        }
80        if (i < nprocs - xprocs) {
81            gp[i].neighbors[DOWN] = i + xprocs;
82        }
83        if ((i % xprocs) > 0) {
84            gp[i].neighbors[LEFT] = i - 1;
85        }
86        if ((i % xprocs) < (xprocs - 1)) {
87            gp[i].neighbors[RIGHT] = i + 1;
88        }
89       
90        j = gp[i].neighbors[UP];
91       
92        if (j != -1) {
93            if ((j % xprocs) > 0) {
94                gp[i].neighbors[UPLEFT] = j - 1;
95            }
96            if ((j % xprocs) < (xprocs - 1)) {
97                gp[i].neighbors[UPRIGHT] = j + 1;
98            }
99        }
100       
101        j = gp[i].neighbors[DOWN];
102       
103        if (j != -1) {
104            if ((j % xprocs) > 0) {
105                gp[i].neighbors[DOWNLEFT] = j - 1;
106            }
107            if ((j % xprocs) < (xprocs - 1)) {
108                gp[i].neighbors[DOWNRIGHT] = j + 1;
109            }
110        }
111    }
112   
113    for (i = 0; i < nprocs; i++) {
114        (*gp[i].rownum) = i / xprocs;
115        (*gp[i].colnum) = i % xprocs;
116    }
117}
Note: See TracBrowser for help on using the repository browser.