[mesa-users] Test Burn Issues

Ryan Orvedahl ryan.orvedahl at stonybrook.edu
Sat Jun 29 20:23:30 EDT 2013


**
 Hi All,

 I am exploring the net/test routines in hopes to use the network in an
external program. When I run the test_net.f program, everything seems to
work okay, but when I compile in debug mode it crashes.

 In the test_burn.f Do_One_Test_Burn subroutine I setup my initial
conditions as:

else if (net_file == 'regular-9.net <http://h_burn.net>' .or. net_file == '
extended-33.net') then

               ! C/O burn
               times(1) = 1.d-2
               h = 1d-14
               logT = 9.d0
               logRho = 8.d0
               xin = 0
               xin(net_iso(ic12)) =  0.5d0
               xin(net_iso(io16)) =  0.5d0

               rtol(:) = 1d-12
               atol(:) = 1d-14

               screening_mode = extended_screening

 and in the mod_test_net.f module, I call Do_One_Test_Burn as:

if (.true.) then
    which_rates_choice = rates_JR_if_available
    call Do_One_Test_Burn('regular-9.net <http://h_burn.net>', .false.)
   !call Do_One_Test_Burn('extended-33.net', .false.)
   call test_net_cleanup
   call net_shutdown
    return
 endif

I added the following debugging flags in the makefile:

%.o:  $(TEST_SRC_DIR)/%.f
    $(TEST_COMPILE) $(FCfixed) -fno-range-check -O1 -fbounds-check \
      -fbacktrace -Wunitialized -Wunused -ffpe-trap=invalid,zero,overflow $<

The following is the error I get:

Program received signal SIGFPE: Floating-point exception - erroneous
arithmetic operation.

Backtrace for this error:
#0  0x7F9D9641CFD7
#1  0x7F9D9641D5A4
#2  0x36BE0359AF
#3  0x3DD4C26443
#4  0x49E86E in __screen5_MOD_screen5_init_az_info
#5  0x44AA88 in __net_screen_MOD_screening_pair
#6  0x44AF67 in __net_screen_MOD_eval_screen_pair
#7  0x44BC8B in __net_screen_MOD_screen_net
#8  0x44C08A in __net_screen_MOD_make_screening_tables
#9  0x43E688 in __net_initialize_MOD_alloc_net_general_info at
net_initialize.f:197
#10  0x411781 in __test_net_support_MOD_test_net_setup at
test_net_support.f:140
#11  0x41ED36 in __test_burn_MOD_do_one_test_burn at test_burn.f:135
#12  0x42CFFF in __mod_test_net_MOD_test at mod_test_net.f:90
#13  0x42D03F in test_net at test_net.f:29
#14  0x36BE021734
Floating point exception

The issue seems to be caused by calling screen5_init_AZ_info with z1 = 0 or
z2 = 0 which leads to trying to evaluate lzav = x53 * log(0)

Does anyone have any ideas on how to fix this?

I am using MESA version 4942 without the SDK. My fortran compiler is gcc
version 4.7.2 20120921. I have included both net files.

Thanks,
Ryan

-- 
Ryan Orvedahl
Undergrad Student
Stony Brook University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130629/612c97bd/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: regular-9.net
Type: application/octet-stream
Size: 234 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130629/612c97bd/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: extended-33.net
Type: application/octet-stream
Size: 660 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130629/612c97bd/attachment-0001.obj>


More information about the Mesa-users mailing list