Thursday, August 30, 2007

Removing values of names set using xset

Jeremy Sanders enquires how to return to the default after having used xset to change APECROOT. Doing
XSPEC12> xset APECROOT
actually works although this is undocumented and is misleading because doing an xset then shows APECROOT still existing but with no value. Fixed so that the name no longer shows up as existing.

Tuesday, August 28, 2007

Abundances in NEI models

Frank Haberl pointed out that the documentation on the NEI models is wrong. The documentation states that abundances are wrt Anders & Grevesse while in fact they are relative to whatever is set by the abund command. This was a change in code from June 2005 which was never documented. Updated all the relevant files (11.3.2ai).

Thursday, August 23, 2007

PGPLOT PS driver

Phil Evans reports occasionally losing a few pixels from text labels due to the bounding box being set too small. Made a simple fix in grtext.f to add 2% in each direction.

rdblur model

The hard lower limit on the Rin parameter for rdblur should be 6, not 0 as currently set. If Rin < 6 then the model returns zero.

Monday, August 20, 2007

xspec updates

Bug fixes/updates from Craig. These are all listed on the bugs page.

12.3.1ak This fixes dummyrsp so it can be applied to any model for a given spectrum, and can now take advantage of XSPEC's recently added multiple sources capability in the same way as the response command. Previously dummyrsp only worked properly for the standard single source cases. This also fixes a Solaris specific bug, which could cause NaN calculation errors when folding a model with a dummy response. Report added on Aug 20, 2007.

The fix is to replace headas/Xspec/src/XSModel/Data/ SpectralData.cxx, src/XSModel/Data/Detector/ UserDummyResponse.cxx, src/XSModel/GlobalContainer/ DataContainer.cxx, src/XSUser/Handler/ xsDummyrsp.cxx, xsResponse.cxx, and rebuild xspec.

12.3.1al Using delcomp to remove a table model component causes a segmentation fault. If editmod is used on a model with a table component, and the component is written with an abbreviated form of "atable", "mtable", etc., editmod erroneously considers it a modified component. Report added on Aug 20, 2007.

The fix is to replace headas/Xspec/src/XSModel/Model/ Model.cxx, src/XSModel/Model/Component/ Component.cxx, and rebuild xspec.

12.3.1am A couple of output formatting improvements: lumin output should be in standard "d.ddde+EE" format and not written as "d.ddd x 10^EE". The margin command output should switch from fixed to scientific format for a wider range of exponent values. Our thanks to Jeremy Sanders for pointing this out. Report added on Aug 20, 2007.

The fix is to replace headas/Xspec/src/XSFit/MCMC/ MarginGrid.cxx, src/XSModel/Model/ Model.cxx, and rebuild xspec.

Thursday, August 16, 2007

building lisatools

Building lisatools using the instructions provided by Michele. A few things to note :

1. The correct svn command is
svn co https://lisatools.svn.sourceforge.net/svnroot/lisatools lisatools

2. The build of synthLISA fails with a C++ error. This is because python2.5 requires the
SWIGged files to have been created by swig 1.3.31, not the swig 1.3.29 in lisatools. Grabbed
the new gzipped swig from sourceforge and replaced in lisatools/Packages. Then edited
master-install.py to use 1.3.31 instead of 1.3.29.

3. This still fails because the SWIGged files are not created in the standard build. Deleted
lisasim/{lisasim-swig_wrap.cpp, lisasim.py} then built synthLISA using the setup.py in
synthLISA-1.3.3. Note that if I rerun the master-install.py it reextracts the contents of the
tarfile and puts back the two files I just deleted.

Monday, August 13, 2007

xspec lumin command

Jeremy Sanders points out that the lumin command doesn't write the number in standard computer notation (x.xxxxxe+yy) unlike every other output number in xspec. So, modified Model.cxx appropriately.

xselect array sizes

Kenji ran into a problem with xselect trying to use a large number of MKFs. Since the cost in memory is low I increased most of the sizes in the xsel.inc file.