From: Lee, Kerry T. \(JSC-SF\)[LMIT] (kerry.t.lee@nasa.gov)
Date: Mon Dec 11 2006 - 16:00:07 CET
This is a multi-part message in MIME format.
------_=_NextPart_001_01C71D35.0BE304A4
Content-Type: text/plain;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
Hi Alfredo,
I have determined that the algorithm is doing things right and even =
though the what1 of the input file does not follow the manual's =
instructions it turns out to give the right results anyway. The =
difference is that instead of making "ia" a positive number only, if one =
makes it negative the algorithm produces the right result anyway. The =
following code does the work right.
AAUSYL =3D 0.01D+00 * DBLE (ITUSYL(NY))
IAUSYL (NY) =3D NINT (AAUSYL)
IEUSYL (NY) =3D ITUSYL (NY) - IAUSYL (NY) * 100
If you put in ITUSYL =3D -123, then IAUSYL =3D -1 and IEUSYL =3D -23, =
but for ITUSYL =3D 77 you get IAUSYL =3D 1 and IEUSYL =3D -23, but I =
assume since IA represents a single bin the sign is ignored and when =
what1 is packed back up before being written to the USRYIELD output the =
algorithm ignores the sign again by doing
ITUSYL(NY) =3D IEUSYL(NY) + ABS(IAUSYL(NY))*100
Although this is a little confusing on how I was getting the right =
results while not following the explicit instructions of the manual, and =
that the written what1 was different than the supplied what1, the fault =
was mine. Even so, it might be better to not ignore the sign of IA when =
writing to the USRYIELD output so that it will agree with what the user =
input.
Thanks
Kerry
-----Original Message-----
From: Alfredo Ferrari [mailto:alfredo.ferrari@cern.ch]
Sent: Mon 12/11/2006 7:35 AM
To: Lee, Kerry T. (JSC-SF)[LMIT]
Cc: fluka-discuss@fluka.org
Subject: RE: USRYIELD binary output
=20
Hi Kerry
the code tries to unpack from what(1) the infos about which two physical =
quantities you want to deal with. Probably a nonsense number is unpacked
into something strange (I should verify).
Ciao
Alfredo
On Fri, 8 Dec 2006, Lee, Kerry T. (JSC-SF)[LMIT] wrote:
> This is a multi-part message in MIME format.
>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0
> Content-Type: text/plain;
> charset=3D"iso-8859-1"
> Content-Transfer-Encoding: quoted-printable
>
> I seem to have found that I was in error when explaining the problem. =
I =3D
> missed the fact that what12 is stored as the next variable after =
what1, =3D
> but the value stored where what1 should be is 77 and the value from =
the =3D
> input file is -123. How is this not being transferred correctly to =
the =3D
> output file?
>
> Thanks
> Kerry
>
>
> -----Original Message-----
> From: Lee, Kerry T. (JSC-SF)[LMIT]
> Sent: Thu 12/7/2006 10:27 AM
> To: fluka-discuss@fluka.org
> Subject: USRYIELD binary output
> =3D20
>
> Dear FLUKA authors,
>
> I am reading a USRYIELD output that was written according to the =3D
> following USRYIELD command
>
> USRYIELD -123.00000 201.00000-50.000000 4.0000000 5.0000000 =3D
> 1.0000000LET+0.1 =3D20
> USRYIELD 1000.0000 0.0000000 100.00000 1.5000000 0.0000000 2603.0000 =
& =3D
> =3D20
> =3D20
> I am trying to parse the file using a C++ routine, so it is imperitive =
I =3D
> understand the details of what is written where. I followed the code =
in =3D
> the manual on reading the USRYIELD output =3D
> (http://www.fluka.org/manual/sect/s090/text.html) and see that the =3D
> variable ITUSYL(NY) is supposed to correspond to the What1. When I =
read =3D
> this file the variable reads a value of 2603, which is my What12 and =
not =3D
> What1. I have taken into account the extra bytes that fortran writes =
=3D
> and all the other values are read in correctly, so is this a bug or am =
I =3D
> missing something?
>
> Thanks
> Kerry
>
>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0
> Content-Type: text/html;
> charset=3D"iso-8859-1"
> Content-Transfer-Encoding: quoted-printable
>
> <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
> <HTML>
> <HEAD>
> <META HTTP-EQUIV=3D3D"Content-Type" CONTENT=3D3D"text/html; =3D
> charset=3D3Diso-8859-1">
> <META NAME=3D3D"Generator" CONTENT=3D3D"MS Exchange Server version =3D
> 6.5.7651.34">
> <TITLE>RE: USRYIELD binary output</TITLE>
> </HEAD>
> <BODY>
> <!-- Converted from text/plain format -->
>
> <P><FONT SIZE=3D3D2>I seem to have found that I was in error when =3D
> explaining the problem. I missed the fact that what12 is stored =
as =3D
> the next variable after what1, but the value stored where what1 should =
=3D
> be is 77 and the value from the input file is -123. How is this =
=3D
> not being transferred correctly to the output file?<BR>
> <BR>
> Thanks<BR>
> Kerry<BR>
> <BR>
> <BR>
> -----Original Message-----<BR>
> From: Lee, Kerry T. (JSC-SF)[LMIT]<BR>
> Sent: Thu 12/7/2006 10:27 AM<BR>
> To: fluka-discuss@fluka.org<BR>
> Subject: USRYIELD binary output<BR>
> <BR>
> <BR>
> Dear FLUKA authors,<BR>
> <BR>
> I am reading a USRYIELD output that was written according to the =3D
> following USRYIELD command<BR>
> <BR>
> USRYIELD -123.00000 201.00000-50.000000 4.0000000 5.0000000 =3D
> 1.0000000LET+0.1 <BR>
> USRYIELD 1000.0000 0.0000000 100.00000 1.5000000 0.0000000 =
=3D
> 2603.0000 & <BR>
> <BR>
> I am trying to parse the file using a C++ routine, so it is imperitive =
I =3D
> understand the details of what is written where. I followed the =
=3D
> code in the manual on reading the USRYIELD output (<A =3D
> =
HREF=3D3D"http://www.fluka.org/manual/sect/s090/text.html">http://www.flu=
ka=3D
> .org/manual/sect/s090/text.html</A>) and see that the variable =3D
> ITUSYL(NY) is supposed to correspond to the What1. When I read =
=3D
> this file the variable reads a value of 2603, which is my What12 and =
not =3D
> What1. I have taken into account the extra bytes that fortran =
=3D
> writes and all the other values are read in correctly, so is this a =
bug =3D
> or am I missing something?<BR>
> <BR>
> Thanks<BR>
> Kerry<BR>
> <BR>
> </FONT>
> </P>
>
> </BODY>
> </HTML>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0--
>
--=20
+------------------------------------------------------------------------=
----+
| Alfredo Ferrari || Tel.: +41.22.767.6119 =
|
| CERN-AB || Fax.: +41.22.767.7555 =
|
| 1211 Geneva 23 || e-mail: Alfredo.Ferrari@cern.ch =
|
| Switzerland || Alfredo.Ferrari@mi.infn.it =
|
+------------------------------------------------------------------------=
----+
------_=_NextPart_001_01C71D35.0BE304A4
Content-Type: text/html;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
charset=3Diso-8859-1">
<META NAME=3D"Generator" CONTENT=3D"MS Exchange Server version =
6.5.7651.34">
<TITLE>RE: USRYIELD binary output</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<P><FONT SIZE=3D2>Hi Alfredo,<BR>
<BR>
I have determined that the algorithm is doing things right and even =
though the what1 of the input file does not follow the manual's =
instructions it turns out to give the right results anyway. The =
difference is that instead of making "ia" a positive number =
only, if one makes it negative the algorithm produces the right result =
anyway. The following code does the work right.<BR>
<BR>
AAUSYL =3D =
0.01D+00 * DBLE (ITUSYL(NY))<BR>
IAUSYL (NY) =3D =
NINT (AAUSYL)<BR>
IEUSYL (NY) =3D =
ITUSYL (NY) - IAUSYL (NY) * 100<BR>
<BR>
If you put in ITUSYL =3D -123, then IAUSYL =3D -1 and IEUSYL =3D -23, =
but for ITUSYL =3D 77 you get IAUSYL =3D 1 and IEUSYL =3D -23, but I =
assume since IA represents a single bin the sign is ignored and when =
what1 is packed back up before being written to the USRYIELD output the =
algorithm ignores the sign again by doing<BR>
<BR>
ITUSYL(NY) =3D IEUSYL(NY) + ABS(IAUSYL(NY))*100<BR>
<BR>
Although this is a little confusing on how I was getting the right =
results while not following the explicit instructions of the manual, and =
that the written what1 was different than the supplied what1, the fault =
was mine. Even so, it might be better to not ignore the sign of IA =
when writing to the USRYIELD output so that it will agree with what the =
user input.<BR>
<BR>
Thanks<BR>
Kerry<BR>
<BR>
-----Original Message-----<BR>
From: Alfredo Ferrari [<A =
HREF=3D"mailto:alfredo.ferrari@cern.ch">mailto:alfredo.ferrari@cern.ch</A=
>]<BR>
Sent: Mon 12/11/2006 7:35 AM<BR>
To: Lee, Kerry T. (JSC-SF)[LMIT]<BR>
Cc: fluka-discuss@fluka.org<BR>
Subject: RE: USRYIELD binary output<BR>
<BR>
Hi Kerry<BR>
<BR>
the code tries to unpack from what(1) the infos about which two =
physical<BR>
quantities you want to deal with. Probably a nonsense number is =
unpacked<BR>
into something strange (I should verify).<BR>
<BR>
&=
nbsp; Ciao<BR>
&=
nbsp; Alfredo<BR>
On Fri, 8 Dec 2006, Lee, Kerry T. (JSC-SF)[LMIT] wrote:<BR>
<BR>
> This is a multi-part message in MIME format.<BR>
><BR>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0<BR>
> Content-Type: text/plain;<BR>
> charset=3D"iso-8859-1"<BR>
> Content-Transfer-Encoding: quoted-printable<BR>
><BR>
> I seem to have found that I was in error when explaining the =
problem. I =3D<BR>
> missed the fact that what12 is stored as the next variable after =
what1, =3D<BR>
> but the value stored where what1 should be is 77 and the value from =
the =3D<BR>
> input file is -123. How is this not being transferred =
correctly to the =3D<BR>
> output file?<BR>
><BR>
> Thanks<BR>
> Kerry<BR>
><BR>
><BR>
> -----Original Message-----<BR>
> From: Lee, Kerry T. (JSC-SF)[LMIT]<BR>
> Sent: Thu 12/7/2006 10:27 AM<BR>
> To: fluka-discuss@fluka.org<BR>
> Subject: USRYIELD binary output<BR>
> =3D20<BR>
><BR>
> Dear FLUKA authors,<BR>
><BR>
> I am reading a USRYIELD output that was written according to the =
=3D<BR>
> following USRYIELD command<BR>
><BR>
> USRYIELD -123.00000 201.00000-50.000000 4.0000000 5.0000000 =
=3D<BR>
> 1.0000000LET+0.1 =3D20<BR>
> USRYIELD 1000.0000 0.0000000 100.00000 1.5000000 =
0.0000000 2603.0000 & =3D<BR>
> =3D20<BR>
> =3D20<BR>
> I am trying to parse the file using a C++ routine, so it is =
imperitive I =3D<BR>
> understand the details of what is written where. I followed =
the code in =3D<BR>
> the manual on reading the USRYIELD output =3D<BR>
> (<A =
HREF=3D"http://www.fluka.org/manual/sect/s090/text.html">http://www.fluka=
.org/manual/sect/s090/text.html</A>) and see that the =3D<BR>
> variable ITUSYL(NY) is supposed to correspond to the What1. =
When I read =3D<BR>
> this file the variable reads a value of 2603, which is my What12 =
and not =3D<BR>
> What1. I have taken into account the extra bytes that fortran =
writes =3D<BR>
> and all the other values are read in correctly, so is this a bug or =
am I =3D<BR>
> missing something?<BR>
><BR>
> Thanks<BR>
> Kerry<BR>
><BR>
><BR>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0<BR>
> Content-Type: text/html;<BR>
> charset=3D"iso-8859-1"<BR>
> Content-Transfer-Encoding: quoted-printable<BR>
><BR>
> <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML =
3.2//EN"><BR>
> <HTML><BR>
> <HEAD><BR>
> <META HTTP-EQUIV=3D3D"Content-Type" =
CONTENT=3D3D"text/html; =3D<BR>
> charset=3D3Diso-8859-1"><BR>
> <META NAME=3D3D"Generator" CONTENT=3D3D"MS =
Exchange Server version =3D<BR>
> 6.5.7651.34"><BR>
> <TITLE>RE: USRYIELD binary output</TITLE><BR>
> </HEAD><BR>
> <BODY><BR>
> <!-- Converted from text/plain format --><BR>
><BR>
> <P><FONT SIZE=3D3D2>I seem to have found that I was in =
error when =3D<BR>
> explaining the problem.&nbsp; I missed the fact that what12 is =
stored as =3D<BR>
> the next variable after what1, but the value stored where what1 =
should =3D<BR>
> be is 77 and the value from the input file is -123.&nbsp; How =
is this =3D<BR>
> not being transferred correctly to the output file?<BR><BR>
> <BR><BR>
> Thanks<BR><BR>
> Kerry<BR><BR>
> <BR><BR>
> <BR><BR>
> -----Original Message-----<BR><BR>
> From: Lee, Kerry T. (JSC-SF)[LMIT]<BR><BR>
> Sent: Thu 12/7/2006 10:27 AM<BR><BR>
> To: fluka-discuss@fluka.org<BR><BR>
> Subject: USRYIELD binary output<BR><BR>
> <BR><BR>
> <BR><BR>
> Dear FLUKA authors,<BR><BR>
> <BR><BR>
> I am reading a USRYIELD output that was written according to the =
=3D<BR>
> following USRYIELD command<BR><BR>
> <BR><BR>
> USRYIELD&nbsp; -123.00000 201.00000-50.000000 4.0000000 =
5.0000000 =3D<BR>
> 1.0000000LET+0.1&nbsp;&nbsp;<BR><BR>
> USRYIELD&nbsp;&nbsp; 1000.0000 0.0000000 100.00000 =
1.5000000 0.0000000 =3D<BR>
> 2603.0000 =
&amp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&=
;nbsp;<BR><BR>
> <BR><BR>
> I am trying to parse the file using a C++ routine, so it is =
imperitive I =3D<BR>
> understand the details of what is written where.&nbsp; I =
followed the =3D<BR>
> code in the manual on reading the USRYIELD output (<A =3D<BR>
> HREF=3D3D"<A =
HREF=3D"http://www.fluka.org/manual/sect/s090/text.html">http://www.fluka=
.org/manual/sect/s090/text.html</A>"><A =
HREF=3D"http://www.fluka=3D">http://www.fluka=3D><BR>
> .org/manual/sect/s090/text.html</A>) and see that the =
variable =3D<BR>
> ITUSYL(NY) is supposed to correspond to the What1.&nbsp; When I =
read =3D<BR>
> this file the variable reads a value of 2603, which is my What12 =
and not =3D<BR>
> What1.&nbsp; I have taken into account the extra bytes that =
fortran =3D<BR>
> writes and all the other values are read in correctly, so is this a =
bug =3D<BR>
> or am I missing something?<BR><BR>
> <BR><BR>
> Thanks<BR><BR>
> Kerry<BR><BR>
> <BR><BR>
> </FONT><BR>
> </P><BR>
><BR>
> </BODY><BR>
> </HTML><BR>
> ------_=3D_NextPart_001_01C71AD8.4F96A1D0--<BR>
><BR>
<BR>
--<BR>
<BR>
+------------------------------------------------------------------------=
----+<BR>
| Alfredo =
Ferrari =
|| Tel.: =
+41.22.767.6119 &nbs=
p; |<BR>
| =
CERN-AB =
=
|| Fax.: =
+41.22.767.7555 &nbs=
p; |<BR>
| 1211 Geneva =
23  =
; || e-mail: =
Alfredo.Ferrari@cern.ch |<BR>
| =
Switzerland &n=
bsp; =
|| =
Alfredo.Ferrari@mi.infn.it |<BR>
+------------------------------------------------------------------------=
----+<BR>
<BR>
</FONT>
</P>
</BODY>
</HTML>
------_=_NextPart_001_01C71D35.0BE304A4--
This archive was generated by hypermail 2.1.6 : Wed Dec 13 2006 - 13:59:05 CET