Fork me on GitHub

source: svn/trunk/Utilities/stdhep/mcf_hepup_xdr.c@ 526

Last change on this file since 526 was 3, checked in by Xavier Rouby, 16 years ago

first commit
File size: 8.0 KB
Line 
1/*******************************************************************************
2* *
3* mcf_hepup_xdr.c -- XDR Utility routines for the Block Stdhep filters *
4* *
5* Copyright (c) 1994 Universities Research Association, Inc. *
6* All rights reserved. *
7* *
8* This material resulted from work developed under a Government Contract and *
9* is subject to the following license: The Government retains a paid-up, *
10* nonexclusive, irrevocable worldwide license to reproduce, prepare derivative *
11* works, perform publicly and display publicly by or for the Government, *
12* including the right to distribute to other Government contractors. Neither *
13* the United States nor the United States Department of Energy, nor any of *
14* their employees, makes any warranty, express or implied, or assumes any *
15* legal liability or responsibility for the accuracy, completeness, or *
16* usefulness of any information, apparatus, product, or process disclosed, or *
17* represents that its use would not infringe privately owned rights. *
18* *
19* *
20* Written by Paul Lebrun, Lynn Garren *
21* *
22* *
23*******************************************************************************/
24#include <stdio.h>
25#include <string.h>
26#include <sys/param.h>
27#include <rpc/types.h>
28#include <sys/types.h>
29#include <rpc/xdr.h>
30#include <limits.h>
31#ifdef SUNOS
32#include <floatingpoint.h>
33#else /* SUNOS */
34#include <float.h>
35#endif /* SUNOS */
36#include <stdlib.h>
37#include <time.h>
38#include "Utilities/mcfio/mcfio_Dict.h"
39#include "Utilities/mcfio/mcf_xdr.h"
40#include "Utilities/stdhep/hepeup.h"
41#include "Utilities/stdhep/heprup.h"
42#include "Utilities/stdhep/stdver.h"
43#ifndef FALSE
44#define FALSE 0
45#endif
46#ifndef TRUE
47#define TRUE 1
48#endif
49
50extern struct stdver stdver_;
51extern struct hepeup hepeup_;
52extern struct heprup heprup_;
53
54bool_t xdr_hepeup_(XDR *xdrs, int *blockid,
55 int *ntot, char** version)
56
57{
58/* Translate the HEPEUP COMMON block from the STDHEP package to/from
59 an XDR stream. Note that we do not allocate memory, because we fill
60 directly the COMMON. Also, mcfio will allocate the space for the
61 string version. */
62
63 unsigned int nn, nn2, nn5;
64 int *idat;
65 double *dat;
66
67 if ((xdrs->x_op == XDR_ENCODE) || (xdrs->x_op == XDR_MCFIOCODE)) {
68 if (*blockid != MCFIO_HEPEUP) {
69 fprintf (stderr, "mcf_hepup_xdr: Inconsistent Blockid %d \n ",
70 (*blockid));
71 return FALSE;
72 }
73
74 nn = (unsigned int) hepeup_.nup; /* Number of elements in idup, istup, vtimup, spinup */
75 nn2 = 2*(unsigned int) hepeup_.nup; /* Number of elements in mothup, icolup */
76 nn5 = 5*(unsigned int) hepeup_.nup; /* Number of elements in pup */
77
78 /* Total length in bytes include blockid, ntot, version, as well
79 as the common block entries */
80 *ntot = sizeof(int)*(5 + 2*nn + 2*nn2) + sizeof(double)*(4 + 2*nn + nn5);
81
82 if (xdrs->x_op == XDR_MCFIOCODE) return TRUE;
83 strncpy(version[0],stdver_.stdhep_ver, 4);
84 }
85
86 if ( (xdr_int(xdrs, blockid) &&
87 xdr_int(xdrs, ntot) &&
88 xdr_string(xdrs, version, MCF_XDR_VERSION_LENGTH) &&
89 xdr_int(xdrs, &(hepeup_.nup)) &&
90 xdr_int(xdrs, &(hepeup_.idprup))) == FALSE) return FALSE;
91
92 if ((xdrs->x_op == XDR_DECODE) && ( *blockid != MCFIO_HEPEUP) ) {
93 fprintf (stderr, "mcf_hepup_xdr: Inconsistent Blockid %d \n ",
94 (*blockid));
95 return FALSE;
96 }
97
98 if ( xdr_double(xdrs, &(hepeup_.xwgtup) ) == FALSE) return FALSE;
99 if ( xdr_double(xdrs, &(hepeup_.scalup) ) == FALSE) return FALSE;
100 if ( xdr_double(xdrs, &(hepeup_.aqedup) ) == FALSE) return FALSE;
101 if ( xdr_double(xdrs, &(hepeup_.aqcdup) ) == FALSE) return FALSE;
102 idat = hepeup_.idup;
103 if ( xdr_array(xdrs, (char **) &idat,
104 &nn, MAXNUP, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
105 idat = hepeup_.istup;
106 if ( xdr_array(xdrs, (char **) &idat,
107 &nn, MAXNUP, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
108 idat = (int *) hepeup_.mothup;
109 if ( xdr_array(xdrs, (char **) &idat,
110 &nn2, 2*MAXNUP, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
111 idat = (int *) hepeup_.icolup;
112 if ( xdr_array(xdrs, (char **) &idat,
113 &nn2, 2*MAXNUP, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
114 dat = (double *) hepeup_.pup;
115 if ( xdr_array(xdrs, (char **) &dat,
116 &nn5, 5*MAXNUP, sizeof(double), (void *) xdr_double) == FALSE) return FALSE;
117 dat = (double *) hepeup_.vtimup;
118 if ( xdr_array(xdrs, (char **) &dat,
119 &nn, MAXNUP, sizeof(double), (void *) xdr_double) == FALSE) return FALSE;
120 dat = (double *) hepeup_.spinup;
121 if ( xdr_array(xdrs, (char **) &dat,
122 &nn, MAXNUP, sizeof(double), (void *) xdr_double) == FALSE) return FALSE;
123 return TRUE;
124}
125
126bool_t xdr_heprup_(XDR *xdrs, int *blockid,
127 int *ntot, char** version)
128
129{
130/* Translate the HEPRUP COMMON block from the STDHEP package to/from
131 an XDR stream. Note that we do not allocate memory, because we fill
132 directly the COMMON. Also, mcfio will allocate the space for the
133 string version. */
134
135 unsigned int nn, n2;
136 int *idat;
137 double *dat;
138
139 if ((xdrs->x_op == XDR_ENCODE) || (xdrs->x_op == XDR_MCFIOCODE)) {
140 if (*blockid != MCFIO_HEPRUP) {
141 fprintf (stderr, "mcf_hepup_xdr: Inconsistent Blockid %d \n ",
142 (*blockid));
143 return FALSE;
144 }
145
146 nn = (unsigned int) heprup_.nprup; /* Number of elements in xsecup, xerrup, xmaxup, lprup */
147 n2 = (unsigned int) 2; /* Number of elements in idbmup, ebmup, pdfgup, pdfsup */
148
149 /* Total length in bytes include blockid, ntot, version, as well
150 as the common block entries */
151 *ntot = sizeof(int)*(5 + 3*n2 + nn) + sizeof(double)*(n2 + 3*nn);
152
153 if (xdrs->x_op == XDR_MCFIOCODE) return TRUE;
154 strncpy(version[0],stdver_.stdhep_ver, 4);
155 }
156
157 if ( (xdr_int(xdrs, blockid) &&
158 xdr_int(xdrs, ntot) &&
159 xdr_string(xdrs, version, MCF_XDR_VERSION_LENGTH) &&
160 xdr_int(xdrs, &(heprup_.idwtup)) &&
161 xdr_int(xdrs, &(heprup_.nprup))) == FALSE) return FALSE;
162
163 if ((xdrs->x_op == XDR_DECODE) && ( *blockid != MCFIO_HEPRUP) ) {
164 fprintf (stderr, "mcf_hepup_xdr: Inconsistent Blockid %d \n ",
165 (*blockid));
166 return FALSE;
167 }
168 idat = heprup_.idbmup;
169 if ( xdr_array(xdrs, (char **) &idat,
170 &n2, 2, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
171 dat = (double *) heprup_.ebmup;
172 if ( xdr_array(xdrs, (char **) &dat,
173 &n2, 2, sizeof(double), (void *) xdr_double) == FALSE) return FALSE;
174 idat = heprup_.pdfgup;
175 if ( xdr_array(xdrs, (char **) &idat,
176 &n2, 2, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
177 idat = heprup_.pdfsup;
178 if ( xdr_array(xdrs, (char **) &idat,
179 &n2, 2, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
180 dat = (double *) heprup_.xsecup;
181 if ( xdr_array(xdrs, (char **) &dat,
182 &nn, MAXPUP, sizeof(double), (void *) xdr_double) == FALSE) return FALSE;
183 dat = (double *) heprup_.xerrup;
184 if ( xdr_array(xdrs, (char **) &dat,
185 &nn, MAXPUP, sizeof(double), (void *) xdr_double) == FALSE) return FALSE;
186 dat = (double *) heprup_.xmaxup;
187 if ( xdr_array(xdrs, (char **) &dat,
188 &nn, MAXPUP, sizeof(double), (void *) xdr_double) == FALSE) return FALSE;
189 idat = heprup_.lprup;
190 if ( xdr_array(xdrs, (char **) &idat,
191 &nn, MAXPUP, sizeof(int), (void *) xdr_int) == FALSE) return FALSE;
192
193 return TRUE;
194}
195
196
Note: See TracBrowser for help on using the repository browser.