00001 /* 00002 * ******** fete: From ENDF To ENDL ********* 00003 * 00004 * Copyright (c) 2006, The Regents of the University of California. 00005 * All rights reserved. 00006 * 00007 * Produced at the Lawrence Livermore National Laboratory. 00008 * Written by David A. Brown, Gerry Hedstrom, Tony Hill 00009 * 00010 * This file is part of fete v1.0 (UCRL-CODE-218718) 00011 * 00012 * Please read the COPYING file for "Our Notice and GNU General 00013 * Public License" in the root of this software distribution. 00014 * 00015 * This program is free software; you can redistribute it and/or modify 00016 * it under the terms of the GNU General Public License (as published by 00017 * the Free Software Foundation) version 2, dated June 1991. 00018 * 00019 * This program is distributed in the hope that it will be useful, 00020 * but WITHOUT ANY WARRANTY; without even the IMPLIED WARRANTY OF 00021 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the terms 00022 * and conditions of the GNU General Public License for more details. 00023 * 00024 * You should have received a copy of the GNU General Public License along 00025 * with this program; if not, write to the Free Software Foundation, Inc., 00026 * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00027 * 00028 * $Revision: 1867 $ 00029 * $Date: 2006-05-15 10:19:04 -0700 (Mon, 15 May 2006) $ 00030 * $Author: dbrown $ 00031 * $Id: endl_precision.hpp 1867 2006-05-15 17:19:04Z dbrown $ 00032 * 00033 * ******** fete: From ENDF To ENDL ********* 00034 */ 00035 00036 #ifndef __ENDL_PRECISION_H 00037 #define __ENDL_PRECISION_H 00038 00039 #include "global_params.hpp" 00040 00041 extern GlobalParameterClass Global; 00042 00043 //! The smallest energy difference, but make it positive 00044 double ENDL_EPSILON( double E ); 00045 00046 //! The width of a jump 00047 inline double ENDL_JUMP_WIDTH( double E ) 00048 { return Global.Value("jump_factor")*ENDL_EPSILON( E ); } 00049 //{ return 2.0*ENDL_EPSILON( E ); } 00050 00051 //! Make a delta function wider than a jump. 00052 inline double DELTA_WIDTH( double E ) 00053 { return 1.1*Global.Value("jump_factor")*ENDL_EPSILON( E ); } 00054 //{ return 2.2*ENDL_EPSILON( E ); } 00055 00056 #endif