Attachment 159770 Details for Bug 249250 – Preprocessed source

Preprocessed source

ccFHwxfj.out (text/plain), 108.46 KB, created by Arun Mahajan on 2007-07-23 10:09:29 UTC

(hide)

Description:

Filename:

MIME Type:

Creator: Arun Mahajan

Created: 2007-07-23 10:09:29 UTC

Size: 108.46 KB

patch

obsolete

>// /usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/cc1 -lang-c -I../util -I../matrices -I/software/mpbi/x86_64/install/libmatheval-1.1.4/include -I/software/mpbi/x86_64/install/fftw-2.1.5/include -I/software/mpbi/x86_64/install/libctl-3.0.2/include -D__GNUC__=3 -D__GNUC_MINOR__=2 -D__GNUC_PATCHLEVEL__=3 -D__GXX_ABI_VERSION=102 -D__ELF__ -Dunix -D__gnu_linux__ -Dlinux -D__ELF__ -D__unix__ -D__gnu_linux__ -D__linux__ -D__unix -D__linux -Asystem(posix) -D__OPTIMIZE__ -D__STDC_HOSTED__=1 -Acpu=x86_64 -Amachine=x86_64 -D__x86_64 -D__x86_64__ -D__SIZE_TYPE__=unsigned long int -D__PTRDIFF_TYPE__=long int -D__tune_athlon__ -D__tune_athlon_sse__ -D__LONG_MAX__=9223372036854775807L -D__LP64__ -D_LP64 -DHAVE_CONFIG_H maxwell_eps.c -quiet -dumpbase maxwell_eps.c -g -O2 -Wall -W -Wbad-function-cast -Wcast-qual -Wpointer-arith -Wcast-align -pedantic -o -
># 1 "maxwell_eps.c"
># 1 "<built-in>"
># 1 "<command line>"
># 1 "maxwell_eps.c"
># 18 "maxwell_eps.c"
># 1 "/usr/include/stdlib.h" 1 3
># 25 "/usr/include/stdlib.h" 3
># 1 "/usr/include/features.h" 1 3
># 296 "/usr/include/features.h" 3
># 1 "/usr/include/sys/cdefs.h" 1 3
># 297 "/usr/include/features.h" 2 3
># 319 "/usr/include/features.h" 3
># 1 "/usr/include/gnu/stubs.h" 1 3
># 320 "/usr/include/features.h" 2 3
># 26 "/usr/include/stdlib.h" 2 3
>
>
>
>
>
>
>
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 213 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 3
>typedef unsigned long int size_t;
># 325 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 3
>typedef int wchar_t;
># 34 "/usr/include/stdlib.h" 2 3
>
>
># 93 "/usr/include/stdlib.h" 3
>
>
>typedef struct
>  {
>    int quot;
>    int rem;
>  } div_t;
>
>
>
>typedef struct
>  {
>    long int quot;
>    long int rem;
>  } ldiv_t;
>
>
>
># 137 "/usr/include/stdlib.h" 3
>extern size_t __ctype_get_mb_cur_max (void) ;
>
>
>
>
>extern double atof (__const char *__nptr) __attribute__ ((__pure__));
>
>extern int atoi (__const char *__nptr) __attribute__ ((__pure__));
>
>extern long int atol (__const char *__nptr) __attribute__ ((__pure__));
>
>
>
>
>
>__extension__ extern long long int atoll (__const char *__nptr)
>     __attribute__ ((__pure__));
>
>
>
>
>
>extern double strtod (__const char *__restrict __nptr,
>                      char **__restrict __endptr) ;
>
># 174 "/usr/include/stdlib.h" 3
>
>
>extern long int strtol (__const char *__restrict __nptr,
>                        char **__restrict __endptr, int __base) ;
>
>extern unsigned long int strtoul (__const char *__restrict __nptr,
>                                  char **__restrict __endptr, int __base)
>     ;
>
>
>
>
>__extension__
>extern long long int strtoq (__const char *__restrict __nptr,
>                             char **__restrict __endptr, int __base) ;
>
>__extension__
>extern unsigned long long int strtouq (__const char *__restrict __nptr,
>                                       char **__restrict __endptr, int __base)
>     ;
>
>
>
>
>
>__extension__
>extern long long int strtoll (__const char *__restrict __nptr,
>                              char **__restrict __endptr, int __base) ;
>
>__extension__
>extern unsigned long long int strtoull (__const char *__restrict __nptr,
>                                        char **__restrict __endptr, int __base)
>     ;
>
># 264 "/usr/include/stdlib.h" 3
>extern double __strtod_internal (__const char *__restrict __nptr,
>                                 char **__restrict __endptr, int __group)
>     ;
>extern float __strtof_internal (__const char *__restrict __nptr,
>                                char **__restrict __endptr, int __group)
>     ;
>extern long double __strtold_internal (__const char *__restrict __nptr,
>                                       char **__restrict __endptr,
>                                       int __group) ;
>
>extern long int __strtol_internal (__const char *__restrict __nptr,
>                                   char **__restrict __endptr,
>                                   int __base, int __group) ;
>
>
>
>extern unsigned long int __strtoul_internal (__const char *__restrict __nptr,
>                                             char **__restrict __endptr,
>                                             int __base, int __group) ;
>
>
>
>
>__extension__
>extern long long int __strtoll_internal (__const char *__restrict __nptr,
>                                         char **__restrict __endptr,
>                                         int __base, int __group) ;
>
>
>
>__extension__
>extern unsigned long long int __strtoull_internal (__const char *
>                                                   __restrict __nptr,
>                                                   char **__restrict __endptr,
>                                                   int __base, int __group)
>     ;
>
>
>
>
>
>
>
>
>extern __inline double
>strtod (__const char *__restrict __nptr, char **__restrict __endptr)
>{
>  return __strtod_internal (__nptr, __endptr, 0);
>}
>extern __inline long int
>strtol (__const char *__restrict __nptr, char **__restrict __endptr,
>        int __base)
>{
>  return __strtol_internal (__nptr, __endptr, __base, 0);
>}
>extern __inline unsigned long int
>strtoul (__const char *__restrict __nptr, char **__restrict __endptr,
>         int __base)
>{
>  return __strtoul_internal (__nptr, __endptr, __base, 0);
>}
>
># 343 "/usr/include/stdlib.h" 3
>__extension__ extern __inline long long int
>strtoq (__const char *__restrict __nptr, char **__restrict __endptr,
>        int __base)
>{
>  return __strtoll_internal (__nptr, __endptr, __base, 0);
>}
>__extension__ extern __inline unsigned long long int
>strtouq (__const char *__restrict __nptr, char **__restrict __endptr,
>         int __base)
>{
>  return __strtoull_internal (__nptr, __endptr, __base, 0);
>}
>
>
>
>
>__extension__ extern __inline long long int
>strtoll (__const char *__restrict __nptr, char **__restrict __endptr,
>         int __base)
>{
>  return __strtoll_internal (__nptr, __endptr, __base, 0);
>}
>__extension__ extern __inline unsigned long long int
>strtoull (__const char * __restrict __nptr, char **__restrict __endptr,
>          int __base)
>{
>  return __strtoull_internal (__nptr, __endptr, __base, 0);
>}
>
>
>
>
>extern __inline double
>atof (__const char *__nptr)
>{
>  return strtod (__nptr, (char **) ((void *)0));
>}
>extern __inline int
>atoi (__const char *__nptr)
>{
>  return (int) strtol (__nptr, (char **) ((void *)0), 10);
>}
>extern __inline long int
>atol (__const char *__nptr)
>{
>  return strtol (__nptr, (char **) ((void *)0), 10);
>}
>
>
>
>
>__extension__ extern __inline long long int
>atoll (__const char *__nptr)
>{
>  return strtoll (__nptr, (char **) ((void *)0), 10);
>}
>
># 408 "/usr/include/stdlib.h" 3
>extern char *l64a (long int __n) ;
>
>
>extern long int a64l (__const char *__s) __attribute__ ((__pure__));
>
>
>
>
># 1 "/usr/include/sys/types.h" 1 3
># 29 "/usr/include/sys/types.h" 3
>
>
># 1 "/usr/include/bits/types.h" 1 3
># 28 "/usr/include/bits/types.h" 3
># 1 "/usr/include/bits/wordsize.h" 1 3
># 29 "/usr/include/bits/types.h" 2 3
>
>
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 32 "/usr/include/bits/types.h" 2 3
>
>
>typedef unsigned char __u_char;
>typedef unsigned short int __u_short;
>typedef unsigned int __u_int;
>typedef unsigned long int __u_long;
>
>
>typedef signed char __int8_t;
>typedef unsigned char __uint8_t;
>typedef signed short int __int16_t;
>typedef unsigned short int __uint16_t;
>typedef signed int __int32_t;
>typedef unsigned int __uint32_t;
>
>typedef signed long int __int64_t;
>typedef unsigned long int __uint64_t;
>
>
>
>
>
>
>
>typedef long int __quad_t;
>typedef unsigned long int __u_quad_t;
># 129 "/usr/include/bits/types.h" 3
># 1 "/usr/include/bits/typesizes.h" 1 3
># 130 "/usr/include/bits/types.h" 2 3
>
>
>
>
>
>
>__extension__ typedef unsigned long int __dev_t;
>__extension__ typedef unsigned int __uid_t;
>__extension__ typedef unsigned int __gid_t;
>__extension__ typedef unsigned long int __ino_t;
>__extension__ typedef unsigned long int __ino64_t;
>__extension__ typedef unsigned int __mode_t;
>__extension__ typedef unsigned long int __nlink_t;
>__extension__ typedef long int __off_t;
>__extension__ typedef long int __off64_t;
>__extension__ typedef int __pid_t;
>__extension__ typedef struct { int __val[2]; } __fsid_t;
>__extension__ typedef long int __clock_t;
>__extension__ typedef unsigned long int __rlim_t;
>__extension__ typedef unsigned long int __rlim64_t;
>__extension__ typedef unsigned int __id_t;
>__extension__ typedef long int __time_t;
>__extension__ typedef unsigned int __useconds_t;
>__extension__ typedef long int __suseconds_t;
>
>__extension__ typedef int __daddr_t;
>__extension__ typedef long int __swblk_t;
>__extension__ typedef int __key_t;
>
>
>__extension__ typedef int __clockid_t;
>
>
>__extension__ typedef int __timer_t;
>
>
>__extension__ typedef long int __blksize_t;
>
>
>
>
>__extension__ typedef long int __blkcnt_t;
>__extension__ typedef long int __blkcnt64_t;
>
>
>__extension__ typedef unsigned long int __fsblkcnt_t;
>__extension__ typedef unsigned long int __fsblkcnt64_t;
>
>
>__extension__ typedef unsigned long int __fsfilcnt_t;
>__extension__ typedef unsigned long int __fsfilcnt64_t;
>
>__extension__ typedef long int __ssize_t;
>
>
>
>typedef __off64_t __loff_t;
>typedef __quad_t *__qaddr_t;
>typedef char *__caddr_t;
>
>
>__extension__ typedef long int __intptr_t;
>
>
>__extension__ typedef unsigned int __socklen_t;
># 32 "/usr/include/sys/types.h" 2 3
>
>
>
>typedef __u_char u_char;
>typedef __u_short u_short;
>typedef __u_int u_int;
>typedef __u_long u_long;
>typedef __quad_t quad_t;
>typedef __u_quad_t u_quad_t;
>typedef __fsid_t fsid_t;
>
>
>
>
>typedef __loff_t loff_t;
>
>
>
>typedef __ino_t ino_t;
># 62 "/usr/include/sys/types.h" 3
>typedef __dev_t dev_t;
>
>
>
>
>typedef __gid_t gid_t;
>
>
>
>
>typedef __mode_t mode_t;
>
>
>
>
>typedef __nlink_t nlink_t;
>
>
>
>
>typedef __uid_t uid_t;
>
>
>
>
>
>typedef __off_t off_t;
># 100 "/usr/include/sys/types.h" 3
>typedef __pid_t pid_t;
>
>
>
>
>typedef __id_t id_t;
>
>
>
>
>typedef __ssize_t ssize_t;
>
>
>
>
>
>typedef __daddr_t daddr_t;
>typedef __caddr_t caddr_t;
>
>
>
>
>
>typedef __key_t key_t;
># 133 "/usr/include/sys/types.h" 3
># 1 "/usr/include/time.h" 1 3
># 74 "/usr/include/time.h" 3
>
>
>typedef __time_t time_t;
>
>
>
># 92 "/usr/include/time.h" 3
>typedef __clockid_t clockid_t;
># 104 "/usr/include/time.h" 3
>typedef __timer_t timer_t;
># 134 "/usr/include/sys/types.h" 2 3
># 147 "/usr/include/sys/types.h" 3
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 148 "/usr/include/sys/types.h" 2 3
>
>
>
>typedef unsigned long int ulong;
>typedef unsigned short int ushort;
>typedef unsigned int uint;
># 191 "/usr/include/sys/types.h" 3
>typedef int int8_t __attribute__ ((__mode__ (__QI__)));
>typedef int int16_t __attribute__ ((__mode__ (__HI__)));
>typedef int int32_t __attribute__ ((__mode__ (__SI__)));
>typedef int int64_t __attribute__ ((__mode__ (__DI__)));
>
>
>typedef unsigned int u_int8_t __attribute__ ((__mode__ (__QI__)));
>typedef unsigned int u_int16_t __attribute__ ((__mode__ (__HI__)));
>typedef unsigned int u_int32_t __attribute__ ((__mode__ (__SI__)));
>typedef unsigned int u_int64_t __attribute__ ((__mode__ (__DI__)));
>
>typedef int register_t __attribute__ ((__mode__ (__word__)));
># 213 "/usr/include/sys/types.h" 3
># 1 "/usr/include/endian.h" 1 3
># 37 "/usr/include/endian.h" 3
># 1 "/usr/include/bits/endian.h" 1 3
># 38 "/usr/include/endian.h" 2 3
># 214 "/usr/include/sys/types.h" 2 3
>
>
># 1 "/usr/include/sys/select.h" 1 3
># 31 "/usr/include/sys/select.h" 3
># 1 "/usr/include/bits/select.h" 1 3
># 32 "/usr/include/sys/select.h" 2 3
>
>
># 1 "/usr/include/bits/sigset.h" 1 3
># 23 "/usr/include/bits/sigset.h" 3
>typedef int __sig_atomic_t;
>
>
>
>
>typedef struct
>  {
>    unsigned long int __val[(1024 / (8 * sizeof (unsigned long int)))];
>  } __sigset_t;
># 35 "/usr/include/sys/select.h" 2 3
>
>
>
>typedef __sigset_t sigset_t;
>
>
>
>
>
># 1 "/usr/include/time.h" 1 3
># 118 "/usr/include/time.h" 3
>struct timespec
>  {
>    __time_t tv_sec;
>    long int tv_nsec;
>  };
># 45 "/usr/include/sys/select.h" 2 3
>
># 1 "/usr/include/bits/time.h" 1 3
># 69 "/usr/include/bits/time.h" 3
>struct timeval
>  {
>    __time_t tv_sec;
>    __suseconds_t tv_usec;
>  };
># 47 "/usr/include/sys/select.h" 2 3
>
>
>typedef __suseconds_t suseconds_t;
>
>
>
>
>
>typedef long int __fd_mask;
># 67 "/usr/include/sys/select.h" 3
>typedef struct
>  {
>
>
>
>
>
>
>    __fd_mask __fds_bits[1024 / (8 * sizeof (__fd_mask))];
>
>
>  } fd_set;
>
>
>
>
>
>
>typedef __fd_mask fd_mask;
># 99 "/usr/include/sys/select.h" 3
>
># 109 "/usr/include/sys/select.h" 3
>extern int select (int __nfds, fd_set *__restrict __readfds,
>                   fd_set *__restrict __writefds,
>                   fd_set *__restrict __exceptfds,
>                   struct timeval *__restrict __timeout);
># 128 "/usr/include/sys/select.h" 3
>
># 217 "/usr/include/sys/types.h" 2 3
>
>
># 1 "/usr/include/sys/sysmacros.h" 1 3
># 29 "/usr/include/sys/sysmacros.h" 3
>__extension__
>extern __inline unsigned int gnu_dev_major (unsigned long long int __dev)
>     ;
>__extension__
>extern __inline unsigned int gnu_dev_minor (unsigned long long int __dev)
>     ;
>__extension__
>extern __inline unsigned long long int gnu_dev_makedev (unsigned int __major,
>                                                        unsigned int __minor)
>     ;
>
>
>__extension__ extern __inline unsigned int
>gnu_dev_major (unsigned long long int __dev)
>{
>  return ((__dev >> 8) & 0xfff) | ((unsigned int) (__dev >> 32) & ~0xfff);
>}
>
>__extension__ extern __inline unsigned int
>gnu_dev_minor (unsigned long long int __dev)
>{
>  return (__dev & 0xff) | ((unsigned int) (__dev >> 12) & ~0xff);
>}
>
>__extension__ extern __inline unsigned long long int
>gnu_dev_makedev (unsigned int __major, unsigned int __minor)
>{
>  return ((__minor & 0xff) | ((__major & 0xfff) << 8)
>          | (((unsigned long long int) (__minor & ~0xff)) << 12)
>          | (((unsigned long long int) (__major & ~0xfff)) << 32));
>}
># 220 "/usr/include/sys/types.h" 2 3
># 231 "/usr/include/sys/types.h" 3
>typedef __blkcnt_t blkcnt_t;
>
>
>
>typedef __fsblkcnt_t fsblkcnt_t;
>
>
>
>typedef __fsfilcnt_t fsfilcnt_t;
># 266 "/usr/include/sys/types.h" 3
># 1 "/usr/include/bits/pthreadtypes.h" 1 3
># 23 "/usr/include/bits/pthreadtypes.h" 3
># 1 "/usr/include/bits/sched.h" 1 3
># 83 "/usr/include/bits/sched.h" 3
>struct __sched_param
>  {
>    int __sched_priority;
>  };
># 24 "/usr/include/bits/pthreadtypes.h" 2 3
>
>
>struct _pthread_fastlock
>{
>  long int __status;
>  int __spinlock;
>
>};
>
>
>
>typedef struct _pthread_descr_struct *_pthread_descr;
>
>
>
>
>
>typedef struct __pthread_attr_s
>{
>  int __detachstate;
>  int __schedpolicy;
>  struct __sched_param __schedparam;
>  int __inheritsched;
>  int __scope;
>  size_t __guardsize;
>  int __stackaddr_set;
>  void *__stackaddr;
>  size_t __stacksize;
>} pthread_attr_t;
>
>
>
>
>
>__extension__ typedef long long __pthread_cond_align_t;
>
>
>
>
>typedef struct
>{
>  struct _pthread_fastlock __c_lock;
>  _pthread_descr __c_waiting;
>  char __padding[48 - sizeof (struct _pthread_fastlock)
>                 - sizeof (_pthread_descr) - sizeof (__pthread_cond_align_t)];
>  __pthread_cond_align_t __align;
>} pthread_cond_t;
>
>
>
>typedef struct
>{
>  int __dummy;
>} pthread_condattr_t;
>
>
>typedef unsigned int pthread_key_t;
>
>
>
>
>
>typedef struct
>{
>  int __m_reserved;
>  int __m_count;
>  _pthread_descr __m_owner;
>  int __m_kind;
>  struct _pthread_fastlock __m_lock;
>} pthread_mutex_t;
>
>
>
>typedef struct
>{
>  int __mutexkind;
>} pthread_mutexattr_t;
>
>
>
>typedef int pthread_once_t;
># 150 "/usr/include/bits/pthreadtypes.h" 3
>typedef unsigned long int pthread_t;
># 267 "/usr/include/sys/types.h" 2 3
>
>
>
># 417 "/usr/include/stdlib.h" 2 3
>
>
>
>
>
>
>extern long int random (void) ;
>
>
>extern void srandom (unsigned int __seed) ;
>
>
>
>
>
>extern char *initstate (unsigned int __seed, char *__statebuf,
>                        size_t __statelen) ;
>
>
>
>extern char *setstate (char *__statebuf) ;
>
>
>
>
>
>
>
>struct random_data
>  {
>    int32_t *fptr;
>    int32_t *rptr;
>    int32_t *state;
>    int rand_type;
>    int rand_deg;
>    int rand_sep;
>    int32_t *end_ptr;
>  };
>
>extern int random_r (struct random_data *__restrict __buf,
>                     int32_t *__restrict __result) ;
>
>extern int srandom_r (unsigned int __seed, struct random_data *__buf) ;
>
>extern int initstate_r (unsigned int __seed, char *__restrict __statebuf,
>                        size_t __statelen,
>                        struct random_data *__restrict __buf) ;
>
>extern int setstate_r (char *__restrict __statebuf,
>                       struct random_data *__restrict __buf) ;
>
>
>
>
>
>
>extern int rand (void) ;
>
>extern void srand (unsigned int __seed) ;
>
>
>
>
>extern int rand_r (unsigned int *__seed) ;
>
>
>
>
>
>
>
>extern double drand48 (void) ;
>extern double erand48 (unsigned short int __xsubi[3]) ;
>
>
>extern long int lrand48 (void) ;
>extern long int nrand48 (unsigned short int __xsubi[3]) ;
>
>
>extern long int mrand48 (void) ;
>extern long int jrand48 (unsigned short int __xsubi[3]) ;
>
>
>extern void srand48 (long int __seedval) ;
>extern unsigned short int *seed48 (unsigned short int __seed16v[3]) ;
>extern void lcong48 (unsigned short int __param[7]) ;
>
>
>
>
>
>struct drand48_data
>  {
>    unsigned short int __x[3];
>    unsigned short int __old_x[3];
>    unsigned short int __c;
>    unsigned short int __init;
>    unsigned long long int __a;
>  };
>
>
>extern int drand48_r (struct drand48_data *__restrict __buffer,
>                      double *__restrict __result) ;
>extern int erand48_r (unsigned short int __xsubi[3],
>                      struct drand48_data *__restrict __buffer,
>                      double *__restrict __result) ;
>
>
>extern int lrand48_r (struct drand48_data *__restrict __buffer,
>                      long int *__restrict __result) ;
>extern int nrand48_r (unsigned short int __xsubi[3],
>                      struct drand48_data *__restrict __buffer,
>                      long int *__restrict __result) ;
>
>
>extern int mrand48_r (struct drand48_data *__restrict __buffer,
>                      long int *__restrict __result) ;
>extern int jrand48_r (unsigned short int __xsubi[3],
>                      struct drand48_data *__restrict __buffer,
>                      long int *__restrict __result) ;
>
>
>extern int srand48_r (long int __seedval, struct drand48_data *__buffer)
>     ;
>
>extern int seed48_r (unsigned short int __seed16v[3],
>                     struct drand48_data *__buffer) ;
>
>extern int lcong48_r (unsigned short int __param[7],
>                      struct drand48_data *__buffer) ;
>
>
>
>
>
>
>
>
>
>extern void *malloc (size_t __size) __attribute__ ((__malloc__));
>
>extern void *calloc (size_t __nmemb, size_t __size)
>     __attribute__ ((__malloc__));
>
>
>
>
>
>
>
>extern void *realloc (void *__ptr, size_t __size) __attribute__ ((__malloc__));
>
>extern void free (void *__ptr) ;
>
>
>
>
>extern void cfree (void *__ptr) ;
>
>
>
># 1 "/usr/include/alloca.h" 1 3
># 25 "/usr/include/alloca.h" 3
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 26 "/usr/include/alloca.h" 2 3
>
>
>
>
>
>
>
>extern void *alloca (size_t __size) ;
>
>
>
>
>
>
># 579 "/usr/include/stdlib.h" 2 3
>
>
>
>
>extern void *valloc (size_t __size) __attribute__ ((__malloc__));
># 592 "/usr/include/stdlib.h" 3
>
>
>extern void abort (void) __attribute__ ((__noreturn__));
>
>
>
>extern int atexit (void (*__func) (void)) ;
>
>
>
>
>
>extern int on_exit (void (*__func) (int __status, void *__arg), void *__arg)
>     ;
>
>
>
>
>
>
>extern void exit (int __status) __attribute__ ((__noreturn__));
>
># 624 "/usr/include/stdlib.h" 3
>
>
>extern char *getenv (__const char *__name) ;
>
>
>
>
>extern char *__secure_getenv (__const char *__name) ;
>
>
>
>
>
>extern int putenv (char *__string) ;
>
>
>
>
>
>extern int setenv (__const char *__name, __const char *__value, int __replace)
>     ;
>
>
>extern int unsetenv (__const char *__name) ;
>
>
>
>
>
>
>extern int clearenv (void) ;
># 663 "/usr/include/stdlib.h" 3
>extern char *mktemp (char *__template) ;
># 674 "/usr/include/stdlib.h" 3
>extern int mkstemp (char *__template);
># 693 "/usr/include/stdlib.h" 3
>extern char *mkdtemp (char *__template) ;
>
>
>
>
>
>
>
>
>extern int system (__const char *__command);
>
># 720 "/usr/include/stdlib.h" 3
>extern char *realpath (__const char *__restrict __name,
>                       char *__restrict __resolved) ;
>
>
>
>
>
>
>typedef int (*__compar_fn_t) (__const void *, __const void *);
>
>
>
>
>
>
>
>
>
>extern void *bsearch (__const void *__key, __const void *__base,
>                      size_t __nmemb, size_t __size, __compar_fn_t __compar);
>
>
>
>extern void qsort (void *__base, size_t __nmemb, size_t __size,
>                   __compar_fn_t __compar);
>
>
>
>extern int abs (int __x) __attribute__ ((__const__));
>extern long int labs (long int __x) __attribute__ ((__const__));
>
>
>
>
>
>
>
>
>
>
>
>
>extern div_t div (int __numer, int __denom)
>     __attribute__ ((__const__));
>extern ldiv_t ldiv (long int __numer, long int __denom)
>     __attribute__ ((__const__));
>
># 784 "/usr/include/stdlib.h" 3
>extern char *ecvt (double __value, int __ndigit, int *__restrict __decpt,
>                   int *__restrict __sign) ;
>
>
>
>
>extern char *fcvt (double __value, int __ndigit, int *__restrict __decpt,
>                   int *__restrict __sign) ;
>
>
>
>
>extern char *gcvt (double __value, int __ndigit, char *__buf) ;
>
>
>
>
>extern char *qecvt (long double __value, int __ndigit,
>                    int *__restrict __decpt, int *__restrict __sign) ;
>extern char *qfcvt (long double __value, int __ndigit,
>                    int *__restrict __decpt, int *__restrict __sign) ;
>extern char *qgcvt (long double __value, int __ndigit, char *__buf) ;
>
>
>
>
>extern int ecvt_r (double __value, int __ndigit, int *__restrict __decpt,
>                   int *__restrict __sign, char *__restrict __buf,
>                   size_t __len) ;
>extern int fcvt_r (double __value, int __ndigit, int *__restrict __decpt,
>                   int *__restrict __sign, char *__restrict __buf,
>                   size_t __len) ;
>
>extern int qecvt_r (long double __value, int __ndigit,
>                    int *__restrict __decpt, int *__restrict __sign,
>                    char *__restrict __buf, size_t __len) ;
>extern int qfcvt_r (long double __value, int __ndigit,
>                    int *__restrict __decpt, int *__restrict __sign,
>                    char *__restrict __buf, size_t __len) ;
>
>
>
>
>
>
>
>extern int mblen (__const char *__s, size_t __n) ;
>
>
>extern int mbtowc (wchar_t *__restrict __pwc,
>                   __const char *__restrict __s, size_t __n) ;
>
>
>extern int wctomb (char *__s, wchar_t __wchar) ;
>
>
>
>extern size_t mbstowcs (wchar_t *__restrict __pwcs,
>                        __const char *__restrict __s, size_t __n) ;
>
>extern size_t wcstombs (char *__restrict __s,
>                        __const wchar_t *__restrict __pwcs, size_t __n)
>     ;
>
>
>
>
>
>
>
>
>extern int rpmatch (__const char *__response) ;
># 916 "/usr/include/stdlib.h" 3
>extern int getloadavg (double __loadavg[], int __nelem) ;
>
>
>
>
>
>
># 19 "maxwell_eps.c" 2
># 1 "/usr/include/stdio.h" 1 3
># 30 "/usr/include/stdio.h" 3
>
>
>
>
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 35 "/usr/include/stdio.h" 2 3
># 44 "/usr/include/stdio.h" 3
>
>
>typedef struct _IO_FILE FILE;
>
>
>
>
>
># 62 "/usr/include/stdio.h" 3
>typedef struct _IO_FILE __FILE;
># 72 "/usr/include/stdio.h" 3
># 1 "/usr/include/libio.h" 1 3
># 32 "/usr/include/libio.h" 3
># 1 "/usr/include/_G_config.h" 1 3
># 14 "/usr/include/_G_config.h" 3
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 354 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 3
>typedef unsigned int wint_t;
># 15 "/usr/include/_G_config.h" 2 3
># 24 "/usr/include/_G_config.h" 3
># 1 "/usr/include/wchar.h" 1 3
># 48 "/usr/include/wchar.h" 3
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 49 "/usr/include/wchar.h" 2 3
>
># 1 "/usr/include/bits/wchar.h" 1 3
># 51 "/usr/include/wchar.h" 2 3
># 76 "/usr/include/wchar.h" 3
>typedef struct
>{
>  int __count;
>  union
>  {
>    wint_t __wch;
>    char __wchb[4];
>  } __value;
>} __mbstate_t;
># 25 "/usr/include/_G_config.h" 2 3
>
>typedef struct
>{
>  __off_t __pos;
>  __mbstate_t __state;
>} _G_fpos_t;
>typedef struct
>{
>  __off64_t __pos;
>  __mbstate_t __state;
>} _G_fpos64_t;
># 44 "/usr/include/_G_config.h" 3
># 1 "/usr/include/gconv.h" 1 3
># 28 "/usr/include/gconv.h" 3
># 1 "/usr/include/wchar.h" 1 3
># 48 "/usr/include/wchar.h" 3
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 49 "/usr/include/wchar.h" 2 3
># 29 "/usr/include/gconv.h" 2 3
>
>
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 32 "/usr/include/gconv.h" 2 3
>
>
>
>
>
>enum
>{
>  __GCONV_OK = 0,
>  __GCONV_NOCONV,
>  __GCONV_NODB,
>  __GCONV_NOMEM,
>
>  __GCONV_EMPTY_INPUT,
>  __GCONV_FULL_OUTPUT,
>  __GCONV_ILLEGAL_INPUT,
>  __GCONV_INCOMPLETE_INPUT,
>
>  __GCONV_ILLEGAL_DESCRIPTOR,
>  __GCONV_INTERNAL_ERROR
>};
>
>
>
>enum
>{
>  __GCONV_IS_LAST = 0x0001,
>  __GCONV_IGNORE_ERRORS = 0x0002
>};
>
>
>
>struct __gconv_step;
>struct __gconv_step_data;
>struct __gconv_loaded_object;
>struct __gconv_trans_data;
>
>
>
>typedef int (*__gconv_fct) (struct __gconv_step *, struct __gconv_step_data *,
>                            __const unsigned char **, __const unsigned char *,
>                            unsigned char **, size_t *, int, int);
>
>
>typedef wint_t (*__gconv_btowc_fct) (struct __gconv_step *, unsigned char);
>
>
>typedef int (*__gconv_init_fct) (struct __gconv_step *);
>typedef void (*__gconv_end_fct) (struct __gconv_step *);
>
>
>
>typedef int (*__gconv_trans_fct) (struct __gconv_step *,
>                                  struct __gconv_step_data *, void *,
>                                  __const unsigned char *,
>                                  __const unsigned char **,
>                                  __const unsigned char *, unsigned char **,
>                                  size_t *);
>
>
>typedef int (*__gconv_trans_context_fct) (void *, __const unsigned char *,
>                                          __const unsigned char *,
>                                          unsigned char *, unsigned char *);
>
>
>typedef int (*__gconv_trans_query_fct) (__const char *, __const char ***,
>                                        size_t *);
>
>
>typedef int (*__gconv_trans_init_fct) (void **, const char *);
>typedef void (*__gconv_trans_end_fct) (void *);
>
>struct __gconv_trans_data
>{
>
>  __gconv_trans_fct __trans_fct;
>  __gconv_trans_context_fct __trans_context_fct;
>  __gconv_trans_end_fct __trans_end_fct;
>  void *__data;
>  struct __gconv_trans_data *__next;
>};
>
>
>
>struct __gconv_step
>{
>  struct __gconv_loaded_object *__shlib_handle;
>  __const char *__modname;
>
>  int __counter;
>
>  char *__from_name;
>  char *__to_name;
>
>  __gconv_fct __fct;
>  __gconv_btowc_fct __btowc_fct;
>  __gconv_init_fct __init_fct;
>  __gconv_end_fct __end_fct;
>
>
>
>  int __min_needed_from;
>  int __max_needed_from;
>  int __min_needed_to;
>  int __max_needed_to;
>
>
>  int __stateful;
>
>  void *__data;
>};
>
>
>
>struct __gconv_step_data
>{
>  unsigned char *__outbuf;
>  unsigned char *__outbufend;
>
>
>
>  int __flags;
>
>
>
>  int __invocation_counter;
>
>
>
>  int __internal_use;
>
>  __mbstate_t *__statep;
>  __mbstate_t __state;
>
>
>
>  struct __gconv_trans_data *__trans;
>};
>
>
>
>typedef struct __gconv_info
>{
>  size_t __nsteps;
>  struct __gconv_step *__steps;
>  __extension__ struct __gconv_step_data __data [];
>} *__gconv_t;
># 45 "/usr/include/_G_config.h" 2 3
>typedef union
>{
>  struct __gconv_info __cd;
>  struct
>  {
>    struct __gconv_info __cd;
>    struct __gconv_step_data __data;
>  } __combined;
>} _G_iconv_t;
>
>typedef int _G_int16_t __attribute__ ((__mode__ (__HI__)));
>typedef int _G_int32_t __attribute__ ((__mode__ (__SI__)));
>typedef unsigned int _G_uint16_t __attribute__ ((__mode__ (__HI__)));
>typedef unsigned int _G_uint32_t __attribute__ ((__mode__ (__SI__)));
># 33 "/usr/include/libio.h" 2 3
># 53 "/usr/include/libio.h" 3
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stdarg.h" 1 3
># 43 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stdarg.h" 3
>typedef __builtin_va_list __gnuc_va_list;
># 54 "/usr/include/libio.h" 2 3
># 163 "/usr/include/libio.h" 3
>struct _IO_jump_t; struct _IO_FILE;
># 173 "/usr/include/libio.h" 3
>typedef void _IO_lock_t;
>
>
>
>
>
>struct _IO_marker {
>  struct _IO_marker *_next;
>  struct _IO_FILE *_sbuf;
>
>
>
>  int _pos;
># 196 "/usr/include/libio.h" 3
>};
>
>
>enum __codecvt_result
>{
>  __codecvt_ok,
>  __codecvt_partial,
>  __codecvt_error,
>  __codecvt_noconv
>};
># 264 "/usr/include/libio.h" 3
>struct _IO_FILE {
>  int _flags;
>
>
>
>
>  char* _IO_read_ptr;
>  char* _IO_read_end;
>  char* _IO_read_base;
>  char* _IO_write_base;
>  char* _IO_write_ptr;
>  char* _IO_write_end;
>  char* _IO_buf_base;
>  char* _IO_buf_end;
>
>  char *_IO_save_base;
>  char *_IO_backup_base;
>  char *_IO_save_end;
>
>  struct _IO_marker *_markers;
>
>  struct _IO_FILE *_chain;
>
>  int _fileno;
>
>
>
>  int _flags2;
>
>  __off_t _old_offset;
>
>
>
>  unsigned short _cur_column;
>  signed char _vtable_offset;
>  char _shortbuf[1];
>
>
>
>  _IO_lock_t *_lock;
># 312 "/usr/include/libio.h" 3
>  __off64_t _offset;
>
>
>
>
>
>  void *__pad1;
>  void *__pad2;
>
>  int _mode;
>
>  char _unused2[15 * sizeof (int) - 2 * sizeof (void *)];
>
>};
>
>
>typedef struct _IO_FILE _IO_FILE;
>
>
>struct _IO_FILE_plus;
>
>extern struct _IO_FILE_plus _IO_2_1_stdin_;
>extern struct _IO_FILE_plus _IO_2_1_stdout_;
>extern struct _IO_FILE_plus _IO_2_1_stderr_;
># 351 "/usr/include/libio.h" 3
>typedef __ssize_t __io_read_fn (void *__cookie, char *__buf, size_t __nbytes);
>
>
>
>
>
>
>
>typedef __ssize_t __io_write_fn (void *__cookie, __const char *__buf,
>                                 size_t __n);
>
>
>
>
>
>
>
>typedef int __io_seek_fn (void *__cookie, __off64_t *__pos, int __w);
>
>
>typedef int __io_close_fn (void *__cookie);
># 403 "/usr/include/libio.h" 3
>extern int __underflow (_IO_FILE *) ;
>extern int __uflow (_IO_FILE *) ;
>extern int __overflow (_IO_FILE *, int) ;
>extern wint_t __wunderflow (_IO_FILE *) ;
>extern wint_t __wuflow (_IO_FILE *) ;
>extern wint_t __woverflow (_IO_FILE *, wint_t) ;
># 433 "/usr/include/libio.h" 3
>extern int _IO_getc (_IO_FILE *__fp) ;
>extern int _IO_putc (int __c, _IO_FILE *__fp) ;
>extern int _IO_feof (_IO_FILE *__fp) ;
>extern int _IO_ferror (_IO_FILE *__fp) ;
>
>extern int _IO_peekc_locked (_IO_FILE *__fp) ;
>
>
>
>
>
>extern void _IO_flockfile (_IO_FILE *) ;
>extern void _IO_funlockfile (_IO_FILE *) ;
>extern int _IO_ftrylockfile (_IO_FILE *) ;
># 463 "/usr/include/libio.h" 3
>extern int _IO_vfscanf (_IO_FILE * __restrict, const char * __restrict,
>                        __gnuc_va_list, int *__restrict) ;
>extern int _IO_vfprintf (_IO_FILE *__restrict, const char *__restrict,
>                         __gnuc_va_list) ;
>extern __ssize_t _IO_padn (_IO_FILE *, int, __ssize_t) ;
>extern size_t _IO_sgetn (_IO_FILE *, void *, size_t) ;
>
>extern __off64_t _IO_seekoff (_IO_FILE *, __off64_t, int, int) ;
>extern __off64_t _IO_seekpos (_IO_FILE *, __off64_t, int) ;
>
>extern void _IO_free_backup_area (_IO_FILE *) ;
># 73 "/usr/include/stdio.h" 2 3
># 86 "/usr/include/stdio.h" 3
>
>
>typedef _G_fpos_t fpos_t;
>
>
>
>
># 138 "/usr/include/stdio.h" 3
># 1 "/usr/include/bits/stdio_lim.h" 1 3
># 139 "/usr/include/stdio.h" 2 3
>
>
>
>extern struct _IO_FILE *stdin;
>extern struct _IO_FILE *stdout;
>extern struct _IO_FILE *stderr;
>
>
>
>
>
>
>
>
>
>extern int remove (__const char *__filename) ;
>
>extern int rename (__const char *__old, __const char *__new) ;
>
>
>
>
>
>
>
>
>
>extern FILE *tmpfile (void);
># 180 "/usr/include/stdio.h" 3
>extern char *tmpnam (char *__s) ;
>
>
>
>
>
>extern char *tmpnam_r (char *__s) ;
># 198 "/usr/include/stdio.h" 3
>extern char *tempnam (__const char *__dir, __const char *__pfx)
>     __attribute__ ((__malloc__));
>
>
>
>
>
>
>
>
>extern int fclose (FILE *__stream);
>
>
>
>
>extern int fflush (FILE *__stream);
>
># 223 "/usr/include/stdio.h" 3
>extern int fflush_unlocked (FILE *__stream);
># 237 "/usr/include/stdio.h" 3
>
>
>
>
>
>
>extern FILE *fopen (__const char *__restrict __filename,
>                    __const char *__restrict __modes);
>
>
>
>
>extern FILE *freopen (__const char *__restrict __filename,
>                      __const char *__restrict __modes,
>                      FILE *__restrict __stream);
># 264 "/usr/include/stdio.h" 3
>
># 275 "/usr/include/stdio.h" 3
>extern FILE *fdopen (int __fd, __const char *__modes) ;
># 296 "/usr/include/stdio.h" 3
>
>
>
>extern void setbuf (FILE *__restrict __stream, char *__restrict __buf) ;
>
>
>
>extern int setvbuf (FILE *__restrict __stream, char *__restrict __buf,
>                    int __modes, size_t __n) ;
>
>
>
>
>
>extern void setbuffer (FILE *__restrict __stream, char *__restrict __buf,
>                       size_t __size) ;
>
>
>extern void setlinebuf (FILE *__stream) ;
>
>
>
>
>
>
>
>
>extern int fprintf (FILE *__restrict __stream,
>                    __const char *__restrict __format, ...);
>
>
>
>
>extern int printf (__const char *__restrict __format, ...);
>
>extern int sprintf (char *__restrict __s,
>                    __const char *__restrict __format, ...) ;
>
>
>
>
>
>extern int vfprintf (FILE *__restrict __s, __const char *__restrict __format,
>                     __gnuc_va_list __arg);
>
>
>
>
>extern int vprintf (__const char *__restrict __format, __gnuc_va_list __arg);
>
>extern int vsprintf (char *__restrict __s, __const char *__restrict __format,
>                     __gnuc_va_list __arg) ;
>
>
>
>
>
>extern int snprintf (char *__restrict __s, size_t __maxlen,
>                     __const char *__restrict __format, ...)
>     __attribute__ ((__format__ (__printf__, 3, 4)));
>
>extern int vsnprintf (char *__restrict __s, size_t __maxlen,
>                      __const char *__restrict __format, __gnuc_va_list __arg)
>     __attribute__ ((__format__ (__printf__, 3, 0)));
>
># 390 "/usr/include/stdio.h" 3
>
>
>
>
>
>extern int fscanf (FILE *__restrict __stream,
>                   __const char *__restrict __format, ...);
>
>
>
>
>extern int scanf (__const char *__restrict __format, ...);
>
>extern int sscanf (__const char *__restrict __s,
>                   __const char *__restrict __format, ...) ;
>
># 432 "/usr/include/stdio.h" 3
>
>
>
>
>
>extern int fgetc (FILE *__stream);
>extern int getc (FILE *__stream);
>
>
>
>
>
>extern int getchar (void);
>
># 456 "/usr/include/stdio.h" 3
>extern int getc_unlocked (FILE *__stream);
>extern int getchar_unlocked (void);
># 467 "/usr/include/stdio.h" 3
>extern int fgetc_unlocked (FILE *__stream);
>
>
>
>
>
>
>
>
>
>
>
>extern int fputc (int __c, FILE *__stream);
>extern int putc (int __c, FILE *__stream);
>
>
>
>
>
>extern int putchar (int __c);
>
># 500 "/usr/include/stdio.h" 3
>extern int fputc_unlocked (int __c, FILE *__stream);
>
>
>
>
>
>
>
>extern int putc_unlocked (int __c, FILE *__stream);
>extern int putchar_unlocked (int __c);
>
>
>
>
>
>
>extern int getw (FILE *__stream);
>
>
>extern int putw (int __w, FILE *__stream);
>
>
>
>
>
>
>
>
>extern char *fgets (char *__restrict __s, int __n, FILE *__restrict __stream);
>
>
>
>
>
>
>extern char *gets (char *__s);
>
># 580 "/usr/include/stdio.h" 3
>
>
>
>
>
>extern int fputs (__const char *__restrict __s, FILE *__restrict __stream);
>
>
>
>
>
>extern int puts (__const char *__s);
>
>
>
>
>
>
>extern int ungetc (int __c, FILE *__stream);
>
>
>
>
>
>
>extern size_t fread (void *__restrict __ptr, size_t __size,
>                     size_t __n, FILE *__restrict __stream);
>
>
>
>
>extern size_t fwrite (__const void *__restrict __ptr, size_t __size,
>                      size_t __n, FILE *__restrict __s);
>
># 633 "/usr/include/stdio.h" 3
>extern size_t fread_unlocked (void *__restrict __ptr, size_t __size,
>                              size_t __n, FILE *__restrict __stream);
>extern size_t fwrite_unlocked (__const void *__restrict __ptr, size_t __size,
>                               size_t __n, FILE *__restrict __stream);
>
>
>
>
>
>
>
>
>extern int fseek (FILE *__stream, long int __off, int __whence);
>
>
>
>
>extern long int ftell (FILE *__stream);
>
>
>
>
>extern void rewind (FILE *__stream);
>
># 688 "/usr/include/stdio.h" 3
>
>
>
>
>
>
>extern int fgetpos (FILE *__restrict __stream, fpos_t *__restrict __pos);
>
>
>
>
>extern int fsetpos (FILE *__stream, __const fpos_t *__pos);
># 711 "/usr/include/stdio.h" 3
>
># 720 "/usr/include/stdio.h" 3
>
>
>extern void clearerr (FILE *__stream) ;
>
>extern int feof (FILE *__stream) ;
>
>extern int ferror (FILE *__stream) ;
>
>
>
>
>extern void clearerr_unlocked (FILE *__stream) ;
>extern int feof_unlocked (FILE *__stream) ;
>extern int ferror_unlocked (FILE *__stream) ;
>
>
>
>
>
>
>
>
>extern void perror (__const char *__s);
>
>
>
>
>
>
># 1 "/usr/include/bits/sys_errlist.h" 1 3
># 27 "/usr/include/bits/sys_errlist.h" 3
>extern int sys_nerr;
>extern __const char *__const sys_errlist[];
># 750 "/usr/include/stdio.h" 2 3
>
>
>
>
>extern int fileno (FILE *__stream) ;
>
>
>
>
>extern int fileno_unlocked (FILE *__stream) ;
># 769 "/usr/include/stdio.h" 3
>extern FILE *popen (__const char *__command, __const char *__modes);
>
>
>
>
>
>extern int pclose (FILE *__stream);
>
>
>
>
>
>extern char *ctermid (char *__s) ;
># 809 "/usr/include/stdio.h" 3
>extern void flockfile (FILE *__stream) ;
>
>
>
>extern int ftrylockfile (FILE *__stream) ;
>
>
>extern void funlockfile (FILE *__stream) ;
># 830 "/usr/include/stdio.h" 3
># 1 "/usr/include/bits/stdio.h" 1 3
># 33 "/usr/include/bits/stdio.h" 3
>extern __inline int
>vprintf (__const char *__restrict __fmt, __gnuc_va_list __arg)
>{
>  return vfprintf (stdout, __fmt, __arg);
>}
>
>
>extern __inline int
>getchar (void)
>{
>  return _IO_getc (stdin);
>}
>
>
>
>
>extern __inline int
>getc_unlocked (FILE *__fp)
>{
>  return ((__fp)->_IO_read_ptr >= (__fp)->_IO_read_end ? __uflow (__fp) : *(unsigned char *) (__fp)->_IO_read_ptr++);
>}
>
>
>extern __inline int
>getchar_unlocked (void)
>{
>  return ((stdin)->_IO_read_ptr >= (stdin)->_IO_read_end ? __uflow (stdin) : *(unsigned char *) (stdin)->_IO_read_ptr++);
>}
>
>
>
>
>extern __inline int
>putchar (int __c)
>{
>  return _IO_putc (__c, stdout);
>}
>
>
>
>
>extern __inline int
>fputc_unlocked (int __c, FILE *__stream)
>{
>  return (((__stream)->_IO_write_ptr >= (__stream)->_IO_write_end) ? __overflow (__stream, (unsigned char) (__c)) : (unsigned char) (*(__stream)->_IO_write_ptr++ = (__c)));
>}
>
>
>
>
>
>extern __inline int
>putc_unlocked (int __c, FILE *__stream)
>{
>  return (((__stream)->_IO_write_ptr >= (__stream)->_IO_write_end) ? __overflow (__stream, (unsigned char) (__c)) : (unsigned char) (*(__stream)->_IO_write_ptr++ = (__c)));
>}
>
>
>extern __inline int
>putchar_unlocked (int __c)
>{
>  return (((stdout)->_IO_write_ptr >= (stdout)->_IO_write_end) ? __overflow (stdout, (unsigned char) (__c)) : (unsigned char) (*(stdout)->_IO_write_ptr++ = (__c)));
>}
># 111 "/usr/include/bits/stdio.h" 3
>extern __inline int
>feof_unlocked (FILE *__stream)
>{
>  return (((__stream)->_flags & 0x10) != 0);
>}
>
>
>extern __inline int
>ferror_unlocked (FILE *__stream)
>{
>  return (((__stream)->_flags & 0x20) != 0);
>}
># 831 "/usr/include/stdio.h" 2 3
>
>
>
># 20 "maxwell_eps.c" 2
># 1 "/usr/include/math.h" 1 3
># 29 "/usr/include/math.h" 3
>
>
>
>
># 1 "/usr/include/bits/huge_val.h" 1 3
># 34 "/usr/include/math.h" 2 3
>
>
>
>
>
>
># 1 "/usr/include/bits/mathdef.h" 1 3
># 41 "/usr/include/math.h" 2 3
># 65 "/usr/include/math.h" 3
># 1 "/usr/include/bits/mathcalls.h" 1 3
># 53 "/usr/include/bits/mathcalls.h" 3
>
>
>extern double acos (double __x) ; extern double __acos (double __x) ;
>
>extern double asin (double __x) ; extern double __asin (double __x) ;
>
>extern double atan (double __x) ; extern double __atan (double __x) ;
>
>extern double atan2 (double __y, double __x) ; extern double __atan2 (double __y, double __x) ;
>
>
>extern double cos (double __x) ; extern double __cos (double __x) ;
>
>extern double sin (double __x) ; extern double __sin (double __x) ;
>
>extern double tan (double __x) ; extern double __tan (double __x) ;
>
>
>
>
>extern double cosh (double __x) ; extern double __cosh (double __x) ;
>
>extern double sinh (double __x) ; extern double __sinh (double __x) ;
>
>extern double tanh (double __x) ; extern double __tanh (double __x) ;
>
># 87 "/usr/include/bits/mathcalls.h" 3
>
>
>extern double acosh (double __x) ; extern double __acosh (double __x) ;
>
>extern double asinh (double __x) ; extern double __asinh (double __x) ;
>
>extern double atanh (double __x) ; extern double __atanh (double __x) ;
>
>
>
>
>
>
>
>extern double exp (double __x) ; extern double __exp (double __x) ;
>
>
>extern double frexp (double __x, int *__exponent) ; extern double __frexp (double __x, int *__exponent) ;
>
>
>extern double ldexp (double __x, int __exponent) ; extern double __ldexp (double __x, int __exponent) ;
>
>
>extern double log (double __x) ; extern double __log (double __x) ;
>
>
>extern double log10 (double __x) ; extern double __log10 (double __x) ;
>
>
>extern double modf (double __x, double *__iptr) ; extern double __modf (double __x, double *__iptr) ;
>
># 127 "/usr/include/bits/mathcalls.h" 3
>
>
>extern double expm1 (double __x) ; extern double __expm1 (double __x) ;
>
>
>extern double log1p (double __x) ; extern double __log1p (double __x) ;
>
>
>extern double logb (double __x) ; extern double __logb (double __x) ;
>
># 152 "/usr/include/bits/mathcalls.h" 3
>
>
>extern double pow (double __x, double __y) ; extern double __pow (double __x, double __y) ;
>
>
>extern double sqrt (double __x) ; extern double __sqrt (double __x) ;
>
>
>
>
>
>extern double hypot (double __x, double __y) ; extern double __hypot (double __x, double __y) ;
>
>
>
>
>
>
>extern double cbrt (double __x) ; extern double __cbrt (double __x) ;
>
>
>
>
>
>
>
>
>extern double ceil (double __x) __attribute__ ((__const__)); extern double __ceil (double __x) __attribute__ ((__const__));
>
>
>extern double fabs (double __x) __attribute__ ((__const__)); extern double __fabs (double __x) __attribute__ ((__const__));
>
>
>extern double floor (double __x) __attribute__ ((__const__)); extern double __floor (double __x) __attribute__ ((__const__));
>
>
>extern double fmod (double __x, double __y) ; extern double __fmod (double __x, double __y) ;
>
>
>
>
>extern int __isinf (double __value) __attribute__ ((__const__));
>
>
>extern int __finite (double __value) __attribute__ ((__const__));
>
>
>
>
>
>extern int isinf (double __value) __attribute__ ((__const__));
>
>
>extern int finite (double __value) __attribute__ ((__const__));
>
>
>extern double drem (double __x, double __y) ; extern double __drem (double __x, double __y) ;
>
>
>
>extern double significand (double __x) ; extern double __significand (double __x) ;
>
>
>
>
>
>extern double copysign (double __x, double __y) __attribute__ ((__const__)); extern double __copysign (double __x, double __y) __attribute__ ((__const__));
>
># 231 "/usr/include/bits/mathcalls.h" 3
>extern int __isnan (double __value) __attribute__ ((__const__));
>
>
>
>extern int isnan (double __value) __attribute__ ((__const__));
>
>
>extern double j0 (double) ; extern double __j0 (double) ;
>extern double j1 (double) ; extern double __j1 (double) ;
>extern double jn (int, double) ; extern double __jn (int, double) ;
>extern double y0 (double) ; extern double __y0 (double) ;
>extern double y1 (double) ; extern double __y1 (double) ;
>extern double yn (int, double) ; extern double __yn (int, double) ;
>
>
>
>
>
>
>extern double erf (double) ; extern double __erf (double) ;
>extern double erfc (double) ; extern double __erfc (double) ;
>extern double lgamma (double) ; extern double __lgamma (double) ;
>
># 265 "/usr/include/bits/mathcalls.h" 3
>extern double gamma (double) ; extern double __gamma (double) ;
>
>
>
>
>
>
>extern double lgamma_r (double, int *__signgamp) ; extern double __lgamma_r (double, int *__signgamp) ;
>
>
>
>
>
>
>
>extern double rint (double __x) ; extern double __rint (double __x) ;
>
>
>extern double nextafter (double __x, double __y) __attribute__ ((__const__)); extern double __nextafter (double __x, double __y) __attribute__ ((__const__));
>
>
>
>
>
>extern double remainder (double __x, double __y) ; extern double __remainder (double __x, double __y) ;
>
>
>
>extern double scalbn (double __x, int __n) ; extern double __scalbn (double __x, int __n) ;
>
>
>
>extern int ilogb (double __x) ; extern int __ilogb (double __x) ;
># 66 "/usr/include/math.h" 2 3
># 88 "/usr/include/math.h" 3
># 1 "/usr/include/bits/mathcalls.h" 1 3
># 53 "/usr/include/bits/mathcalls.h" 3
>
>
>extern float acosf (float __x) ; extern float __acosf (float __x) ;
>
>extern float asinf (float __x) ; extern float __asinf (float __x) ;
>
>extern float atanf (float __x) ; extern float __atanf (float __x) ;
>
>extern float atan2f (float __y, float __x) ; extern float __atan2f (float __y, float __x) ;
>
>
>extern float cosf (float __x) ; extern float __cosf (float __x) ;
>
>extern float sinf (float __x) ; extern float __sinf (float __x) ;
>
>extern float tanf (float __x) ; extern float __tanf (float __x) ;
>
>
>
>
>extern float coshf (float __x) ; extern float __coshf (float __x) ;
>
>extern float sinhf (float __x) ; extern float __sinhf (float __x) ;
>
>extern float tanhf (float __x) ; extern float __tanhf (float __x) ;
>
># 87 "/usr/include/bits/mathcalls.h" 3
>
>
>extern float acoshf (float __x) ; extern float __acoshf (float __x) ;
>
>extern float asinhf (float __x) ; extern float __asinhf (float __x) ;
>
>extern float atanhf (float __x) ; extern float __atanhf (float __x) ;
>
>
>
>
>
>
>
>extern float expf (float __x) ; extern float __expf (float __x) ;
>
>
>extern float frexpf (float __x, int *__exponent) ; extern float __frexpf (float __x, int *__exponent) ;
>
>
>extern float ldexpf (float __x, int __exponent) ; extern float __ldexpf (float __x, int __exponent) ;
>
>
>extern float logf (float __x) ; extern float __logf (float __x) ;
>
>
>extern float log10f (float __x) ; extern float __log10f (float __x) ;
>
>
>extern float modff (float __x, float *__iptr) ; extern float __modff (float __x, float *__iptr) ;
>
># 127 "/usr/include/bits/mathcalls.h" 3
>
>
>extern float expm1f (float __x) ; extern float __expm1f (float __x) ;
>
>
>extern float log1pf (float __x) ; extern float __log1pf (float __x) ;
>
>
>extern float logbf (float __x) ; extern float __logbf (float __x) ;
>
># 152 "/usr/include/bits/mathcalls.h" 3
>
>
>extern float powf (float __x, float __y) ; extern float __powf (float __x, float __y) ;
>
>
>extern float sqrtf (float __x) ; extern float __sqrtf (float __x) ;
>
>
>
>
>
>extern float hypotf (float __x, float __y) ; extern float __hypotf (float __x, float __y) ;
>
>
>
>
>
>
>extern float cbrtf (float __x) ; extern float __cbrtf (float __x) ;
>
>
>
>
>
>
>
>
>extern float ceilf (float __x) __attribute__ ((__const__)); extern float __ceilf (float __x) __attribute__ ((__const__));
>
>
>extern float fabsf (float __x) __attribute__ ((__const__)); extern float __fabsf (float __x) __attribute__ ((__const__));
>
>
>extern float floorf (float __x) __attribute__ ((__const__)); extern float __floorf (float __x) __attribute__ ((__const__));
>
>
>extern float fmodf (float __x, float __y) ; extern float __fmodf (float __x, float __y) ;
>
>
>
>
>extern int __isinff (float __value) __attribute__ ((__const__));
>
>
>extern int __finitef (float __value) __attribute__ ((__const__));
>
>
>
>
>
>extern int isinff (float __value) __attribute__ ((__const__));
>
>
>extern int finitef (float __value) __attribute__ ((__const__));
>
>
>extern float dremf (float __x, float __y) ; extern float __dremf (float __x, float __y) ;
>
>
>
>extern float significandf (float __x) ; extern float __significandf (float __x) ;
>
>
>
>
>
>extern float copysignf (float __x, float __y) __attribute__ ((__const__)); extern float __copysignf (float __x, float __y) __attribute__ ((__const__));
>
># 231 "/usr/include/bits/mathcalls.h" 3
>extern int __isnanf (float __value) __attribute__ ((__const__));
>
>
>
>extern int isnanf (float __value) __attribute__ ((__const__));
>
>
>extern float j0f (float) ; extern float __j0f (float) ;
>extern float j1f (float) ; extern float __j1f (float) ;
>extern float jnf (int, float) ; extern float __jnf (int, float) ;
>extern float y0f (float) ; extern float __y0f (float) ;
>extern float y1f (float) ; extern float __y1f (float) ;
>extern float ynf (int, float) ; extern float __ynf (int, float) ;
>
>
>
>
>
>
>extern float erff (float) ; extern float __erff (float) ;
>extern float erfcf (float) ; extern float __erfcf (float) ;
>extern float lgammaf (float) ; extern float __lgammaf (float) ;
>
># 265 "/usr/include/bits/mathcalls.h" 3
>extern float gammaf (float) ; extern float __gammaf (float) ;
>
>
>
>
>
>
>extern float lgammaf_r (float, int *__signgamp) ; extern float __lgammaf_r (float, int *__signgamp) ;
>
>
>
>
>
>
>
>extern float rintf (float __x) ; extern float __rintf (float __x) ;
>
>
>extern float nextafterf (float __x, float __y) __attribute__ ((__const__)); extern float __nextafterf (float __x, float __y) __attribute__ ((__const__));
>
>
>
>
>
>extern float remainderf (float __x, float __y) ; extern float __remainderf (float __x, float __y) ;
>
>
>
>extern float scalbnf (float __x, int __n) ; extern float __scalbnf (float __x, int __n) ;
>
>
>
>extern int ilogbf (float __x) ; extern int __ilogbf (float __x) ;
># 89 "/usr/include/math.h" 2 3
># 109 "/usr/include/math.h" 3
># 1 "/usr/include/bits/mathcalls.h" 1 3
># 53 "/usr/include/bits/mathcalls.h" 3
>
>
>extern long double acosl (long double __x) ; extern long double __acosl (long double __x) ;
>
>extern long double asinl (long double __x) ; extern long double __asinl (long double __x) ;
>
>extern long double atanl (long double __x) ; extern long double __atanl (long double __x) ;
>
>extern long double atan2l (long double __y, long double __x) ; extern long double __atan2l (long double __y, long double __x) ;
>
>
>extern long double cosl (long double __x) ; extern long double __cosl (long double __x) ;
>
>extern long double sinl (long double __x) ; extern long double __sinl (long double __x) ;
>
>extern long double tanl (long double __x) ; extern long double __tanl (long double __x) ;
>
>
>
>
>extern long double coshl (long double __x) ; extern long double __coshl (long double __x) ;
>
>extern long double sinhl (long double __x) ; extern long double __sinhl (long double __x) ;
>
>extern long double tanhl (long double __x) ; extern long double __tanhl (long double __x) ;
>
># 87 "/usr/include/bits/mathcalls.h" 3
>
>
>extern long double acoshl (long double __x) ; extern long double __acoshl (long double __x) ;
>
>extern long double asinhl (long double __x) ; extern long double __asinhl (long double __x) ;
>
>extern long double atanhl (long double __x) ; extern long double __atanhl (long double __x) ;
>
>
>
>
>
>
>
>extern long double expl (long double __x) ; extern long double __expl (long double __x) ;
>
>
>extern long double frexpl (long double __x, int *__exponent) ; extern long double __frexpl (long double __x, int *__exponent) ;
>
>
>extern long double ldexpl (long double __x, int __exponent) ; extern long double __ldexpl (long double __x, int __exponent) ;
>
>
>extern long double logl (long double __x) ; extern long double __logl (long double __x) ;
>
>
>extern long double log10l (long double __x) ; extern long double __log10l (long double __x) ;
>
>
>extern long double modfl (long double __x, long double *__iptr) ; extern long double __modfl (long double __x, long double *__iptr) ;
>
># 127 "/usr/include/bits/mathcalls.h" 3
>
>
>extern long double expm1l (long double __x) ; extern long double __expm1l (long double __x) ;
>
>
>extern long double log1pl (long double __x) ; extern long double __log1pl (long double __x) ;
>
>
>extern long double logbl (long double __x) ; extern long double __logbl (long double __x) ;
>
># 152 "/usr/include/bits/mathcalls.h" 3
>
>
>extern long double powl (long double __x, long double __y) ; extern long double __powl (long double __x, long double __y) ;
>
>
>extern long double sqrtl (long double __x) ; extern long double __sqrtl (long double __x) ;
>
>
>
>
>
>extern long double hypotl (long double __x, long double __y) ; extern long double __hypotl (long double __x, long double __y) ;
>
>
>
>
>
>
>extern long double cbrtl (long double __x) ; extern long double __cbrtl (long double __x) ;
>
>
>
>
>
>
>
>
>extern long double ceill (long double __x) __attribute__ ((__const__)); extern long double __ceill (long double __x) __attribute__ ((__const__));
>
>
>extern long double fabsl (long double __x) __attribute__ ((__const__)); extern long double __fabsl (long double __x) __attribute__ ((__const__));
>
>
>extern long double floorl (long double __x) __attribute__ ((__const__)); extern long double __floorl (long double __x) __attribute__ ((__const__));
>
>
>extern long double fmodl (long double __x, long double __y) ; extern long double __fmodl (long double __x, long double __y) ;
>
>
>
>
>extern int __isinfl (long double __value) __attribute__ ((__const__));
>
>
>extern int __finitel (long double __value) __attribute__ ((__const__));
>
>
>
>
>
>extern int isinfl (long double __value) __attribute__ ((__const__));
>
>
>extern int finitel (long double __value) __attribute__ ((__const__));
>
>
>extern long double dreml (long double __x, long double __y) ; extern long double __dreml (long double __x, long double __y) ;
>
>
>
>extern long double significandl (long double __x) ; extern long double __significandl (long double __x) ;
>
>
>
>
>
>extern long double copysignl (long double __x, long double __y) __attribute__ ((__const__)); extern long double __copysignl (long double __x, long double __y) __attribute__ ((__const__));
>
># 231 "/usr/include/bits/mathcalls.h" 3
>extern int __isnanl (long double __value) __attribute__ ((__const__));
>
>
>
>extern int isnanl (long double __value) __attribute__ ((__const__));
>
>
>extern long double j0l (long double) ; extern long double __j0l (long double) ;
>extern long double j1l (long double) ; extern long double __j1l (long double) ;
>extern long double jnl (int, long double) ; extern long double __jnl (int, long double) ;
>extern long double y0l (long double) ; extern long double __y0l (long double) ;
>extern long double y1l (long double) ; extern long double __y1l (long double) ;
>extern long double ynl (int, long double) ; extern long double __ynl (int, long double) ;
>
>
>
>
>
>
>extern long double erfl (long double) ; extern long double __erfl (long double) ;
>extern long double erfcl (long double) ; extern long double __erfcl (long double) ;
>extern long double lgammal (long double) ; extern long double __lgammal (long double) ;
>
># 265 "/usr/include/bits/mathcalls.h" 3
>extern long double gammal (long double) ; extern long double __gammal (long double) ;
>
>
>
>
>
>
>extern long double lgammal_r (long double, int *__signgamp) ; extern long double __lgammal_r (long double, int *__signgamp) ;
>
>
>
>
>
>
>
>extern long double rintl (long double __x) ; extern long double __rintl (long double __x) ;
>
>
>extern long double nextafterl (long double __x, long double __y) __attribute__ ((__const__)); extern long double __nextafterl (long double __x, long double __y) __attribute__ ((__const__));
>
>
>
>
>
>extern long double remainderl (long double __x, long double __y) ; extern long double __remainderl (long double __x, long double __y) ;
>
>
>
>extern long double scalbnl (long double __x, int __n) ; extern long double __scalbnl (long double __x, int __n) ;
>
>
>
>extern int ilogbl (long double __x) ; extern int __ilogbl (long double __x) ;
># 110 "/usr/include/math.h" 2 3
># 125 "/usr/include/math.h" 3
>extern int signgam;
># 252 "/usr/include/math.h" 3
>typedef enum
>{
>  _IEEE_ = -1,
>  _SVID_,
>  _XOPEN_,
>  _POSIX_,
>  _ISOC_
>} _LIB_VERSION_TYPE;
>
>
>
>
>extern _LIB_VERSION_TYPE _LIB_VERSION;
># 277 "/usr/include/math.h" 3
>struct exception
>
>  {
>    int type;
>    char *name;
>    double arg1;
>    double arg2;
>    double retval;
>  };
>
>
>
>
>extern int matherr (struct exception *__exc);
># 362 "/usr/include/math.h" 3
># 1 "/usr/include/bits/mathinline.h" 1 3
># 25 "/usr/include/bits/mathinline.h" 3
># 1 "/usr/include/bits/wordsize.h" 1 3
># 26 "/usr/include/bits/mathinline.h" 2 3
># 363 "/usr/include/math.h" 2 3
># 423 "/usr/include/math.h" 3
>
># 21 "maxwell_eps.c" 2
>
># 1 "../config.h" 1
># 23 "maxwell_eps.c" 2
># 1 "../util/check.h" 1
># 27 "../util/check.h"
># 1 "../util/mpi_utils.h" 1
># 23 "../util/mpi_utils.h"
>extern void mpi_die(const char *template, ...)
>
>     __attribute__ ((format (printf, 1, 2)))
>
>;
>
>extern void mpi_one_fprintf(FILE *f, const char *template, ...)
>
>     __attribute__ ((format (printf, 2, 3)))
>
>;
>
>extern void mpi_one_printf(const char *template, ...)
>
>     __attribute__ ((format (printf, 1, 2)))
>
>;
>
>extern int mpi_is_master(void);
>
>extern void mpi_assert_equal(double x);
>
>extern void mpi_begin_critical_section(int tag);
>extern void mpi_end_critical_section(int tag);
># 28 "../util/check.h" 2
># 46 "../util/check.h"
>extern void* (*my_malloc_hook) (size_t);
># 59 "../util/check.h"
>extern void debug_output_malloc_count(void);
>extern void debug_check_memory_leaks(void);
># 24 "maxwell_eps.c" 2
># 1 "../util/mpiglue.h" 1
># 38 "../util/mpiglue.h"
># 1 "/usr/include/string.h" 1 3
># 28 "/usr/include/string.h" 3
>
>
>
>
>
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 34 "/usr/include/string.h" 2 3
>
>
>
>
>extern void *memcpy (void *__restrict __dest,
>                     __const void *__restrict __src, size_t __n) ;
>
>
>extern void *memmove (void *__dest, __const void *__src, size_t __n)
>     ;
>
>
>
>
>
>
>extern void *memccpy (void *__restrict __dest, __const void *__restrict __src,
>                      int __c, size_t __n)
>     ;
>
>
>
>
>
>extern void *memset (void *__s, int __c, size_t __n) ;
>
>
>extern int memcmp (__const void *__s1, __const void *__s2, size_t __n)
>     __attribute__ ((__pure__));
>
>
>extern void *memchr (__const void *__s, int __c, size_t __n)
>      __attribute__ ((__pure__));
>
># 80 "/usr/include/string.h" 3
>
>
>extern char *strcpy (char *__restrict __dest, __const char *__restrict __src)
>     ;
>
>extern char *strncpy (char *__restrict __dest,
>                      __const char *__restrict __src, size_t __n) ;
>
>
>extern char *strcat (char *__restrict __dest, __const char *__restrict __src)
>     ;
>
>extern char *strncat (char *__restrict __dest, __const char *__restrict __src,
>                      size_t __n) ;
>
>
>extern int strcmp (__const char *__s1, __const char *__s2)
>     __attribute__ ((__pure__));
>
>extern int strncmp (__const char *__s1, __const char *__s2, size_t __n)
>     __attribute__ ((__pure__));
>
>
>extern int strcoll (__const char *__s1, __const char *__s2)
>     __attribute__ ((__pure__));
>
>extern size_t strxfrm (char *__restrict __dest,
>                       __const char *__restrict __src, size_t __n) ;
>
># 126 "/usr/include/string.h" 3
>extern char *strdup (__const char *__s) __attribute__ ((__malloc__));
># 160 "/usr/include/string.h" 3
>
>
>extern char *strchr (__const char *__s, int __c) __attribute__ ((__pure__));
>
>extern char *strrchr (__const char *__s, int __c) __attribute__ ((__pure__));
>
>
>
>
>
>
>
>
>
>
>
>extern size_t strcspn (__const char *__s, __const char *__reject)
>     __attribute__ ((__pure__));
>
>
>extern size_t strspn (__const char *__s, __const char *__accept)
>     __attribute__ ((__pure__));
>
>extern char *strpbrk (__const char *__s, __const char *__accept)
>     __attribute__ ((__pure__));
>
>extern char *strstr (__const char *__haystack, __const char *__needle)
>     __attribute__ ((__pure__));
>
>
>
>extern char *strtok (char *__restrict __s, __const char *__restrict __delim)
>     ;
>
>
>
>
>extern char *__strtok_r (char *__restrict __s,
>                         __const char *__restrict __delim,
>                         char **__restrict __save_ptr) ;
>
>extern char *strtok_r (char *__restrict __s, __const char *__restrict __delim,
>                       char **__restrict __save_ptr) ;
># 228 "/usr/include/string.h" 3
>
>
>extern size_t strlen (__const char *__s) __attribute__ ((__pure__));
>
># 241 "/usr/include/string.h" 3
>
>
>extern char *strerror (int __errnum) ;
>
>
>
>
>extern char *strerror_r (int __errnum, char *__buf, size_t __buflen) ;
>
>
>
>
>extern void __bzero (void *__s, size_t __n) ;
>
>
>
>extern void bcopy (__const void *__src, void *__dest, size_t __n) ;
>
>
>extern void bzero (void *__s, size_t __n) ;
>
>
>extern int bcmp (__const void *__s1, __const void *__s2, size_t __n)
>     __attribute__ ((__pure__));
>
>
>extern char *index (__const char *__s, int __c) __attribute__ ((__pure__));
>
>
>extern char *rindex (__const char *__s, int __c) __attribute__ ((__pure__));
>
>
>
>extern int ffs (int __i) __attribute__ ((__const__));
># 287 "/usr/include/string.h" 3
>extern int strcasecmp (__const char *__s1, __const char *__s2)
>     __attribute__ ((__pure__));
>
>
>extern int strncasecmp (__const char *__s1, __const char *__s2, size_t __n)
>     __attribute__ ((__pure__));
># 309 "/usr/include/string.h" 3
>extern char *strsep (char **__restrict __stringp,
>                     __const char *__restrict __delim) ;
># 372 "/usr/include/string.h" 3
># 1 "/usr/include/bits/string.h" 1 3
># 373 "/usr/include/string.h" 2 3
>
>
># 1 "/usr/include/bits/string2.h" 1 3
># 389 "/usr/include/bits/string2.h" 3
>extern void *__rawmemchr (const void *__s, int __c);
># 919 "/usr/include/bits/string2.h" 3
>extern __inline size_t __strcspn_c1 (__const char *__s, int __reject);
>extern __inline size_t
>__strcspn_c1 (__const char *__s, int __reject)
>{
>  register size_t __result = 0;
>  while (__s[__result] != '\0' && __s[__result] != __reject)
>    ++__result;
>  return __result;
>}
>
>extern __inline size_t __strcspn_c2 (__const char *__s, int __reject1,
>                                     int __reject2);
>extern __inline size_t
>__strcspn_c2 (__const char *__s, int __reject1, int __reject2)
>{
>  register size_t __result = 0;
>  while (__s[__result] != '\0' && __s[__result] != __reject1
>         && __s[__result] != __reject2)
>    ++__result;
>  return __result;
>}
>
>extern __inline size_t __strcspn_c3 (__const char *__s, int __reject1,
>                                     int __reject2, int __reject3);
>extern __inline size_t
>__strcspn_c3 (__const char *__s, int __reject1, int __reject2,
>              int __reject3)
>{
>  register size_t __result = 0;
>  while (__s[__result] != '\0' && __s[__result] != __reject1
>         && __s[__result] != __reject2 && __s[__result] != __reject3)
>    ++__result;
>  return __result;
>}
># 976 "/usr/include/bits/string2.h" 3
>extern __inline size_t __strspn_c1 (__const char *__s, int __accept);
>extern __inline size_t
>__strspn_c1 (__const char *__s, int __accept)
>{
>  register size_t __result = 0;
>
>  while (__s[__result] == __accept)
>    ++__result;
>  return __result;
>}
>
>extern __inline size_t __strspn_c2 (__const char *__s, int __accept1,
>                                    int __accept2);
>extern __inline size_t
>__strspn_c2 (__const char *__s, int __accept1, int __accept2)
>{
>  register size_t __result = 0;
>
>  while (__s[__result] == __accept1 || __s[__result] == __accept2)
>    ++__result;
>  return __result;
>}
>
>extern __inline size_t __strspn_c3 (__const char *__s, int __accept1,
>                                    int __accept2, int __accept3);
>extern __inline size_t
>__strspn_c3 (__const char *__s, int __accept1, int __accept2, int __accept3)
>{
>  register size_t __result = 0;
>
>  while (__s[__result] == __accept1 || __s[__result] == __accept2
>         || __s[__result] == __accept3)
>    ++__result;
>  return __result;
>}
># 1033 "/usr/include/bits/string2.h" 3
>extern __inline char *__strpbrk_c2 (__const char *__s, int __accept1,
>                                     int __accept2);
>extern __inline char *
>__strpbrk_c2 (__const char *__s, int __accept1, int __accept2)
>{
>
>  while (*__s != '\0' && *__s != __accept1 && *__s != __accept2)
>    ++__s;
>  return *__s == '\0' ? ((void *)0) : (char *) (size_t) __s;
>}
>
>extern __inline char *__strpbrk_c3 (__const char *__s, int __accept1,
>                                     int __accept2, int __accept3);
>extern __inline char *
>__strpbrk_c3 (__const char *__s, int __accept1, int __accept2,
>              int __accept3)
>{
>
>  while (*__s != '\0' && *__s != __accept1 && *__s != __accept2
>         && *__s != __accept3)
>    ++__s;
>  return *__s == '\0' ? ((void *)0) : (char *) (size_t) __s;
>}
># 1085 "/usr/include/bits/string2.h" 3
>extern __inline char *__strtok_r_1c (char *__s, char __sep, char **__nextp);
>extern __inline char *
>__strtok_r_1c (char *__s, char __sep, char **__nextp)
>{
>  char *__result;
>  if (__s == ((void *)0))
>    __s = *__nextp;
>  while (*__s == __sep)
>    ++__s;
>  __result = ((void *)0);
>  if (*__s != '\0')
>    {
>      __result = __s++;
>      while (*__s != '\0')
>        if (*__s++ == __sep)
>          {
>            __s[-1] = '\0';
>            break;
>          }
>      *__nextp = __s;
>    }
>  return __result;
>}
># 1117 "/usr/include/bits/string2.h" 3
>extern char *__strsep_g (char **__stringp, __const char *__delim);
># 1135 "/usr/include/bits/string2.h" 3
>extern __inline char *__strsep_1c (char **__s, char __reject);
>extern __inline char *
>__strsep_1c (char **__s, char __reject)
>{
>  register char *__retval = *__s;
>  if (__retval != ((void *)0) && (*__s = (__extension__ (__builtin_constant_p (__reject) && (__reject) == '\0' ? (char *) __rawmemchr (__retval, __reject) : strchr (__retval, __reject)))) != ((void *)0))
>    *(*__s)++ = '\0';
>  return __retval;
>}
>
>extern __inline char *__strsep_2c (char **__s, char __reject1, char __reject2);
>extern __inline char *
>__strsep_2c (char **__s, char __reject1, char __reject2)
>{
>  register char *__retval = *__s;
>  if (__retval != ((void *)0))
>    {
>      register char *__cp = __retval;
>      while (1)
>        {
>          if (*__cp == '\0')
>            {
>              __cp = ((void *)0);
>          break;
>            }
>          if (*__cp == __reject1 || *__cp == __reject2)
>            {
>              *__cp++ = '\0';
>              break;
>            }
>          ++__cp;
>        }
>      *__s = __cp;
>    }
>  return __retval;
>}
>
>extern __inline char *__strsep_3c (char **__s, char __reject1, char __reject2,
>                                   char __reject3);
>extern __inline char *
>__strsep_3c (char **__s, char __reject1, char __reject2, char __reject3)
>{
>  register char *__retval = *__s;
>  if (__retval != ((void *)0))
>    {
>      register char *__cp = __retval;
>      while (1)
>        {
>          if (*__cp == '\0')
>            {
>              __cp = ((void *)0);
>          break;
>            }
>          if (*__cp == __reject1 || *__cp == __reject2 || *__cp == __reject3)
>            {
>              *__cp++ = '\0';
>              break;
>            }
>          ++__cp;
>        }
>      *__s = __cp;
>    }
>  return __retval;
>}
># 1216 "/usr/include/bits/string2.h" 3
>extern char *__strdup (__const char *__string) __attribute__ ((__malloc__));
># 1235 "/usr/include/bits/string2.h" 3
>extern char *__strndup (__const char *__string, size_t __n)
>     __attribute__ ((__malloc__));
># 376 "/usr/include/string.h" 2 3
>
>
>
>
># 39 "../util/mpiglue.h" 2
># 1 "/usr/include/time.h" 1 3
># 30 "/usr/include/time.h" 3
>
>
>
>
>
>
>
>
># 1 "/usr/lib/gcc-lib/x86_64-redhat-linux/3.2.3/include/stddef.h" 1 3
># 39 "/usr/include/time.h" 2 3
>
>
>
># 1 "/usr/include/bits/time.h" 1 3
># 40 "/usr/include/bits/time.h" 3
>extern long int __sysconf (int);
># 43 "/usr/include/time.h" 2 3
># 58 "/usr/include/time.h" 3
>
>
>typedef __clock_t clock_t;
>
>
>
># 129 "/usr/include/time.h" 3
>
>
>struct tm
>{
>  int tm_sec;
>  int tm_min;
>  int tm_hour;
>  int tm_mday;
>  int tm_mon;
>  int tm_year;
>  int tm_wday;
>  int tm_yday;
>  int tm_isdst;
>
>
>  long int tm_gmtoff;
>  __const char *tm_zone;
>
>
>
>
>};
>
>
>
>
>
>
>
>
>struct itimerspec
>  {
>    struct timespec it_interval;
>    struct timespec it_value;
>  };
>
>
>struct sigevent;
># 178 "/usr/include/time.h" 3
>
>
>
>extern clock_t clock (void) ;
>
>
>extern time_t time (time_t *__timer) ;
>
>
>extern double difftime (time_t __time1, time_t __time0)
>     __attribute__ ((__const__));
>
>
>extern time_t mktime (struct tm *__tp) ;
>
>
>
>
>
>extern size_t strftime (char *__restrict __s, size_t __maxsize,
>                        __const char *__restrict __format,
>                        __const struct tm *__restrict __tp) ;
>
># 226 "/usr/include/time.h" 3
>
>
>
>extern struct tm *gmtime (__const time_t *__timer) ;
>
>
>
>extern struct tm *localtime (__const time_t *__timer) ;
>
>
>
>
>
>extern struct tm *gmtime_r (__const time_t *__restrict __timer,
>                            struct tm *__restrict __tp) ;
>
>
>
>extern struct tm *localtime_r (__const time_t *__restrict __timer,
>                               struct tm *__restrict __tp) ;
>
>
>
>
>
>extern char *asctime (__const struct tm *__tp) ;
>
>
>extern char *ctime (__const time_t *__timer) ;
>
>
>
>
>
>
>
>extern char *asctime_r (__const struct tm *__restrict __tp,
>                        char *__restrict __buf) ;
>
>
>extern char *ctime_r (__const time_t *__restrict __timer,
>                      char *__restrict __buf) ;
>
>
>
>
>extern char *__tzname[2];
>extern int __daylight;
>extern long int __timezone;
>
>
>
>
>extern char *tzname[2];
>
>
>
>extern void tzset (void) ;
>
>
>
>extern int daylight;
>extern long int timezone;
>
>
>
>
>
>extern int stime (__const time_t *__when) ;
># 309 "/usr/include/time.h" 3
>extern time_t timegm (struct tm *__tp) ;
>
>
>extern time_t timelocal (struct tm *__tp) ;
>
>
>extern int dysize (int __year) __attribute__ ((__const__));
># 324 "/usr/include/time.h" 3
>extern int nanosleep (__const struct timespec *__requested_time,
>                      struct timespec *__remaining);
>
>
>
>extern int clock_getres (clockid_t __clock_id, struct timespec *__res) ;
>
>
>extern int clock_gettime (clockid_t __clock_id, struct timespec *__tp) ;
>
>
>extern int clock_settime (clockid_t __clock_id, __const struct timespec *__tp)
>     ;
># 353 "/usr/include/time.h" 3
>extern int timer_create (clockid_t __clock_id,
>                         struct sigevent *__restrict __evp,
>                         timer_t *__restrict __timerid) ;
>
>
>extern int timer_delete (timer_t __timerid) ;
>
>
>extern int timer_settime (timer_t __timerid, int __flags,
>                          __const struct itimerspec *__restrict __value,
>                          struct itimerspec *__restrict __ovalue) ;
>
>
>extern int timer_gettime (timer_t __timerid, struct itimerspec *__value)
>     ;
>
>
>extern int timer_getoverrun (timer_t __timerid) ;
># 413 "/usr/include/time.h" 3
>
># 40 "../util/mpiglue.h" 2
># 73 "../util/mpiglue.h"
>typedef int MPI_Status;
>
>typedef clock_t mpiglue_clock_t;
># 25 "maxwell_eps.c" 2
>
>
># 1 "maxwell.h" 1
># 21 "maxwell.h"
># 1 "../matrices/scalar.h" 1
># 25 "../matrices/scalar.h"
>typedef double real;
>
>
>
>
>typedef struct {
>     real re, im;
>} scalar_complex;
># 110 "../matrices/scalar.h"
>typedef scalar_complex scalar;
># 22 "maxwell.h" 2
># 1 "../matrices/matrices.h" 1
># 23 "../matrices/matrices.h"
>typedef struct {
>     int N, localN, Nstart, allocN;
>     int c;
>     int n, p, alloc_p;
>     scalar *data;
>} evectmatrix;
>
>typedef struct {
>     int p, alloc_p;
>     scalar *data;
>} sqmatrix;
>
>
>extern double evectmatrix_flops;
>
>
>
>extern evectmatrix create_evectmatrix(int N, int c, int p,
>                                      int localN, int Nstart, int allocN);
>extern void destroy_evectmatrix(evectmatrix X);
>extern sqmatrix create_sqmatrix(int p);
>extern void destroy_sqmatrix(sqmatrix X);
>
>
>
>extern void matrix_XtY_diag(scalar *X, scalar *Y, int n, int p, scalar *diag);
>extern void matrix_XtY_diag_real(scalar *X, scalar *Y, int n, int p,
>                                 real *diag);
>extern void matrix_XtX_diag_real(scalar *X, int n, int p, real *diag);
>extern void matrix_XpaY_diag(scalar *X, real a, scalar *Y,
>                             scalar *diag, int n, int p);
>extern void matrix_XpaY_diag_real(scalar *X, real a, scalar *Y,
>                                  real *diag, int n, int p);
>extern void matrix_X_diag_real_pY_diag_real(scalar *X, real *diag1,
>                                            scalar *Y, real *diag2,
>                                            int n, int p);
>extern real matrix_re_trace_A_diag_real(scalar *A, real *diag, int p);
>extern scalar matrix_diag_trace(scalar *diag, int p);
>extern real matrix_diag_real_trace(real *diag, int p);
>
>
>
>extern void evectmatrix_resize(evectmatrix *A, int p, short preserve_data);
>extern void evectmatrix_copy(evectmatrix X, evectmatrix Y);
>extern void evectmatrix_copy_slice(evectmatrix X, evectmatrix Y,
>                                   int ix, int iy, int p);
>extern void evectmatrix_aXpbY(real a, evectmatrix X, real b, evectmatrix Y);
>extern void evectmatrix_aXpbYS_sub(real a, evectmatrix X,
>                                   real b, evectmatrix Y,
>                                   sqmatrix S, int Soffset, short sdagger);
>extern void evectmatrix_XeYS(evectmatrix X, evectmatrix Y,
>                             sqmatrix S, short sherm);
>extern void evectmatrix_XpaYS(evectmatrix X, real a, evectmatrix Y,
>                              sqmatrix S, short sdagger);
>extern void evectmatrix_XtX(sqmatrix U, evectmatrix X, sqmatrix S);
>extern void evectmatrix_XtY(sqmatrix U, evectmatrix X, evectmatrix Y,
>                            sqmatrix S);
>extern void evectmatrix_XtY_slice(sqmatrix U, evectmatrix X, evectmatrix Y,
>                                  int ix, int iy, int p, sqmatrix S);
>extern void evectmatrixXtY_sub(sqmatrix U, int Uoffset,
>                               evectmatrix X, evectmatrix Y, sqmatrix S);
>extern void evectmatrix_XtY_diag(evectmatrix X, evectmatrix Y, scalar *diag,
>                                 scalar *scratch_diag);
>extern void evectmatrix_XtY_diag_real(evectmatrix X, evectmatrix Y,
>                                      real *diag, real *scratch_diag);
>extern void evectmatrix_XtX_diag_real(evectmatrix X, real *diag,
>                                      real *scratch_diag);
>extern scalar evectmatrix_traceXtY(evectmatrix X, evectmatrix Y);
>
>
>
>extern void sqmatrix_copy(sqmatrix A, sqmatrix B);
>extern void sqmatrix_resize(sqmatrix *A, int p, short preserve_data);
>extern void sqmatrix_copy_upper2full(sqmatrix F, sqmatrix U);
>extern void sqmatrix_symmetrize(sqmatrix Asym, sqmatrix A);
>extern scalar sqmatrix_trace(sqmatrix U);
>extern scalar sqmatrix_traceAtB(sqmatrix A, sqmatrix B);
>extern void sqmatrix_AeBC(sqmatrix A, sqmatrix B, short bdagger,
>                          sqmatrix C, short cdagger);
>extern void sqmatrix_ApaBC(sqmatrix A, real a, sqmatrix B, short bdagger,
>                           sqmatrix C, short cdagger);
>extern void sqmatrix_ApaB(sqmatrix A, real a, sqmatrix B);
>extern void sqmatrix_aApbB(real a, sqmatrix A, real b, sqmatrix B);
>extern void sqmatrix_invert(sqmatrix U, short positive_definite,
>                            sqmatrix Work);
>extern void sqmatrix_eigensolve(sqmatrix U, real *eigenvals, sqmatrix W);
>extern void sqmatrix_sqrt(sqmatrix Usqrt, sqmatrix U, sqmatrix W);
># 23 "maxwell.h" 2
>
>
># 1 "/software/mpbi/x86_64/install/fftw-2.1.5/include/fftw.h" 1
># 46 "/software/mpbi/x86_64/install/fftw-2.1.5/include/fftw.h"
>typedef double fftw_real;
>
>
>
>
>
>typedef struct {
>     fftw_real re, im;
>} fftw_complex;
>
>
>
>typedef enum {
>     FFTW_FORWARD = -1, FFTW_BACKWARD = 1
>} fftw_direction;
>
>
>typedef fftw_complex FFTW_COMPLEX;
>typedef fftw_real FFTW_REAL;
># 80 "/software/mpbi/x86_64/install/fftw-2.1.5/include/fftw.h"
>typedef enum {
>     FFTW_SUCCESS = 0, FFTW_FAILURE = -1
>} fftw_status;
>
>
>
>
>typedef void (fftw_notw_codelet)
>     (const fftw_complex *, fftw_complex *, int, int);
>typedef void (fftw_twiddle_codelet)
>     (fftw_complex *, const fftw_complex *, int,
>      int, int);
>typedef void (fftw_generic_codelet)
>     (fftw_complex *, const fftw_complex *, int,
>      int, int, int);
>typedef void (fftw_real2hc_codelet)
>     (const fftw_real *, fftw_real *, fftw_real *,
>      int, int, int);
>typedef void (fftw_hc2real_codelet)
>     (const fftw_real *, const fftw_real *,
>      fftw_real *, int, int, int);
>typedef void (fftw_hc2hc_codelet)
>     (fftw_real *, const fftw_complex *,
>      int, int, int);
>typedef void (fftw_rgeneric_codelet)
>     (fftw_real *, const fftw_complex *, int,
>      int, int, int);
># 115 "/software/mpbi/x86_64/install/fftw-2.1.5/include/fftw.h"
>enum fftw_node_type {
>     FFTW_NOTW, FFTW_TWIDDLE, FFTW_GENERIC, FFTW_RADER,
>     FFTW_REAL2HC, FFTW_HC2REAL, FFTW_HC2HC, FFTW_RGENERIC
>};
>
>
>typedef struct {
>     const char *name;
>     void (*codelet) ();
>     int size;
>     fftw_direction dir;
>     enum fftw_node_type type;
>     int signature;
>     int ntwiddle;
>     const int *twiddle_order;
>
>
>
>
>} fftw_codelet_desc;
># 150 "/software/mpbi/x86_64/install/fftw-2.1.5/include/fftw.h"
>extern const char * fftw_version;
># 165 "/software/mpbi/x86_64/install/fftw-2.1.5/include/fftw.h"
>typedef struct fftw_twiddle_struct {
>     int n;
>     const fftw_codelet_desc *cdesc;
>     fftw_complex *twarray;
>     struct fftw_twiddle_struct *next;
>     int refcnt;
>} fftw_twiddle;
>
>typedef struct fftw_rader_data_struct {
>     struct fftw_plan_struct *plan;
>     fftw_complex *omega;
>     int g, ginv;
>     int p, flags, refcount;
>     struct fftw_rader_data_struct *next;
>     fftw_codelet_desc *cdesc;
>} fftw_rader_data;
>
>typedef void (fftw_rader_codelet)
>     (fftw_complex *, const fftw_complex *, int,
>      int, int, fftw_rader_data *);
>
>
>typedef struct fftw_plan_node_struct {
>     enum fftw_node_type type;
>
>     union {
>
>          struct {
>               int size;
>               fftw_notw_codelet *codelet;
>               const fftw_codelet_desc *codelet_desc;
>          } notw;
>
>
>          struct {
>               int size;
>               fftw_twiddle_codelet *codelet;
>               fftw_twiddle *tw;
>               struct fftw_plan_node_struct *recurse;
>               const fftw_codelet_desc *codelet_desc;
>          } twiddle;
>
>
>          struct {
>               int size;
>               fftw_generic_codelet *codelet;
>               fftw_twiddle *tw;
>               struct fftw_plan_node_struct *recurse;
>          } generic;
>
>
>          struct {
>               int size;
>               fftw_rader_codelet *codelet;
>               fftw_rader_data *rader_data;
>               fftw_twiddle *tw;
>               struct fftw_plan_node_struct *recurse;
>          } rader;
>
>
>          struct {
>               int size;
>               fftw_real2hc_codelet *codelet;
>               const fftw_codelet_desc *codelet_desc;
>          } real2hc;
>
>
>          struct {
>               int size;
>               fftw_hc2real_codelet *codelet;
>               const fftw_codelet_desc *codelet_desc;
>          } hc2real;
>
>
>          struct {
>               int size;
>               fftw_direction dir;
>               fftw_hc2hc_codelet *codelet;
>               fftw_twiddle *tw;
>               struct fftw_plan_node_struct *recurse;
>               const fftw_codelet_desc *codelet_desc;
>          } hc2hc;
>
>
>          struct {
>               int size;
>               fftw_direction dir;
>               fftw_rgeneric_codelet *codelet;
>               fftw_twiddle *tw;
>               struct fftw_plan_node_struct *recurse;
>          } rgeneric;
>     } nodeu;
>
>     int refcnt;
>} fftw_plan_node;
>
>typedef enum {
>     FFTW_NORMAL_RECURSE = 0,
>     FFTW_VECTOR_RECURSE = 1
>} fftw_recurse_kind;
>
>struct fftw_plan_struct {
>     int n;
>     int refcnt;
>     fftw_direction dir;
>     int flags;
>     int wisdom_signature;
>     enum fftw_node_type wisdom_type;
>     struct fftw_plan_struct *next;
>     fftw_plan_node *root;
>     double cost;
>     fftw_recurse_kind recurse_kind;
>     int vector_size;
>};
>
>typedef struct fftw_plan_struct *fftw_plan;
># 300 "/software/mpbi/x86_64/install/fftw-2.1.5/include/fftw.h"
>extern fftw_plan fftw_create_plan_specific(int n, fftw_direction dir,
>                                           int flags,
>                                           fftw_complex *in, int istride,
>                                         fftw_complex *out, int ostride);
>
>extern fftw_plan fftw_create_plan(int n, fftw_direction dir, int flags);
>extern void fftw_print_plan(fftw_plan plan);
>extern void fftw_destroy_plan(fftw_plan plan);
>extern void fftw(fftw_plan plan, int howmany, fftw_complex *in, int istride,
>                 int idist, fftw_complex *out, int ostride, int odist);
>extern void fftw_one(fftw_plan plan, fftw_complex *in, fftw_complex *out);
>extern void fftw_die(const char *s);
>extern void *fftw_malloc(size_t n);
>extern void fftw_free(void *p);
>extern void fftw_check_memory_leaks(void);
>extern void fftw_print_max_memory_usage(void);
>
>typedef void *(*fftw_malloc_type_function) (size_t n);
>typedef void (*fftw_free_type_function) (void *p);
>typedef void (*fftw_die_type_function) (const char *errString);
>extern fftw_malloc_type_function fftw_malloc_hook;
>extern fftw_free_type_function fftw_free_hook;
>extern fftw_die_type_function fftw_die_hook;
>
>extern size_t fftw_sizeof_fftw_real(void);
>
>
>
>
>
>
>
>extern void fftw_forget_wisdom(void);
>extern void fftw_export_wisdom(void (*emitter) (char c, void *), void *data);
>extern fftw_status fftw_import_wisdom(int (*g) (void *), void *data);
>extern void fftw_export_wisdom_to_file(FILE *output_file);
>extern fftw_status fftw_import_wisdom_from_file(FILE *input_file);
>extern char *fftw_export_wisdom_to_string(void);
>extern fftw_status fftw_import_wisdom_from_string(const char *input_string);
>
>
>
>
>
>
>extern void fftw_fprint_plan(FILE *f, fftw_plan plan);
>
>
>
>
>typedef struct {
>     int is_in_place;
>
>     int rank;
>
>
>
>     int *n;
>
>
>
>     fftw_direction dir;
>
>     int *n_before;
>
>
>     int *n_after;
>
>     fftw_plan *plans;
>
>     int nbuffers, nwork;
>     fftw_complex *work;
>
>
>
>
>} fftwnd_data;
>
>typedef fftwnd_data *fftwnd_plan;
>
>
>extern fftwnd_plan fftw2d_create_plan(int nx, int ny, fftw_direction dir,
>                                      int flags);
>extern fftwnd_plan fftw3d_create_plan(int nx, int ny, int nz,
>                                      fftw_direction dir, int flags);
>extern fftwnd_plan fftwnd_create_plan(int rank, const int *n,
>                                      fftw_direction dir,
>                                      int flags);
>
>extern fftwnd_plan fftw2d_create_plan_specific(int nx, int ny,
>                                               fftw_direction dir,
>                                               int flags,
>                                           fftw_complex *in, int istride,
>                                         fftw_complex *out, int ostride);
>extern fftwnd_plan fftw3d_create_plan_specific(int nx, int ny, int nz,
>                                           fftw_direction dir, int flags,
>                                           fftw_complex *in, int istride,
>                                         fftw_complex *out, int ostride);
>extern fftwnd_plan fftwnd_create_plan_specific(int rank, const int *n,
>                                               fftw_direction dir,
>                                               int flags,
>                                           fftw_complex *in, int istride,
>                                         fftw_complex *out, int ostride);
>
>
>extern void fftwnd_destroy_plan(fftwnd_plan plan);
>
>
>extern void fftwnd_fprint_plan(FILE *f, fftwnd_plan p);
>extern void fftwnd_print_plan(fftwnd_plan p);
>
>
>
>extern void fftwnd(fftwnd_plan plan, int howmany,
>                   fftw_complex *in, int istride, int idist,
>                   fftw_complex *out, int ostride, int odist);
>extern void fftwnd_one(fftwnd_plan p, fftw_complex *in, fftw_complex *out);
># 26 "maxwell.h" 2
># 1 "/software/mpbi/x86_64/install/fftw-2.1.5/include/rfftw.h" 1
># 34 "/software/mpbi/x86_64/install/fftw-2.1.5/include/rfftw.h"
>typedef fftw_plan rfftw_plan;
>typedef fftwnd_plan rfftwnd_plan;
>
>
>
>
>extern void rfftw(rfftw_plan plan, int howmany, fftw_real *in, int istride,
>                  int idist, fftw_real *out, int ostride, int odist);
>extern void rfftw_one(rfftw_plan plan, fftw_real *in, fftw_real *out);
>
>extern rfftw_plan rfftw_create_plan_specific(int n, fftw_direction dir,
>                                            int flags,
>                                            fftw_real *in, int istride,
>                                            fftw_real *out, int ostride);
>
>extern rfftw_plan rfftw_create_plan(int n, fftw_direction dir, int flags);
>extern void rfftw_destroy_plan(rfftw_plan plan);
>
>extern void rfftw_fprint_plan(FILE *f, rfftw_plan p);
>extern void rfftw_print_plan(rfftw_plan p);
>
>extern void rfftw_executor_simple(int n, fftw_real *in,
>                                  fftw_real *out,
>                                  fftw_plan_node *p,
>                                  int istride,
>                                  int ostride,
>                                  fftw_recurse_kind recurse_kind);
>
>extern rfftwnd_plan rfftwnd_create_plan_specific(int rank, const int *n,
>                                                fftw_direction dir, int flags,
>                                                fftw_real *in, int istride,
>                                                fftw_real *out, int ostride);
>extern rfftwnd_plan rfftw2d_create_plan_specific(int nx, int ny,
>                                           fftw_direction dir, int flags,
>                                              fftw_real *in, int istride,
>                                            fftw_real *out, int ostride);
>extern rfftwnd_plan rfftw3d_create_plan_specific(int nx, int ny, int nz,
>                                           fftw_direction dir, int flags,
>                                              fftw_real *in, int istride,
>                                            fftw_real *out, int ostride);
>extern rfftwnd_plan rfftwnd_create_plan(int rank, const int *n,
>                                          fftw_direction dir, int flags);
>extern rfftwnd_plan rfftw2d_create_plan(int nx, int ny,
>                                          fftw_direction dir, int flags);
>extern rfftwnd_plan rfftw3d_create_plan(int nx, int ny, int nz,
>                                          fftw_direction dir, int flags);
>extern void rfftwnd_destroy_plan(rfftwnd_plan plan);
>extern void rfftwnd_fprint_plan(FILE *f, rfftwnd_plan plan);
>extern void rfftwnd_print_plan(rfftwnd_plan plan);
>extern void rfftwnd_real_to_complex(rfftwnd_plan p, int howmany,
>                                   fftw_real *in, int istride, int idist,
>                              fftw_complex *out, int ostride, int odist);
>extern void rfftwnd_complex_to_real(rfftwnd_plan p, int howmany,
>                                fftw_complex *in, int istride, int idist,
>                                 fftw_real *out, int ostride, int odist);
>extern void rfftwnd_one_real_to_complex(rfftwnd_plan p,
>                                        fftw_real *in, fftw_complex *out);
>extern void rfftwnd_one_complex_to_real(rfftwnd_plan p,
>                                        fftw_complex *in, fftw_real *out);
># 27 "maxwell.h" 2
># 55 "maxwell.h"
>typedef struct {
>     real kmag;
>     real mx, my, mz;
>     real nx, ny, nz;
>} k_data;
>
>
>
>
>typedef struct {
>
>
>
>
>     real m00, m01, m02,
>               m11, m12,
>                    m22;
>
>} symmetric_matrix;
># 91 "maxwell.h"
>typedef struct {
>     int nx, ny, nz;
>     int local_nx, local_ny;
>     int local_x_start, local_y_start;
>     int last_dim, last_dim_size, other_dims;
>
>     int num_bands;
>     int N, local_N, N_start, alloc_N;
>
>     int fft_output_size;
>
>     int max_fft_bands, num_fft_bands;
>
>     real current_k[3];
>     int parity;
># 116 "maxwell.h"
>     fftwnd_plan plan, iplan;
>
>
>
>
>
>
>     scalar *fft_data;
>
>     int zero_k;
>     k_data *k_plus_G;
>     real *k_plus_G_normsqr;
>
>     symmetric_matrix *eps_inv;
>     real eps_inv_mean;
>} maxwell_data;
>
>extern maxwell_data *create_maxwell_data(int nx, int ny, int nz,
>                                         int *local_N, int *N_start,
>                                         int *alloc_N,
>                                         int num_bands,
>                                         int num_fft_bands);
>extern void destroy_maxwell_data(maxwell_data *d);
>
>extern void maxwell_set_num_bands(maxwell_data *d, int num_bands);
>
>extern void update_maxwell_data_k(maxwell_data *d, real k[3],
>                                  real G1[3], real G2[3], real G3[3]);
>
>extern void set_maxwell_data_parity(maxwell_data *d, int parity);
>
>typedef void (*maxwell_dielectric_function) (symmetric_matrix *eps,
>                                             symmetric_matrix *eps_inv,
>                                             real r[3], void *epsilon_data);
>
>extern void set_maxwell_dielectric(maxwell_data *md,
>                                   const int mesh_size[3],
>                                   real R[3][3], real G[3][3],
>                                   maxwell_dielectric_function epsilon,
>                                   void *epsilon_data);
>
>extern void maxwell_sym_matrix_eigs(real eigs[3], const symmetric_matrix *V);
>extern void maxwell_sym_matrix_invert(symmetric_matrix *Vinv,
>                                      const symmetric_matrix *V);
>
>extern void maxwell_compute_fft(int dir, maxwell_data *d, scalar *array,
>                                int howmany, int stride, int dist);
>extern void maxwell_compute_d_from_H(maxwell_data *d, evectmatrix Xin,
>                                     scalar_complex *dfield,
>                                     int cur_band_start, int cur_num_bands);
>extern void maxwell_compute_h_from_H(maxwell_data *d, evectmatrix Hin,
>                                     scalar_complex *hfield,
>                                     int cur_band_start, int cur_num_bands);
>extern void maxwell_compute_e_from_d(maxwell_data *d,
>                                     scalar_complex *dfield,
>                                     int cur_num_bands);
>
>extern void maxwell_vectorfield_otherhalf(maxwell_data *d,
>                                          scalar_complex *field,
>                                          real phasex,real phasey,real phasez);
>extern void maxwell_scalarfield_otherhalf(maxwell_data *d, real *field);
>
>void assign_symmatrix_vector(scalar_complex *newv,
>                             const symmetric_matrix matrix,
>                             const scalar_complex *oldv);
>
>extern void maxwell_operator(evectmatrix Xin, evectmatrix Xout, void *data,
>                             int is_current_eigenvector, evectmatrix Work);
>extern void maxwell_simple_precondition(evectmatrix X,
>                                        void *data, real *eigenvals);
>extern void maxwell_preconditioner(evectmatrix Xin, evectmatrix Xout,
>                                   void *data,
>                                   evectmatrix Y, real *eigenvals,
>                                   sqmatrix YtY);
>extern void maxwell_preconditioner2(evectmatrix Xin, evectmatrix Xout,
>                                    void *data,
>                                    evectmatrix Y, real *eigenvals,
>                                    sqmatrix YtY);
>
>extern void maxwell_ucross_op(evectmatrix Xin, evectmatrix Xout,
>                              maxwell_data *d, const real u[3]);
>
>extern void maxwell_parity_constraint(evectmatrix X, void *data);
>extern void maxwell_zparity_constraint(evectmatrix X, void *data);
>extern void maxwell_yparity_constraint(evectmatrix X, void *data);
>
>extern int maxwell_zero_k_num_const_bands(evectmatrix X, maxwell_data *d);
>extern void maxwell_zero_k_set_const_bands(evectmatrix X, maxwell_data *d);
>extern void maxwell_zero_k_constraint(evectmatrix X, void *data);
>
>extern real *maxwell_zparity(evectmatrix X, maxwell_data *d);
>extern real *maxwell_yparity(evectmatrix X, maxwell_data *d);
>
>typedef struct {
>     maxwell_data *d;
>     real target_frequency;
>} maxwell_target_data;
>
>extern maxwell_target_data *create_maxwell_target_data(maxwell_data *d,
>                                                       real target_frequency);
>extern void destroy_maxwell_target_data(maxwell_target_data *d);
>extern void maxwell_target_operator1(evectmatrix Xin, evectmatrix Xout,
>                                     void *data,
>                                     int is_current_eigenvector,
>                                     evectmatrix Work);
>extern void maxwell_target_operator(evectmatrix Xin, evectmatrix Xout,
>                                    void *data, int is_current_eigenvector,
>                                    evectmatrix Work);
>extern void maxwell_target_preconditioner(evectmatrix Xin, evectmatrix Xout,
>                                          void *data,
>                                          evectmatrix Y, real *eigenvals,
>                                          sqmatrix YtY);
>extern void maxwell_target_preconditioner2(evectmatrix Xin, evectmatrix Xout,
>                                           void *data,
>                                           evectmatrix Y, real *eigenvals,
>                                           sqmatrix YtY);
>
>extern void spherical_quadrature_points(real *x, real *y, real *z,
>                                        real *weight, int num_sq_pts);
>
>extern int check_maxwell_dielectric(maxwell_data *d,
>                                    int negative_epsilon_okp);
># 28 "maxwell_eps.c" 2
># 39 "maxwell_eps.c"
>extern void zheev_(char *, char *, int *, scalar *, int *, real *,
>                 scalar *, int *, real *, int *);
>extern void dsyev_(char *, char *, int *, real *, int *, real *,
>                 real *, int *, int *);
>
>
>void maxwell_sym_matrix_eigs(real eigs[3], const symmetric_matrix *V)
>{
>     int n = 3, nw = 9, info;
># 58 "maxwell_eps.c"
>     real Vm[3][3], W[9];
>     Vm[0][0] = V->m00;
>     Vm[1][1] = V->m11;
>     Vm[2][2] = V->m22;
>     Vm[0][1] = Vm[1][0] = V->m01;
>     Vm[0][2] = Vm[2][0] = V->m02;
>     Vm[1][2] = Vm[2][1] = V->m12;
>     dsyev_("V", "U", &n, &Vm[0][0], &n, eigs, W, &nw, &info);
>
>     do { if (!(info >= 0)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "invalid argument in heev" "\n", 67); } } while (0);
>     do { if (!(info <= 0)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "failure to converge in heev" "\n", 68); } } while (0);
>}
>
>
>
>void maxwell_sym_matrix_invert(symmetric_matrix *Vinv,
>                               const symmetric_matrix *V)
>{
>     real m00 = V->m00, m11 = V->m11, m22 = V->m22;
># 123 "maxwell_eps.c"
>     real m01 = V->m01, m02 = V->m02, m12 = V->m12;
>
>     if (m01 == 0.0 && m02 == 0.0 && m12 == 0.0) {
>
>          Vinv->m00 = 1.0 / m00;
>          Vinv->m11 = 1.0 / m11;
>          Vinv->m22 = 1.0 / m22;
>          Vinv->m01 = Vinv->m02 = Vinv->m12 = 0.0;
>     }
>     else {
>          double detinv;
>
>
>          detinv = m00*m11*m22 - m02*m11*m02 + 2.0 * m01*m12*m02 -
>               m01*m01*m22 - m12*m12*m00;
>
>          do { if (!(detinv != 0.0)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "singular 3x3 matrix" "\n", 139); } } while (0);
>
>          detinv = 1.0/detinv;
>
>          Vinv->m00 = detinv * (m11*m22 - m12*m12);
>          Vinv->m11 = detinv * (m00*m22 - m02*m02);
>          Vinv->m22 = detinv * (m11*m00 - m01*m01);
>
>          Vinv->m02 = detinv * (m01*m12 - m11*m02);
>          Vinv->m01 = detinv * (m12*m02 - m01*m22);
>          Vinv->m12 = detinv * (m01*m02 - m00*m12);
>     }
>
>}
>
>
>static int sym_matrix_positive_definite(symmetric_matrix *V)
>{
>     real det2, det3;
>     real m00 = V->m00, m11 = V->m11, m22 = V->m22;
># 168 "maxwell_eps.c"
>     real m01 = V->m01, m02 = V->m02, m12 = V->m12;
>
>     det2 = m00*m11 - m01*m01;
>     det3 = det2*m22 - m02*m11*m02 + 2.0 * m01*m12*m02 - m12*m12*m00;
>
>
>     return (m00 > 0.0 && det2 > 0.0 && det3 > 0.0);
>}
>
>
>
>int check_maxwell_dielectric(maxwell_data *d,
>                             int negative_epsilon_okp)
>{
>     int i, require_2d;
>
>     require_2d = d->nz == 1 && (d->parity & ((1<<0) | (1<<1)));
>
>     for (i = 0; i < d->fft_output_size; ++i) {
>          if (!negative_epsilon_okp &&
>              !sym_matrix_positive_definite(d->eps_inv + i))
>               return 1;
>          if (require_2d) {
>
>
>
>
>
>
>
>               if (d->eps_inv[i].m02 != 0.0 || d->eps_inv[i].m12 != 0.0)
>                    return 2;
>
>          }
>     }
>     return 0;
>}
># 214 "maxwell_eps.c"
>void spherical_quadrature_points(real *x, real *y, real *z,
>                                 real *weight, int num_sq_pts)
>{
>     int i,j,k,l, n = 0;
>     real x0, y0, z0, w;
>
>     if (num_sq_pts == 50) {
># 237 "maxwell_eps.c"
>          x0 = 1; y0 = z0 = 0;
>          w = 9216 / 725760.0;
>          for (i = 0; i < 2; ++i) {
>               x0 = -x0;
>               for (j = 0; j < 3; ++j) {
>                    {real SHIFT3_dummy = z0; z0 = y0; y0 = x0; x0 = SHIFT3_dummy;};
>                    x[n] = x0; y[n] = y0; z[n] = z0; weight[n++] = w;
>               }
>          }
>
>          x0 = y0 = sqrt(0.5); z0 = 0;
>          w = 16384 / 725760.0;
>          for (i = 0; i < 2; ++i) {
>               x0 = -x0;
>               for (j = 0; j < 2; ++j) {
>                    y0 = -y0;
>                    for (k = 0; k < 3; ++k) {
>                         {real SHIFT3_dummy = z0; z0 = y0; y0 = x0; x0 = SHIFT3_dummy;};
>                         x[n] = x0; y[n] = y0; z[n] = z0; weight[n++] = w;
>                    }
>               }
>          }
>
>          x0 = y0 = z0 = sqrt(1.0 / 3.0);
>          w = 15309 / 725760.0;
>          for (i = 0; i < 2; ++i) {
>               x0 = -x0;
>               for (j = 0; j < 2; ++j) {
>                    y0 = -y0;
>                    for (k = 0; k < 2; ++k) {
>                         z0 = -z0;
>                         x[n] = x0; y[n] = y0; z[n] = z0; weight[n++] = w;
>                    }
>               }
>          }
>
>          x0 = y0 = sqrt(1.0 / 11.0); z0 = 3 * x0;
>          w = 14641 / 725760.0;
>          for (i = 0; i < 2; ++i) {
>               x0 = -x0;
>               for (j = 0; j < 2; ++j) {
>                    y0 = -y0;
>                    for (k = 0; k < 2; ++k) {
>                         z0 = -z0;
>                         for (l = 0; l < 3; ++l) {
>                              {real SHIFT3_dummy = z0; z0 = y0; y0 = x0; x0 = SHIFT3_dummy;};
>                              x[n] = x0; y[n] = y0; z[n] = z0; weight[n++] = w;
>                         }
>                    }
>               }
>          }
>     }
>     else if (num_sq_pts == 72 || num_sq_pts == 12) {
>
>
>
>
>
>
>
>          x0 = sqrt(0.5 - sqrt(0.05));
>          y0 = sqrt(0.5 + sqrt(0.05));
>          z0 = 0;
>          if (num_sq_pts == 72)
>               w = 125 / 10080.0;
>          else
>               w = 1 / 12.0;
>          for (i = 0; i < 2; ++i) {
>               x0 = -x0;
>               for (j = 0; j < 2; ++j) {
>                    y0 = -y0;
>                    for (k = 0; k < 3; ++k) {
>                         {real SHIFT3_dummy = z0; z0 = y0; y0 = x0; x0 = SHIFT3_dummy;};
>                         x[n] = x0; y[n] = y0; z[n] = z0; weight[n++] = w;
>                    }
>               }
>          }
>
>          if (num_sq_pts == 72) {
>
>
>               real coords[3][5] = {
>                    { -0.151108275, 0.315838353, 0.346307112, -0.101808787, -0.409228403 },
>                    { 0.155240600, 0.257049387, 0.666277790, 0.817386065, 0.501547712 },
>                    { 0.976251323, 0.913330032, 0.660412970, 0.567022920, 0.762221757 }
>               };
>
>               w = 143 / 10080.0;
>               for (l = 0; l < 5; ++l) {
>                    x0 = coords[0][l]; y0 = coords[1][l]; z0 = coords[2][l];
>                    for (i = 0; i < 3; ++i) {
>                         double dummy = x0;
>                         x0 = z0;
>                         z0 = -y0;
>                         y0 = -dummy;
>                         for (j = 0; j < 3; ++j) {
>                              {real SHIFT3_dummy = z0; z0 = y0; y0 = x0; x0 = SHIFT3_dummy;};
>                              x[n] = x0; y[n] = y0; z[n] = z0; weight[n++] = w;
>                         }
>                         y0 = -y0;
>                         z0 = -z0;
>                         x[n] = x0; y[n] = y0; z[n] = z0; weight[n++] = w;
>                    }
>               }
>          }
>     }
>     else
>          mpi_die("spherical_quadrature_points: passed unknown # points!");
>
>     do { if (!(n == num_sq_pts)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "bug in spherical_quadrature_points: wrong number of points!" "\n", 347); } } while (0);
>
>}
># 374 "maxwell_eps.c"
>static void get_mesh(int nx, int ny, int nz, const int mesh_size[3],
>                     real R[3][3], real G[3][3],
>                     real mesh_center[3], int *mesh_prod,
>                     real moment_mesh[50][3],
>                     real moment_mesh_weights[50],
>                     int *size_moment_mesh)
>{
>     int i,j;
>     real min_diam = 1e20;
>     real mesh_total[3] = { 0, 0, 0 };
>     int rank = nz > 1 ? 3 : (ny > 1 ? 2 : 1);
>     real weight_sum = 0.0;
>
>     *mesh_prod = 1;
>     for (i = 0; i < 3; ++i) {
>          int ms = ((mesh_size[i]) > (1) ? (mesh_size[i]) : (1));
>          mesh_center[i] = (ms - 1) * 0.5;
>          *mesh_prod *= ms;
>     }
>
>
>
>
>
>     if (rank == 1) {
>
>
>          *size_moment_mesh = 2;
>          moment_mesh[0][0] = -0.5;
>          moment_mesh[0][1] = 0.0;
>          moment_mesh[0][2] = 0.0;
>          moment_mesh[1][0] = 0.5;
>          moment_mesh[1][1] = 0.0;
>          moment_mesh[1][2] = 0.0;
>          moment_mesh_weights[0] = moment_mesh_weights[1] = 0.5;
>     }
>     else if (rank == 2) {
>
>          *size_moment_mesh = 12;
>          for (i = 0; i < *size_moment_mesh; ++i) {
>               moment_mesh[i][0] =
>                    cos(2*i * 3.141592653589793238462643383279502884197 / *size_moment_mesh) * 0.5;
>               moment_mesh[i][1] =
>                    sin(2*i * 3.141592653589793238462643383279502884197 / *size_moment_mesh) * 0.5;
>               moment_mesh[i][2] = 0.0;
>               moment_mesh_weights[i] = 1.0 / *size_moment_mesh;
>          }
>     }
>     else {
>          real x[50], y[50], z[50];
>
>          *size_moment_mesh = 50;
>          spherical_quadrature_points(x,y,z, moment_mesh_weights, 50);
>          for (i = 0; i < *size_moment_mesh; ++i) {
>               moment_mesh[i][0] = x[i] * 0.5;
>               moment_mesh[i][1] = y[i] * 0.5;
>               moment_mesh[i][2] = z[i] * 0.5;
>          }
>     }
>
>     do { if (!(*size_moment_mesh <= 50)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "yikes, moment mesh too big" "\n", 434); } } while (0);
>
>     for (i = 0; i < *size_moment_mesh; ++i)
>          weight_sum += moment_mesh_weights[i];
>     do { if (!(fabs(weight_sum - 1.0) < 1.0e-6)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "bug, incorrect moment weights" "\n", 438); } } while (0);
>
>
>
>
>
>
>
>     for (i = 0; i < rank; ++i) {
>          real ri = R[i][0] * R[i][0] + R[i][1] * R[i][1] + R[i][2] * R[i][2];
>          ri = sqrt(ri) / (i == 0 ? nx : (i == 1 ? ny : nz));
>          min_diam = ((min_diam) < (ri) ? (min_diam) : (ri));
>     }
>
>
>     for (i = 0; i < *size_moment_mesh; ++i) {
>          real len = 0;
>          for (j = 0; j < 3; ++j) {
>               moment_mesh[i][j] *= min_diam;
>               len += moment_mesh[i][j] * moment_mesh[i][j];
>               mesh_total[j] += moment_mesh[i][j];
>          }
>          do { if (!(fabs(len - min_diam*min_diam*(0.5*0.5)) < 1.0e-6)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "bug in get_mesh: moment_mesh vector is wrong length" "\n", 462); } } while (0);
>
>
>     }
>     do { if (!(fabs(mesh_total[0]) + fabs(mesh_total[1]) + fabs(mesh_total[2]) < 1.0e-6)) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "bug in get_mesh: moment_mesh does not average to zero" "\n", 465); } } while (0);
>
>
>
>
>     for (i = 0; i < *size_moment_mesh; ++i) {
>          real v[3];
>          for (j = 0; j < 3; ++j)
>               v[j] = moment_mesh[i][j];
>          for (j = 0; j < 3; ++j)
>               moment_mesh[i][j] = G[j][0]*v[0] + G[j][1]*v[1] + G[j][2]*v[2];
>     }
>}
># 502 "maxwell_eps.c"
>void set_maxwell_dielectric(maxwell_data *md,
>                            const int mesh_size[3],
>                            real R[3][3], real G[3][3],
>                            maxwell_dielectric_function epsilon,
>                            void *epsilon_data)
>{
>     real s1, s2, s3, m1, m2, m3;
>     real mesh_center[3];
>     real moment_mesh[50][3];
>     real moment_mesh_weights[50];
>     real eps_inv_total = 0.0;
>     int i, j, k;
>     int mesh_prod;
>     real mesh_prod_inv;
>     int size_moment_mesh = 0;
>     int n1, n2, n3;
>
>
>
>
>
>
>
>     n1 = md->nx; n2 = md->ny; n3 = md->nz;
>
>     get_mesh(n1, n2, n3, mesh_size, R, G,
>              mesh_center, &mesh_prod, moment_mesh, moment_mesh_weights,
>              &size_moment_mesh);
>     mesh_prod_inv = 1.0 / mesh_prod;
># 540 "maxwell_eps.c"
>     s1 = 1.0 / n1;
>     s2 = 1.0 / n2;
>     s3 = 1.0 / n3;
>     m1 = s1 / ((1) > (mesh_size[0] - 1) ? (1) : (mesh_size[0] - 1));
>     m2 = s2 / ((1) > (mesh_size[1] - 1) ? (1) : (mesh_size[1] - 1));
>     m3 = s3 / ((1) > (mesh_size[2] - 1) ? (1) : (mesh_size[2] - 1));
># 558 "maxwell_eps.c"
>     for (i = 0; i < n1; ++i)
>          for (j = 0; j < n2; ++j)
>               for (k = 0; k < n3; ++k)
>     {
>
>
>
>          int eps_index = ((i * n2 + j) * n3 + k);
># 630 "maxwell_eps.c"
>     {
>          int mi, mj, mk;
>
>
>
>
>
>          symmetric_matrix eps_mean = {0,0,0,0,0,0},
>                           eps_inv_mean = {0,0,0,0,0,0}, eps_mean_inv;
>
>          real norm_len;
>          real norm0, norm1, norm2;
>          short means_different_p, diag_eps_p;
>
>          for (mi = 0; mi < mesh_size[0]; ++mi)
>               for (mj = 0; mj < mesh_size[1]; ++mj)
>                    for (mk = 0; mk < mesh_size[2]; ++mk) {
>                         real r[3];
>                         symmetric_matrix eps, eps_inv;
>                         r[0] = i * s1 + (mi - mesh_center[0]) * m1;
>                         r[1] = j * s2 + (mj - mesh_center[1]) * m2;
>                         r[2] = k * s3 + (mk - mesh_center[2]) * m3;
>                         epsilon(&eps, &eps_inv, r, epsilon_data);
>                         eps_mean.m00 += eps.m00;
>                         eps_mean.m11 += eps.m11;
>                         eps_mean.m22 += eps.m22;
>                         eps_inv_mean.m00 += eps_inv.m00;
>                         eps_inv_mean.m11 += eps_inv.m11;
>                         eps_inv_mean.m22 += eps_inv.m22;
># 667 "maxwell_eps.c"
>                         eps_mean.m01 += eps.m01;
>                         eps_mean.m02 += eps.m02;
>                         eps_mean.m12 += eps.m12;
>                         eps_inv_mean.m01 += eps_inv.m01;
>                         eps_inv_mean.m02 += eps_inv.m02;
>                         eps_inv_mean.m12 += eps_inv.m12;
>
>                    }
>
>          diag_eps_p = ((eps_mean).m01 == 0.0 && (eps_mean).m02 == 0.0 && (eps_mean).m12 == 0.0);
>          if (diag_eps_p) {
>               eps_mean_inv.m00 = mesh_prod / eps_mean.m00;
>               eps_mean_inv.m11 = mesh_prod / eps_mean.m11;
>               eps_mean_inv.m22 = mesh_prod / eps_mean.m22;
>
>
>
>
>
>               eps_mean_inv.m01 = eps_mean_inv.m02 = eps_mean_inv.m12 = 0.0;
>
>               eps_inv_mean.m00 *= mesh_prod_inv;
>               eps_inv_mean.m11 *= mesh_prod_inv;
>               eps_inv_mean.m22 *= mesh_prod_inv;
>
>               means_different_p =
>                    fabs(eps_mean_inv.m00 - eps_inv_mean.m00) > 1.0e-6 ||
>                    fabs(eps_mean_inv.m11 - eps_inv_mean.m11) > 1.0e-6 ||
>                    fabs(eps_mean_inv.m22 - eps_inv_mean.m22) > 1.0e-6;
>          }
>          else {
>               eps_inv_mean.m00 *= mesh_prod_inv;
>               eps_inv_mean.m11 *= mesh_prod_inv;
>               eps_inv_mean.m22 *= mesh_prod_inv;
>               eps_mean.m00 *= mesh_prod_inv;
>               eps_mean.m11 *= mesh_prod_inv;
>               eps_mean.m22 *= mesh_prod_inv;
># 718 "maxwell_eps.c"
>               eps_mean.m01 *= mesh_prod_inv;
>               eps_mean.m02 *= mesh_prod_inv;
>               eps_mean.m12 *= mesh_prod_inv;
>               eps_inv_mean.m01 *= mesh_prod_inv;
>               eps_inv_mean.m02 *= mesh_prod_inv;
>               eps_inv_mean.m12 *= mesh_prod_inv;
>
>               maxwell_sym_matrix_invert(&eps_mean_inv, &eps_mean);
>
>               means_different_p =
>                    fabs(eps_mean_inv.m00 - eps_inv_mean.m00) > 1.0e-6 ||
>                    fabs(eps_mean_inv.m11 - eps_inv_mean.m11) > 1.0e-6 ||
>                    fabs(eps_mean_inv.m22 - eps_inv_mean.m22) > 1.0e-6;
># 740 "maxwell_eps.c"
>               means_different_p = means_different_p ||
>                    fabs(eps_mean_inv.m01 - eps_inv_mean.m01) > 1.0e-6 ||
>                    fabs(eps_mean_inv.m02 - eps_inv_mean.m02) > 1.0e-6 ||
>                    fabs(eps_mean_inv.m12 - eps_inv_mean.m12) > 1.0e-6;
>
>          }
>
>
>
>
>          if (means_different_p) {
>               real moment0 = 0, moment1 = 0, moment2 = 0;
>
>               for (mi = 0; mi < size_moment_mesh; ++mi) {
>                    real r[3], eps_trace;
>                    symmetric_matrix eps, eps_inv;
>                    r[0] = i * s1 + moment_mesh[mi][0];
>                    r[1] = j * s2 + moment_mesh[mi][1];
>                    r[2] = k * s3 + moment_mesh[mi][2];
>                    epsilon(&eps, &eps_inv, r, epsilon_data);
>                    eps_trace = eps.m00 + eps.m11 + eps.m22;
>                    eps_trace *= moment_mesh_weights[mi];
>                    moment0 += eps_trace * moment_mesh[mi][0];
>                    moment1 += eps_trace * moment_mesh[mi][1];
>                    moment2 += eps_trace * moment_mesh[mi][2];
>               }
>
>               norm0 = R[0][0]*moment0 + R[1][0]*moment1 + R[2][0]*moment2;
>               norm1 = R[0][1]*moment0 + R[1][1]*moment1 + R[2][1]*moment2;
>               norm2 = R[0][2]*moment0 + R[1][2]*moment1 + R[2][2]*moment2;
>
>               norm_len = sqrt(norm0*norm0 + norm1*norm1 + norm2*norm2);
>          }
>
>          if (means_different_p && norm_len > 1.0e-6) {
>               real x0, x1, x2;
>
>               norm_len = 1.0/norm_len;
>               norm0 *= norm_len;
>               norm1 *= norm_len;
>               norm2 *= norm_len;
># 801 "maxwell_eps.c"
>               x0 = (eps_inv_mean.m00 - eps_mean_inv.m00) * norm0;
>               x1 = (eps_inv_mean.m11 - eps_mean_inv.m11) * norm1;
>               x2 = (eps_inv_mean.m22 - eps_mean_inv.m22) * norm2;
>               if (diag_eps_p) {
># 813 "maxwell_eps.c"
>                    md->eps_inv[eps_index].m01 = 0.5*(x0*norm1 + x1*norm0);
>                    md->eps_inv[eps_index].m02 = 0.5*(x0*norm2 + x2*norm0);
>                    md->eps_inv[eps_index].m12 = 0.5*(x1*norm2 + x2*norm1);
>
>               }
>               else {
># 847 "maxwell_eps.c"
>                    x0 += ((eps_inv_mean.m01 - eps_mean_inv.m01) * norm1 +
>                           (eps_inv_mean.m02 - eps_mean_inv.m02) * norm2);
>                    x1 += ((eps_inv_mean.m01 - eps_mean_inv.m01) * norm0 +
>                           (eps_inv_mean.m12 - eps_mean_inv.m12) * norm2);
>                    x2 += ((eps_inv_mean.m02 - eps_mean_inv.m02) * norm0 +
>                           (eps_inv_mean.m12 - eps_mean_inv.m12) * norm1);
>
>                    md->eps_inv[eps_index].m01 = (0.5*(x0*norm1 + x1*norm0)
>                                                  + eps_mean_inv.m01);
>                    md->eps_inv[eps_index].m02 = (0.5*(x0*norm2 + x2*norm0)
>                                                  + eps_mean_inv.m02);
>                    md->eps_inv[eps_index].m12 = (0.5*(x1*norm2 + x2*norm1)
>                                                  + eps_mean_inv.m12);
>
>               }
>               md->eps_inv[eps_index].m00 = x0*norm0 + eps_mean_inv.m00;
>               md->eps_inv[eps_index].m11 = x1*norm1 + eps_mean_inv.m11;
>               md->eps_inv[eps_index].m22 = x2*norm2 + eps_mean_inv.m22;
>          }
>          else {
>               md->eps_inv[eps_index] = eps_mean_inv;
>          }
>
>          eps_inv_total += (md->eps_inv[eps_index].m00 +
>                            md->eps_inv[eps_index].m11 +
>                            md->eps_inv[eps_index].m22);
>     }}
>
>     { real bbbb = *(&eps_inv_total); { do { if (!((&bbbb) != ((&eps_inv_total)))) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "MPI_Allreduce doesn't work for sendbuf == recvbuf" "\n", 876); } } while (0); memcpy(((&eps_inv_total)), (&bbbb), (1) * sizeof(real)); }; };
>
>     n1 = md->fft_output_size;
>     { int bbbb = *(&n1); { do { if (!((&bbbb) != ((&n1)))) { ; mpi_die("CHECK failure on line %d of " "maxwell_eps.c" ": " "MPI_Allreduce doesn't work for sendbuf == recvbuf" "\n", 878); } } while (0); memcpy(((&n1)), (&bbbb), (1) * sizeof(int)); }; };
>     md->eps_inv_mean = eps_inv_total / (3 * n1);
>}

Actions: View

Attachments on bug 249250: 159770