src/gmxlib/nonbonded/nb_kernel_bluegene/nb_kernel_bluegene.c

   1 /* -*- mode: c; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4; c-file-style: "stroustrup"; -*-
   2  *
   3  *
   4  * This file is part of Gromacs        Copyright (c) 1991-2004
   5  * David van der Spoel, Erik Lindahl, University of Groningen.
   6  *
   7  * This program is free software; you can redistribute it and/or
   8  * modify it under the terms of the GNU General Public License
   9  * as published by the Free Software Foundation; either version 2
  10  * of the License, or (at your option) any later version.
  11  *
  12  * To help us fund GROMACS development, we humbly ask that you cite
  13  * the research papers on the package. Check out http://www.gromacs.org
  14  *
  15  * And Hey:
  16  * Gnomes, ROck Monsters And Chili Sauce
  17  */
  18 #ifdef HAVE_CONFIG_H
  19 #include <config.h>
  20 #endif
  21
  22 /* Must come directly after config.h */
  23 #ifdef GMX_THREAD_SHM_FDECOMP
  24 #include <thread_mpi.h>
  25 #endif
  26
  27 #include <types/simple.h>
  28 #include <types/nrnb.h>
  29
  30 #include "nb_kernel_bluegene.h"
  31
  32 /* Include Bluegene kernel headers in local directory */
  33 #include "nb_kernel010_bluegene.h"
  34 #include "nb_kernel020_bluegene.h"
  35 #include "nb_kernel030_bluegene.h"
  36 #include "nb_kernel100_bluegene.h"
  37 #include "nb_kernel101_bluegene.h"
  38 #include "nb_kernel102_bluegene.h"
  39 #include "nb_kernel103_bluegene.h"
  40 #include "nb_kernel104_bluegene.h"
  41 #include "nb_kernel110_bluegene.h"
  42 #include "nb_kernel111_bluegene.h"
  43 #include "nb_kernel112_bluegene.h"
  44 #include "nb_kernel113_bluegene.h"
  45 #include "nb_kernel114_bluegene.h"
  46 #include "nb_kernel120_bluegene.h"
  47 #include "nb_kernel121_bluegene.h"
  48 #include "nb_kernel122_bluegene.h"
  49 #include "nb_kernel123_bluegene.h"
  50 #include "nb_kernel124_bluegene.h"
  51 #include "nb_kernel130_bluegene.h"
  52 #include "nb_kernel131_bluegene.h"
  53 #include "nb_kernel132_bluegene.h"
  54 #include "nb_kernel133_bluegene.h"
  55 #include "nb_kernel134_bluegene.h"
  56 #include "nb_kernel200_bluegene.h"
  57 #include "nb_kernel201_bluegene.h"
  58 #include "nb_kernel202_bluegene.h"
  59 #include "nb_kernel203_bluegene.h"
  60 #include "nb_kernel204_bluegene.h"
  61 #include "nb_kernel210_bluegene.h"
  62 #include "nb_kernel211_bluegene.h"
  63 #include "nb_kernel212_bluegene.h"
  64 #include "nb_kernel213_bluegene.h"
  65 #include "nb_kernel214_bluegene.h"
  66 #include "nb_kernel220_bluegene.h"
  67 #include "nb_kernel221_bluegene.h"
  68 #include "nb_kernel222_bluegene.h"
  69 #include "nb_kernel223_bluegene.h"
  70 #include "nb_kernel224_bluegene.h"
  71 #include "nb_kernel230_bluegene.h"
  72 #include "nb_kernel231_bluegene.h"
  73 #include "nb_kernel232_bluegene.h"
  74 #include "nb_kernel233_bluegene.h"
  75 #include "nb_kernel234_bluegene.h"
  76 #include "nb_kernel300_bluegene.h"
  77 #include "nb_kernel301_bluegene.h"
  78 #include "nb_kernel302_bluegene.h"
  79 #include "nb_kernel303_bluegene.h"
  80 #include "nb_kernel304_bluegene.h"
  81 #include "nb_kernel310_bluegene.h"
  82 #include "nb_kernel311_bluegene.h"
  83 #include "nb_kernel312_bluegene.h"
  84 #include "nb_kernel313_bluegene.h"
  85 #include "nb_kernel314_bluegene.h"
  86 #include "nb_kernel320_bluegene.h"
  87 #include "nb_kernel321_bluegene.h"
  88 #include "nb_kernel322_bluegene.h"
  89 #include "nb_kernel323_bluegene.h"
  90 #include "nb_kernel324_bluegene.h"
  91 #include "nb_kernel330_bluegene.h"
  92 #include "nb_kernel331_bluegene.h"
  93 #include "nb_kernel332_bluegene.h"
  94 #include "nb_kernel333_bluegene.h"
  95 #include "nb_kernel334_bluegene.h"
  96 #include "nb_kernel400_bluegene.h"
  97 #include "nb_kernel410_bluegene.h"
  98 #include "nb_kernel420_bluegene.h"
  99 #include "nb_kernel430_bluegene.h"
 100
 101
 102
 103 #include <stdlib.h>
 104 #include <stdio.h>
 105 /* Necessary headers for POSIX-style long jumps. */
 106 #include <signal.h>
 107 #include <setjmp.h>
 108
 109
 110
 111 #include "../nb_kerneltype.h"
 112 #include "nb_kernel_bluegene.h"
 113
 114
 115 static nb_kernel_t *
 116 kernellist_bluegene[eNR_NBKERNEL_NR] =
 117 {
 118     nb_kernel010_bluegene,
 119     nb_kernel020_bluegene,
 120     nb_kernel030_bluegene,
 121     nb_kernel100_bluegene,
 122     nb_kernel101_bluegene,
 123     nb_kernel102_bluegene,
 124     nb_kernel103_bluegene,
 125     nb_kernel104_bluegene,
 126     nb_kernel110_bluegene,
 127     nb_kernel111_bluegene,
 128     nb_kernel112_bluegene,
 129     nb_kernel113_bluegene,
 130     nb_kernel114_bluegene,
 131     nb_kernel120_bluegene,
 132     nb_kernel121_bluegene,
 133     nb_kernel122_bluegene,
 134     nb_kernel123_bluegene,
 135     nb_kernel124_bluegene,
 136     nb_kernel130_bluegene,
 137     nb_kernel131_bluegene,
 138     nb_kernel132_bluegene,
 139     nb_kernel133_bluegene,
 140     nb_kernel134_bluegene,
 141     nb_kernel200_bluegene,
 142     nb_kernel201_bluegene,
 143     nb_kernel202_bluegene,
 144     nb_kernel203_bluegene,
 145     nb_kernel204_bluegene,
 146     nb_kernel210_bluegene,
 147     nb_kernel211_bluegene,
 148     nb_kernel212_bluegene,
 149     nb_kernel213_bluegene,
 150     nb_kernel214_bluegene,
 151     nb_kernel220_bluegene,
 152     nb_kernel221_bluegene,
 153     nb_kernel222_bluegene,
 154     nb_kernel223_bluegene,
 155     nb_kernel224_bluegene,
 156     nb_kernel230_bluegene,
 157     nb_kernel231_bluegene,
 158     nb_kernel232_bluegene,
 159     nb_kernel233_bluegene,
 160     nb_kernel234_bluegene,
 161     nb_kernel300_bluegene,
 162     nb_kernel301_bluegene,
 163     nb_kernel302_bluegene,
 164     nb_kernel303_bluegene,
 165     nb_kernel304_bluegene,
 166     nb_kernel310_bluegene,
 167     nb_kernel311_bluegene,
 168     nb_kernel312_bluegene,
 169     nb_kernel313_bluegene,
 170     nb_kernel314_bluegene,
 171     nb_kernel320_bluegene,
 172     nb_kernel321_bluegene,
 173     nb_kernel322_bluegene,
 174     nb_kernel323_bluegene,
 175     nb_kernel324_bluegene,
 176     nb_kernel330_bluegene,
 177     nb_kernel331_bluegene,
 178     nb_kernel332_bluegene,
 179     nb_kernel333_bluegene,
 180     nb_kernel334_bluegene,
 181     nb_kernel400_bluegene,
 182     nb_kernel410_bluegene,
 183     nb_kernel430_bluegene
 184 };
 185
 186
 187 #ifdef GMX_THREAD_SHM_FDECOMP
 188 static tMPI_Thread_mutex_t
 189 nb_kernel_bluegene_test_mutex = TMPI_THREAD_MUTEX_INITIALIZER;
 190 #endif
 191
 192
 193 /*! Posix long jump label */
 194 static jmp_buf
 195 nb_kernel_bluegene_testprog;
 196
 197 /*! Result of bluegene test */
 198 static gmx_bool
 199 nb_kernel_bluegene_present;
 200
 201
 202 static void
 203 nb_kernel_bluegene_sigill_handler(int n)
 204 {
 205   nb_kernel_bluegene_present=FALSE;
 206   longjmp(nb_kernel_bluegene_testprog,n);
 207 }
 208
 209
 210
 211 /* Return 0 if Bluegene support is present, or
 212  * non-zero on failure.
 213  */
 214 int
 215 nb_kernel_bluegene_test(FILE *                log)
 216 {
 217         /*
 218          * This should NOT be called from threads,
 219          * but just in case you still try to do it...
 220          */
 221 #ifdef GMX_THREAD_SHM_FDECOMP
 222         tMPI_Thread_mutex_lock(&nb_kernel_bluegene_test_mutex);
 223 #endif
 224
 225     if(log)
 226         fprintf(log,"Testing BlueGene support...");
 227
 228         nb_kernel_bluegene_present = TRUE;
 229         signal(SIGILL,nb_kernel_bluegene_sigill_handler);
 230
 231         /* return to this point after executing the signal handler
 232          * if we catch a SIGILL
 233          */
 234         setjmp(nb_kernel_bluegene_testprog);
 235
 236         if(nb_kernel_bluegene_present)
 237     {
 238         double _Complex a,b;
 239         double c = 2.0;
 240
 241         a = __cmplx(c,c);
 242         b = __fpmadd(a,a,a);
 243     }
 244
 245         /* If Bluegene worked, then success is still 1.
 246      * If we got SIGILL, it was set to 0 in sigill_handler().
 247      */
 248
 249         if(log)
 250                 fprintf(log," %spresent.\n",
 251                                 nb_kernel_bluegene_present ? "":"not ");
 252
 253 #ifdef GMX_THREAD_SHM_FDECOMP
 254         tMPI_Thread_mutex_unlock(&nb_kernel_bluegene_test_mutex);
 255 #endif
 256
 257         return ((nb_kernel_bluegene_present) ? 0 : -1);
 258 }
 259
 260
 261
 262
 263 void
 264 nb_kernel_setup_bluegene(FILE *log,nb_kernel_t **list)
 265 {
 266     int i;
 267     nb_kernel_t *p;
 268
 269     if(nb_kernel_bluegene_test(log) != 0)
 270         return;
 271
 272     for(i=0;i<eNR_NBKERNEL_NR;i++)
 273     {
 274         p = kernellist_bluegene[i];
 275         if(p!=NULL)
 276             list[i] = p;
 277     }
 278 }
 279
 280
 281
 282
 283
 284
 285
 286
 287
 288