src/gmxlib/nonbonded/nb_kernel_ppc_altivec/nb_kernel_ppc_altivec.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
  23 #include <types/nrnb.h>
  24 #include "types/simple.h"
  25
  26 #include "../nb_kerneltype.h"
  27 #include "nb_kernel_ppc_altivec.h"
  28 #include "nb_kernel_ppc_altivec_test.h"
  29
  30 /* Include altivec kernel headers in local directory */
  31 #include "nb_kernel010_ppc_altivec.h"
  32 #include "nb_kernel030_ppc_altivec.h"
  33 #include "nb_kernel100_ppc_altivec.h"
  34 #include "nb_kernel101_ppc_altivec.h"
  35 #include "nb_kernel102_ppc_altivec.h"
  36 #include "nb_kernel103_ppc_altivec.h"
  37 #include "nb_kernel104_ppc_altivec.h"
  38 #include "nb_kernel110_ppc_altivec.h"
  39 #include "nb_kernel111_ppc_altivec.h"
  40 #include "nb_kernel112_ppc_altivec.h"
  41 #include "nb_kernel113_ppc_altivec.h"
  42 #include "nb_kernel114_ppc_altivec.h"
  43 #include "nb_kernel130_ppc_altivec.h"
  44 #include "nb_kernel131_ppc_altivec.h"
  45 #include "nb_kernel132_ppc_altivec.h"
  46 #include "nb_kernel133_ppc_altivec.h"
  47 #include "nb_kernel134_ppc_altivec.h"
  48 #include "nb_kernel200_ppc_altivec.h"
  49 #include "nb_kernel201_ppc_altivec.h"
  50 #include "nb_kernel202_ppc_altivec.h"
  51 #include "nb_kernel203_ppc_altivec.h"
  52 #include "nb_kernel204_ppc_altivec.h"
  53 #include "nb_kernel210_ppc_altivec.h"
  54 #include "nb_kernel211_ppc_altivec.h"
  55 #include "nb_kernel212_ppc_altivec.h"
  56 #include "nb_kernel213_ppc_altivec.h"
  57 #include "nb_kernel214_ppc_altivec.h"
  58 #include "nb_kernel230_ppc_altivec.h"
  59 #include "nb_kernel231_ppc_altivec.h"
  60 #include "nb_kernel232_ppc_altivec.h"
  61 #include "nb_kernel233_ppc_altivec.h"
  62 #include "nb_kernel234_ppc_altivec.h"
  63 #include "nb_kernel300_ppc_altivec.h"
  64 #include "nb_kernel301_ppc_altivec.h"
  65 #include "nb_kernel302_ppc_altivec.h"
  66 #include "nb_kernel303_ppc_altivec.h"
  67 #include "nb_kernel304_ppc_altivec.h"
  68 #include "nb_kernel310_ppc_altivec.h"
  69 #include "nb_kernel311_ppc_altivec.h"
  70 #include "nb_kernel312_ppc_altivec.h"
  71 #include "nb_kernel313_ppc_altivec.h"
  72 #include "nb_kernel314_ppc_altivec.h"
  73 #include "nb_kernel330_ppc_altivec.h"
  74 #include "nb_kernel331_ppc_altivec.h"
  75 #include "nb_kernel332_ppc_altivec.h"
  76 #include "nb_kernel333_ppc_altivec.h"
  77 #include "nb_kernel334_ppc_altivec.h"
  78 #include "nb_kernel400_ppc_altivec.h"
  79 #include "nb_kernel410_ppc_altivec.h"
  80 #include "nb_kernel430_ppc_altivec.h"
  81
  82 /* Necessary headers for POSIX-style long jumps. */
  83 #include <signal.h>
  84 #include <setjmp.h>
  85
  86
  87 static nb_kernel_t *
  88 kernellist_ppc_altivec[eNR_NBKERNEL_NR] =
  89 {
  90     nb_kernel010_ppc_altivec,
  91     NULL,
  92     nb_kernel030_ppc_altivec,
  93     nb_kernel100_ppc_altivec,
  94     nb_kernel101_ppc_altivec,
  95     nb_kernel102_ppc_altivec,
  96     nb_kernel103_ppc_altivec,
  97     nb_kernel104_ppc_altivec,
  98     nb_kernel110_ppc_altivec,
  99     nb_kernel111_ppc_altivec,
 100     nb_kernel112_ppc_altivec,
 101     nb_kernel113_ppc_altivec,
 102     nb_kernel114_ppc_altivec,
 103     NULL,
 104     NULL,
 105     NULL,
 106     NULL,
 107     NULL,
 108     nb_kernel130_ppc_altivec,
 109     nb_kernel131_ppc_altivec,
 110     nb_kernel132_ppc_altivec,
 111     nb_kernel133_ppc_altivec,
 112     nb_kernel134_ppc_altivec,
 113     nb_kernel200_ppc_altivec,
 114     nb_kernel201_ppc_altivec,
 115     nb_kernel202_ppc_altivec,
 116     nb_kernel203_ppc_altivec,
 117     nb_kernel204_ppc_altivec,
 118     nb_kernel210_ppc_altivec,
 119     nb_kernel211_ppc_altivec,
 120     nb_kernel212_ppc_altivec,
 121     nb_kernel213_ppc_altivec,
 122     nb_kernel214_ppc_altivec,
 123     NULL,
 124     NULL,
 125     NULL,
 126     NULL,
 127     NULL,
 128     nb_kernel230_ppc_altivec,
 129     nb_kernel231_ppc_altivec,
 130     nb_kernel232_ppc_altivec,
 131     nb_kernel233_ppc_altivec,
 132     nb_kernel234_ppc_altivec,
 133     nb_kernel300_ppc_altivec,
 134     nb_kernel301_ppc_altivec,
 135     nb_kernel302_ppc_altivec,
 136     nb_kernel303_ppc_altivec,
 137     nb_kernel304_ppc_altivec,
 138     nb_kernel310_ppc_altivec,
 139     nb_kernel311_ppc_altivec,
 140     nb_kernel312_ppc_altivec,
 141     nb_kernel313_ppc_altivec,
 142     nb_kernel314_ppc_altivec,
 143     NULL,
 144     NULL,
 145     NULL,
 146     NULL,
 147     NULL,
 148     nb_kernel330_ppc_altivec,
 149     nb_kernel331_ppc_altivec,
 150     nb_kernel332_ppc_altivec,
 151     nb_kernel333_ppc_altivec,
 152     nb_kernel334_ppc_altivec,
 153     nb_kernel400_ppc_altivec,
 154     nb_kernel410_ppc_altivec,
 155     nb_kernel430_ppc_altivec
 156 };
 157
 158
 159 static jmp_buf
 160 nb_kernel_ppc_altivec_testprog;
 161
 162
 163 /*! Result of Altivec test */
 164 static gmx_bool
 165 nb_kernel_ppc_altivec_present;
 166
 167
 168 static void
 169 nb_kernel_ppc_altivec_sigill_handler(int n)
 170 {
 171     nb_kernel_ppc_altivec_present=FALSE;
 172     longjmp(nb_kernel_ppc_altivec_testprog,n);
 173 }
 174
 175
 176
 177 /* Return 0 if Altivec support is present, or
 178  *  non-zero upon error.
 179  */
 180 int
 181 nb_kernel_ppc_altivec_test(FILE * log)
 182 {
 183
 184     if(log)
 185         fprintf(log,"Testing Altivec/VMX support...");
 186
 187     nb_kernel_ppc_altivec_present = TRUE;
 188     signal(SIGILL,nb_kernel_ppc_altivec_sigill_handler);
 189
 190     /* return to this point after executing the signal handler
 191      * if we catch a SIGILL
 192      */
 193     setjmp(nb_kernel_ppc_altivec_testprog);
 194
 195     if(nb_kernel_ppc_altivec_present)
 196         nb_kernel_ppc_altivec_issue_instructions();
 197
 198     /* If altivec worked, then success is still 1.
 199         * If we got SIGILL, it was set to 0 in sigill_handler().
 200         */
 201
 202     if(log)
 203         fprintf(log," %spresent.\n",
 204                 nb_kernel_ppc_altivec_present ? "":"not ");
 205
 206
 207     return ((nb_kernel_ppc_altivec_present) ? 0 : -1);
 208 }
 209
 210
 211
 212
 213
 214 void
 215 nb_kernel_setup_ppc_altivec(FILE *log,nb_kernel_t **list)
 216 {
 217     int i;
 218     nb_kernel_t *p;
 219
 220     if(nb_kernel_ppc_altivec_test(log) != 0)
 221         return;
 222
 223         if(log)
 224                 fprintf(log,"Configuring PPC/Altivec nonbonded kernels...\n");
 225
 226     for(i=0;i<eNR_NBKERNEL_NR;i++)
 227     {
 228         p = kernellist_ppc_altivec[i];
 229         if(p!=NULL)
 230             list[i] = p;
 231     }
 232 }
 233
 234
 235
 236
 237