Commits June 2010

commits@lists.einsteintoolkit.org

12 participants
242 discussions

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 129)
by cott＠tapir.caltech.edu 29 Jun '10

29 Jun '10

User: cott Date: 2010/06/30 01:49 AM Modified: /trunk/src/ GRHydro_PPM.F90 Log: * remove unused EOS code from GRHydro_PPM.F90 File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_PPM.F90 Delta lines: +0 -18 =================================================================== --- trunk/src/GRHydro_PPM.F90 2010-06-26 13:47:35 UTC (rev 128) +++ trunk/src/GRHydro_PPM.F90 2010-06-30 06:49:09 UTC (rev 129) @@ -64,11 +64,6 @@ DECLARE_CCTK_PARAMETERS -#ifdef _EOS_BASE_INC_ -#undef _EOS_BASE_INC_ -#endif -#include "EOS_Base.inc" - CCTK_INT :: handle,poly,nx CCTK_REAL :: dx CCTK_REAL, dimension(nx) :: rho,velx,vely,velz,eps @@ -109,19 +104,6 @@ trivial_rp = .true. -!!$ Loop to set the pressure and eps -!!$ poly=0 -!!$ if (poly .eq. 0) then -!!$ do i = 1, nx -!!$ press(i) = EOS_Pressure(handle, rho(i), eps(i)) -!!$ end do -!!$ else -!!$ do i = 1, nx -!!$ press(i) = EOS_Pressure(handle, rho(i), eps(i)) -!!$ eps(i) = EOS_SpecificIntEnergy(handle, rho(i), press(i)) -!!$ end do -!!$ end if - !!$ Average slopes delta_m(a). See (1.7) of Colella and Woodward, p.178 !!$ This is the expression for an even grid.

1 0

[svn:einsteintoolkit] incoming/EOS_Omni/ (Rev. 1)
by cott＠tapir.caltech.edu 29 Jun '10

29 Jun '10

User: cott Date: 2010/06/29 09:59 AM Added: /EOS_Omni/ README, interface.ccl, param.ccl, schedule.ccl /EOS_Omni/doc/ documentation.tex /EOS_Omni/par/ /EOS_Omni/src/ EOS_Omni_Module.F90, EOS_Omni_SingleVarCalls.F90, EOS_Omni_Startup.F90, make.code.defn /EOS_Omni/test/ Log: * EOS_Omni thorn: Initial implementation providing a polytrope and a gamma-law EOS for testing. File Changes: Directory: /EOS_Omni/ ===================== File [added]: README Delta lines: +24 -0 =================================================================== --- EOS_Omni/README (rev 0) +++ EOS_Omni/README 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,24 @@ +Cactus Code Thorn EOS_Omni +Author(s) : Christian D. Ott <cott(a)tapir.caltech.edu> +Maintainer(s): Christian D. Ott <cott(a)tapir.caltech.edu> +Licence : LGPL +-------------------------------------------------------------------------- + +This Thorn provides a one-fits-all EOS handler. +All EOS are actually implemented right inside this single +thorn without crazy interfaces and lots of function calls. + +eoskey controls the type of EOS being used: + +eoskey: + +1 --- polytropic EOS +2 --- gamma-law EOS +3 --- hybrid EOS +4 --- finite-T microphysical NSE EOS + + +keyerr settings: + +-1 -- called EOS routine with keytemp=1, but + routine does not allow such a call File [added]: interface.ccl Delta lines: +89 -0 =================================================================== --- EOS_Omni/interface.ccl (rev 0) +++ EOS_Omni/interface.ccl 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,89 @@ +# Interface definition for thorn EOS_Omni + +implements: EOS_Omni + +inherits: + +void FUNCTION EOS_Omni_press(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY press, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_press WITH EOS_Omni_EOS_Press LANGUAGE Fortran + + +void FUNCTION EOS_Omni_DPressByDEps(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY DPressByDEps, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_DPressByDEps WITH EOS_Omni_EOS_DPressByDEps LANGUAGE Fortran + + +void FUNCTION EOS_Omni_DPressByDRho(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY DPressByDRho, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_DPressByDRho WITH EOS_Omni_EOS_DPressByDRho LANGUAGE Fortran + + +void FUNCTION EOS_Omni_cs2(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL OUT ARRAY cs2, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_cs2 WITH EOS_Omni_EOS_cs2 LANGUAGE Fortran + + +void FUNCTION EOS_Omni_EpsFromPress(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL IN ARRAY rho, \ + CCTK_REAL INOUT ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL IN ARRAY press, \ + CCTK_REAL OUT ARRAY xeps, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_EpsFromPress WITH EOS_Omni_EOS_eps_from_press LANGUAGE Fortran + + +void FUNCTION EOS_Omni_RestMassDensityFromEpsPress(CCTK_INT IN eoskey, \ + CCTK_INT IN havetemp, \ + CCTK_INT IN npoints, \ + CCTK_REAL OUT ARRAY rho, \ + CCTK_REAL IN ARRAY eps, \ + CCTK_REAL INOUT ARRAY temp, \ + CCTK_REAL IN ARRAY ye, \ + CCTK_REAL IN ARRAY press, \ + CCTK_INT OUT ARRAY keyerr, \ + CCTK_INT OUT anyerr) + +PROVIDES FUNCTION EOS_Omni_RestMassDensityFromEpsPress WITH EOS_Omni_EOS_RestMassDensityFromEpsPress LANGUAGE Fortran File [added]: param.ccl Delta lines: +26 -0 =================================================================== --- EOS_Omni/param.ccl (rev 0) +++ EOS_Omni/param.ccl 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,26 @@ +# Parameter definitions for thorn EOS_Omni + + +restricted: + +# poly EOS +CCTK_REAL poly_gamma "Adiabatic Index for poly EOS" +{ + : :: "" +} 2.0 + +REAL poly_k "Polytropic constant in c=G=Msun=1" +{ + : :: "" +} 100.0 + +# gamma-law EOS +CCTK_REAL gl_gamma "Adiabatic Index for gamma-law EOS" +{ + : :: "" +} 2.0 + +REAL gl_k "Polytropic constant in c=G=Msun=1 for gamma-law EOS" +{ + : :: "" +} 100.0 File [added]: schedule.ccl Delta lines: +8 -0 =================================================================== --- EOS_Omni/schedule.ccl (rev 0) +++ EOS_Omni/schedule.ccl 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,8 @@ +# Schedule definitions for thorn EOS_Omni + + +schedule EOS_Omni_Startup AT WRAGH +{ + LANG: Fortran + OPTIONS: global +} "Set up conversion factors and other fun stuff" Directory: /EOS_Omni/doc/ ========================= File [added]: documentation.tex Delta lines: +144 -0 =================================================================== --- EOS_Omni/doc/documentation.tex (rev 0) +++ EOS_Omni/doc/documentation.tex 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,144 @@ +% *======================================================================* +% Cactus Thorn template for ThornGuide documentation +% Author: Ian Kelley +% Date: Sun Jun 02, 2002 +% $Header$ +% +% Thorn documentation in the latex file doc/documentation.tex +% will be included in ThornGuides built with the Cactus make system. +% The scripts employed by the make system automatically include +% pages about variables, parameters and scheduling parsed from the +% relevant thorn CCL files. +% +% This template contains guidelines which help to assure that your +% documentation will be correctly added to ThornGuides. More +% information is available in the Cactus UsersGuide. +% +% Guidelines: +% - Do not change anything before the line +% % START CACTUS THORNGUIDE", +% except for filling in the title, author, date, etc. fields. +% - Each of these fields should only be on ONE line. +% - Author names should be separated with a \\ or a comma. +% - You can define your own macros, but they must appear after +% the START CACTUS THORNGUIDE line, and must not redefine standard +% latex commands. +% - To avoid name clashes with other thorns, 'labels', 'citations', +% 'references', and 'image' names should conform to the following +% convention: +% ARRANGEMENT_THORN_LABEL +% For example, an image wave.eps in the arrangement CactusWave and +% thorn WaveToyC should be renamed to CactusWave_WaveToyC_wave.eps +% - Graphics should only be included using the graphicx package. +% More specifically, with the "\includegraphics" command. Do +% not specify any graphic file extensions in your .tex file. This +% will allow us to create a PDF version of the ThornGuide +% via pdflatex. +% - References should be included with the latex "\bibitem" command. +% - Use \begin{abstract}...\end{abstract} instead of \abstract{...} +% - Do not use \appendix, instead include any appendices you need as +% standard sections. +% - For the benefit of our Perl scripts, and for future extensions, +% please use simple latex. +% +% *======================================================================* +% +% Example of including a graphic image: +% \begin{figure}[ht] +% \begin{center} +% \includegraphics[width=6cm]{MyArrangement_MyThorn_MyFigure} +% \end{center} +% \caption{Illustration of this and that} +% \label{MyArrangement_MyThorn_MyLabel} +% \end{figure} +% +% Example of using a label: +% \label{MyArrangement_MyThorn_MyLabel} +% +% Example of a citation: +% \cite{MyArrangement_MyThorn_Author99} +% +% Example of including a reference +% \bibitem{MyArrangement_MyThorn_Author99} +% {J. Author, {\em The Title of the Book, Journal, or periodical}, 1 (1999), +% 1--16. {\tt http://www.nowhere.com/}} +% +% *======================================================================* + +% If you are using CVS use this line to give version information +% $Header$ + +\documentclass{article} + +% Use the Cactus ThornGuide style file +% (Automatically used from Cactus distribution, if you have a +% thorn without the Cactus Flesh download this from the Cactus +% homepage at www.cactuscode.org) +\usepackage{../../../../doc/latex/cactus} + +\begin{document} + +% The author of the documentation +\author{Christian D. Ott \textless cott(a)tapir.caltech.edu\textgreater} + +% The title of the document (not necessarily the name of the Thorn) +\title{EOS\_Omni} + +% the date your document was last changed, if your document is in CVS, +% please use: +% \date{$ $Date: 2004-01-07 12:12:39 -0800 (Wed, 07 Jan 2004) $ $} +\date{June 27 2010} + +\maketitle + +% Do not delete next line +% START CACTUS THORNGUIDE + +% Add all definitions used in this documentation here +% \def\mydef etc + +% Add an abstract for this thorn's documentation +\begin{abstract} + +\end{abstract} + +% The following sections are suggestive only. +% Remove them or add your own. + +\section{Introduction} + +\section{Physical System} + +\section{Numerical Implementation} + +\section{Using This Thorn} + +\subsection{Obtaining This Thorn} + +\subsection{Basic Usage} + +\subsection{Special Behaviour} + +\subsection{Interaction With Other Thorns} + +\subsection{Examples} + +\subsection{Support and Feedback} + +\section{History} + +\subsection{Thorn Source Code} + +\subsection{Thorn Documentation} + +\subsection{Acknowledgements} + + +\begin{thebibliography}{9} + +\end{thebibliography} + +% Do not delete next line +% END CACTUS THORNGUIDE + +\end{document} Directory: /EOS_Omni/src/ ========================= File [added]: EOS_Omni_Module.F90 Delta lines: +22 -0 =================================================================== --- EOS_Omni/src/EOS_Omni_Module.F90 (rev 0) +++ EOS_Omni/src/EOS_Omni_Module.F90 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,22 @@ +module EOS_Omni_Module + + implicit none + + real*8,parameter :: rho_gf = 1.61930347d-18 + real*8,parameter :: press_gf = 1.80171810d-39 + real*8,parameter :: eps_gf = 1.11265006d-21 + real*8,parameter :: time_gf = 2.03001708d+05 + real*8,parameter :: mass_gf = 5.02765209d-34 + real*8,parameter :: length_gf = 6.77140812d-06 + + real*8,parameter :: inv_rho_gf = 6.17549470205236d17 + real*8,parameter :: inv_press_gf = 5.55025783445257d38 + real*8,parameter :: inv_eps_gf = 8.98755175549085d20 + real*8,parameter :: inv_time_gf = 4.92606692747629d-6 + real*8,parameter :: inv_mass_gf = 1.98899999860571d33 + real*8,parameter :: inv_length_gf = 147679.77092481d0 + + real*8 :: poly_k_cgs = 0.0d0 + real*8 :: gl_k_cgs = 0.0d0 + +end module EOS_Omni_Module File [added]: EOS_Omni_SingleVarCalls.F90 Delta lines: +333 -0 =================================================================== --- EOS_Omni/src/EOS_Omni_SingleVarCalls.F90 (rev 0) +++ EOS_Omni/src/EOS_Omni_SingleVarCalls.F90 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,333 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + +! eoskey: +! 1 --- polytropic EOS +! 2 --- gamma-law EOS +! 3 --- hybrid EOS +! 4 --- finite-T microphysical NSE EOS + +subroutine EOS_Omni_EOS_Press(eoskey,keytemp,npoints,& + rho,eps,temp,ye,press,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: press(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + press(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**poly_gamma + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + press(i) = (gl_gamma - 1.0d0) * rho(i) * eps(i) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_Press + + + +subroutine EOS_Omni_EOS_DPressByDRho(eoskey,keytemp,npoints,& + rho,eps,temp,ye,dpdrhoe,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: dpdrhoe(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdrhoe(i) = press_gf * poly_k_cgs * & + poly_gamma * inv_rho_gf * & + (rho(i)*inv_rho_gf) ** (poly_gamma - 1.d0) + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdrhoe(i) = (gl_gamma-1.0d0) * & + eps(i) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_DPressByDRho + +subroutine EOS_Omni_EOS_DPressByDEps(eoskey,keytemp,npoints,& + rho,eps,temp,ye,dpdepsrho,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: dpdepsrho(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdepsrho(i) = 0.0d0 + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + dpdepsrho(i) = (gl_gamma - 1.0d0) * & + rho(i) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_DPressByDEps + +subroutine EOS_Omni_EOS_cs2(eoskey,keytemp,npoints,& + rho,eps,temp,ye,cs2,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: cs2(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + real*8 :: xpress,xdpdrhoe,xdpderho + + anyerr = 0 + keyerr(:) = 0 + + select case (eoskey) + case (1) + ! polytropic EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) / & + (poly_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + xpress = press_gf*poly_k_cgs * & + (rho(i)*inv_rho_gf)**(poly_gamma) + cs2(i) = poly_gamma * xpress / rho(i) / & + (1 + eps(i) + xpress/rho(i)) + enddo + case (2) + ! gamma-law EOS + if(keytemp.eq.1) then + do i=1,npoints + eps(i) = press_gf * gl_k_cgs * & + (rho(i)*inv_rho_gf)**(gl_gamma) / & + (gl_gamma - 1.0d0) / rho(i) + enddo + endif + do i=1,npoints + xpress = (gl_gamma-1.0d0)*rho(i)*eps(i) + xdpdrhoe = (gl_gamma-1.0d0)*eps(i) + xdpderho = (gl_gamma-1.0d0)*rho(i) + cs2(i) = (xdpdrhoe + xpress * xdpderho / (rho(i)**2)) / & + (1.0d0 + eps(i) + xpress/rho(i)) + enddo + + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_cs2 + + +subroutine EOS_Omni_EOS_eps_from_press(eoskey,keytemp,npoints,& + rho,eps,temp,ye,press,xeps,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: rho(npoints),ye(npoints),press(npoints) + CCTK_REAL, intent(inout) :: eps(npoints), temp(npoints) + CCTK_REAL, intent(out) :: xeps(npoints) + + ! local vars + integer :: i + character(256) :: warnstring + + if(keytemp.eq.1) then + anyerr = 1 + keyerr(:) = -1 + else + anyerr = 0 + keyerr(:) = 0 + endif + + select case (eoskey) + case (1) + ! polytropic EOS + do i=1,npoints + xeps(i) = press(i) / (poly_gamma - 1.0d0) / rho(i) + eps(i) = xeps(i) + enddo + case (2) + ! gamma-law EOS + do i=1,npoints + xeps(i) = press(i) / (gl_gamma - 1.0d0) / rho(i) + eps(i) = xeps(i) + enddo + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_eps_from_press + + +subroutine EOS_Omni_EOS_RestMassDensityFromEpsPress(eoskey,keytemp,npoints,& + rho,eps,temp,ye,press,keyerr,anyerr) + + use EOS_Omni_Module + implicit none + DECLARE_CCTK_PARAMETERS + + CCTK_INT, intent(in) :: eoskey,keytemp,npoints + CCTK_INT, intent(out) :: keyerr(npoints) + CCTK_INT, intent(out) :: anyerr + CCTK_REAL, intent(in) :: ye(npoints),press(npoints),eps(npoints) + CCTK_REAL, intent(inout) :: rho(npoints),temp(npoints) + + + ! local vars + integer :: i + character(256) :: warnstring + + if(keytemp.eq.1) then + anyerr = 1 + keyerr(:) = -1 + else + anyerr = 0 + keyerr(:) = 0 + endif + + select case (eoskey) + case (1) + ! polytropic EOS + do i=1,npoints + rho(i) = press(i) / ((poly_gamma - 1.0d0)*eps(i)) + enddo + case (2) + ! gamma-law EOS + do i=1,npoints + rho(i) = press(i) / ((gl_gamma - 1.0d0)*eps(i)) + enddo + case DEFAULT + write(warnstring,*) "eoskey ",eoskey," not implemented!" + call CCTK_WARN(0,warnstring) + end select + +end subroutine EOS_Omni_EOS_RestMassDensityFromEpsPress File [added]: EOS_Omni_Startup.F90 Delta lines: +18 -0 =================================================================== --- EOS_Omni/src/EOS_Omni_Startup.F90 (rev 0) +++ EOS_Omni/src/EOS_Omni_Startup.F90 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,18 @@ +#include "cctk.h" +#include "cctk_Parameters.h" +#include "cctk_Arguments.h" +#include "cctk_Functions.h" + +subroutine EOS_Omni_Startup(CCTK_ARGUMENTS) + + use EOS_Omni_Module + implicit none + + DECLARE_CCTK_PARAMETERS + DECLARE_CCTK_ARGUMENTS + + poly_k_cgs = poly_k * rho_gf**poly_gamma / press_gf + + gl_k_cgs = gl_k * rho_gf**gl_gamma / press_gf + +end subroutine EOS_Omni_Startup File [added]: make.code.defn Delta lines: +7 -0 =================================================================== --- EOS_Omni/src/make.code.defn (rev 0) +++ EOS_Omni/src/make.code.defn 2010-06-29 14:59:02 UTC (rev 1) @@ -0,0 +1,7 @@ +# Main make.code.defn file for thorn EOS_Omni + +# Source files in this directory +SRCS = EOS_Omni_Module.F90 EOS_Omni_Startup.F90 EOS_Omni_SingleVarCalls.F90 + +# Subdirectories containing source files +SUBDIRS =

1 0

[svn:einsteintoolkit] www/info/support/ (Rev. 308)
by gallen＠cct.lsu.edu 28 Jun '10

28 Jun '10

User: gallen Date: 2010/06/28 08:55 AM Modified: /info/support/ index.php Log: Tweak to phone details File Changes: Directory: /info/support/ ========================= File [modified]: index.php Delta lines: +2 -2 =================================================================== --- info/support/index.php 2010-06-27 04:34:14 UTC (rev 307) +++ info/support/index.php 2010-06-28 13:55:19 UTC (rev 308) @@ -38,8 +38,8 @@ call is 10am (central US time) every Monday. Instructions how to dial in are given in the following: <h4>Dial-in Instructions</h4> -Each participant of the teleconference will call 225-578-4942 -(8-4942 from LSU-campus or 866-573-0359 toll free). +Each participant of the teleconference will call (+1) 225-578-4942 +((+1) 866-573-0359 toll free from inside the USA). They will hear: <pre>

1 0

[svn:einsteintoolkit] www/release-info/ (Rev. 307)
by knarf＠cct.lsu.edu 26 Jun '10

26 Jun '10

User: knarf Date: 2010/06/26 11:34 PM Modified: /release-info/ pelican__1.log, pelican__2.log Log: new testsuite results for pelican File Changes: Directory: /release-info/ ========================= File [modified]: pelican__1.log Delta lines: +54 -41 =================================================================== --- release-info/pelican__1.log 2010-06-24 14:43:05 UTC (rev 306) +++ release-info/pelican__1.log 2010-06-27 04:34:14 UTC (rev 307) @@ -172,37 +172,6 @@ Missing: IDPUNCTURE TATJACOBI TATMG TATPETSC - Tests missed for different number of processors required: - - nonstaggered in CarpetReduce - (Carpet/CarpetReduce/test/nonstaggered.par) - Requires 1 processors - - staggered in CarpetReduce - (Carpet/CarpetReduce/test/staggered.par) - Requires 1 processors - - bhns_eval in TwoPunctures - (EinsteinInitialData/TwoPunctures/test/bhns_eval.par) - Requires 1 processors - - bhns_interp in TwoPunctures - (EinsteinInitialData/TwoPunctures/test/bhns_interp.par) - Requires 1 processors - - twopunctures in TwoPunctures - (EinsteinInitialData/TwoPunctures/test/twopunctures.par) - Requires 1 processors - - teukolsky in WeylScal4 - (EinsteinAnalysis/WeylScal4/test/teukolsky.par) - Requires 1 processors - - teukolskyID in WeylScal4 - (EinsteinAnalysis/WeylScal4/test/teukolskyID.par) - Requires 1 processors - - Tests with different number of test files: Multipole (test_22) @@ -232,6 +201,8 @@ AHFinderDirect [15] + CarpetReduce [2] + DistortedBHIVP [1] Exact [3] @@ -262,8 +233,12 @@ TOVSolver [5] + TwoPunctures [3] + WeylScal4 [2] + + Details: @@ -325,6 +300,12 @@ recoverML + CarpetReduce: + + nonstaggered + + staggered + DistortedBHIVP: E2xeon_test_dbh @@ -471,8 +452,22 @@ test_two_max + TwoPunctures: + bhns_eval + bhns_interp + + twopunctures + + WeylScal4: + + teukolsky + + teukolskyID + + + Thorns with no valid testsuite parameter files: @@ -494,7 +489,6 @@ CarpetInterp2 CarpetLib CarpetMask - CarpetReduce CarpetRegrid CarpetRegrid2 CarpetSlab @@ -535,6 +529,7 @@ IOUtil InitBase LAPACK + LORENE LegoExcision LocalInterp LocalReduce @@ -545,6 +540,9 @@ ML_BSSN_Helper ML_BSSN_O2 ML_BSSN_O2_Helper + Meudon_Bin_BH + Meudon_Bin_NS + Meudon_Mag_NS MoL NaNChecker NewRad @@ -574,8 +572,6 @@ Time TimerReport TmunuBase - TwoPunctures - WeylScal4 jpeg6b @@ -634,6 +630,10 @@ Success: 55 files compared, 9 differ in the last digits AHFinderDirect: recoverML Success: 5 files compared, 4 differ in the last digits + CarpetReduce: nonstaggered + Failure: 5 files compared, 1 differ, 1 differ significantly + CarpetReduce: staggered + Failure: 5 files compared, 1 differ, 1 differ significantly DistortedBHIVP: E2xeon_test_dbh Failure: 100 files compared, 44 differ, 44 differ significantly Exact: Schwarzschild_EF @@ -750,6 +750,16 @@ Success: 232 files compared, 8 differ in the last digits TOVSolver: test_two_max Success: 232 files compared, 8 differ in the last digits + TwoPunctures: bhns_eval + Success: 232 files compared, 55 differ in the last digits + TwoPunctures: bhns_interp + Success: 232 files compared, 55 differ in the last digits + TwoPunctures: twopunctures + Success: 144 files compared, 93 differ in the last digits + WeylScal4: teukolsky + Failure: 8 files compared, 8 differ, 8 differ significantly + WeylScal4: teukolskyID + Success: 8 files identical ======================================================================== @@ -757,19 +767,22 @@ Summary for configuration sim Total available tests -> 131 - Unrunnable tests -> 48 - Runnable tests -> 83 - Total number of thorns -> 120 - Number of tested thorns -> 19 - Number of tests passed -> 81 + Unrunnable tests -> 41 + Runnable tests -> 90 + Total number of thorns -> 124 + Number of tested thorns -> 22 + Number of tests passed -> 85 Number passed only to - set tolerance -> 66 - Number failed -> 2 + set tolerance -> 69 + Number failed -> 5 Tests failed: + nonstaggered (from CarpetReduce) + staggered (from CarpetReduce) E2xeon_test_dbh (from DistortedBHIVP) Schwarzschild_EF (from Exact) + teukolsky (from WeylScal4) ======================================================================== File [modified]: pelican__2.log Delta lines: +5 -1 =================================================================== --- release-info/pelican__2.log 2010-06-24 14:43:05 UTC (rev 306) +++ release-info/pelican__2.log 2010-06-27 04:34:14 UTC (rev 307) @@ -535,6 +535,7 @@ IOUtil InitBase LAPACK + LORENE LegoExcision LocalInterp LocalReduce @@ -545,6 +546,9 @@ ML_BSSN_Helper ML_BSSN_O2 ML_BSSN_O2_Helper + Meudon_Bin_BH + Meudon_Bin_NS + Meudon_Mag_NS MoL NaNChecker NewRad @@ -759,7 +763,7 @@ Total available tests -> 131 Unrunnable tests -> 48 Runnable tests -> 83 - Total number of thorns -> 120 + Total number of thorns -> 124 Number of tested thorns -> 19 Number of tests passed -> 66 Number passed only to

1 0

[svn:einsteintoolkit] GRHydro/trunk/doc/ (Rev. 128)
by roland.haas＠physics.gatech.edu 26 Jun '10

26 Jun '10

User: rhaas Date: 2010/06/26 08:47 AM Modified: /trunk/doc/ documentation.tex Log: GRHydro: update documentation glitches Thanks to Bruno Giacomazzo for pointing these out. Fix some naming issues with GRHydro vs. Whisky in the GRHydro documentation. The Whisky maintainers are still listed as GRHydro maintainers. File Changes: Directory: /trunk/doc/ ====================== File [modified]: documentation.tex Delta lines: +14 -7 =================================================================== --- trunk/doc/documentation.tex 2010-06-22 01:05:13 UTC (rev 127) +++ trunk/doc/documentation.tex 2010-06-26 13:47:35 UTC (rev 128) @@ -114,9 +114,10 @@ \section{Introduction} \label{sec:intro} -The {\tt GRHydro}\footnote{The name is +The {\tt GRHydro}\footnote{GRHydro started as a branch of Whisky. + The name Whisky is due to Tom Goodale, who pointed out that the in the original Gaelic - {\it uisge beatha} (from which {\tt GRHydro} is derived) meant {\it water of + {\it uisge beatha} (from which {\tt Whisky} is derived) meant {\it water of life}. This name was chosen at a free and fair democratic ballot at the EU Network meeting in Southampton that gave the right answer thanks to a little bit of help from the authors of the code after @@ -135,11 +136,15 @@ you know the required physics and numerical methods, and also the basics of Cactus\footnote{http://www.cactuscode.org}. If you don't, then skip this section and come back to it after reading the rest of this ThornGuide of Cactus. For more details such -as thornlists and parameter files, take a look at the {\tt GRHydro} web page +as thornlists and parameter files, take a look at the Einstein Toolkit web page which is currently stored at \begin{verbatim} -http://www.GRHydrocode.org +http://www.einsteintoolkit.org \end{verbatim} +as well as the original {\tt Whisky} website at +\begin{verbatim} +http://www.whiskycode.org +\end{verbatim} {\tt GRHydro} provides the hydro variables and methods to evolve them. It does not provide any information about initial data or equations of @@ -217,7 +222,9 @@ \subsection{Obtaining This Thorn} The public version of GRHydro can be found on the -website {\tt http://www.GRHydrocode.org}. +website {\tt http://www.einsteintoolkit.org}. +The latest public version of Whisky can be found on the +website {\tt http://www.whiskycode.org}. \subsection{Basic Usage} @@ -293,7 +300,7 @@ The {\tt GRHydro} web page is located at \begin{verbatim} -http://www.GRHydrocode.org +http://www.einsteintoolkit.org \end{verbatim} and contains information on obtaining the code, together with thornlists and sample parameter files. There is also a section on @@ -1974,7 +1981,7 @@ \item ~1998-: Developed as Cactus thorn {\tt MAHC} inside the GR{\tt Astro\_Hydro} arrangement at Washington University, primarily by Mark Miller. -\item 2002-: {\tt GRHydro} written based on {\tt GR3D}. +\item 2002-: {\tt Whisky} written based on {\tt GR3D}. \end{itemize} This is necessarily only a sketch; many people have contributed to

1 0

[svn:einsteintoolkit] TwoPunctures/trunk/src/ (Rev. 110)
by knarf＠cct.lsu.edu 24 Jun '10

24 Jun '10

User: knarf Date: 2010/06/24 01:24 PM Modified: /trunk/src/ TwoPunctures.c Log: remove (leftover) permanent change of the grid functions x y and z File Changes: Directory: /trunk/src/ ====================== File [modified]: TwoPunctures.c Delta lines: +0 -10 =================================================================== --- trunk/src/TwoPunctures.c 2010-06-17 19:59:38 UTC (rev 109) +++ trunk/src/TwoPunctures.c 2010-06-24 18:24:29 UTC (rev 110) @@ -613,16 +613,6 @@ } /* for j */ } /* for k */ - for (int k = imin[2]; k < imax[2]; ++k) - for (int j = imin[1]; j < imax[1]; ++j) - for (int i = imin[0]; i < imax[0]; ++i) - { - int ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); - x[ijk] += center_offset[0]; - y[ijk] += center_offset[1]; - z[ijk] += center_offset[2]; - } - if (use_sources && rescale_sources) { Rescale_Sources(cctkGH,

1 0

[svn:einsteintoolkit] www/toolkit/ (Rev. 306)
by schnetter＠cct.lsu.edu 24 Jun '10

24 Jun '10

User: eschnett Date: 2010/06/24 09:43 AM Modified: /toolkit/ index.php Log: Add pointer to tutorial near list of ET components File Changes: Directory: /toolkit/ ==================== File [modified]: index.php Delta lines: +5 -6 =================================================================== --- toolkit/index.php 2010-06-24 14:40:01 UTC (rev 305) +++ toolkit/index.php 2010-06-24 14:43:05 UTC (rev 306) @@ -28,14 +28,13 @@ with partnerships with other developers who contribute open software and coordinate together on development. - The components of the Einstein Toolkit include:

1 0

[svn:einsteintoolkit] www/toolkit/ (Rev. 305)
by schnetter＠cct.lsu.edu 24 Jun '10

24 Jun '10

User: eschnett Date: 2010/06/24 09:40 AM Modified: /toolkit/ index.php Log: Correct broken links File Changes: Directory: /toolkit/ ==================== File [modified]: index.php Delta lines: +8 -6 =================================================================== --- toolkit/index.php 2010-06-24 14:33:34 UTC (rev 304) +++ toolkit/index.php 2010-06-24 14:40:01 UTC (rev 305) @@ -34,7 +34,8 @@ thornlist: <a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…">Einstein Toolkit Thornlist</a>. A <a href="../info/documentation/">worked example</a> describes how to download, compile and run the - toolkit. --> + toolkit. +--> The components of the Einstein Toolkit include: @@ -64,11 +65,10 @@ Retrieval Language and can be checked out using the <a href="https://svn.cactuscode.org/Utilities/branches/ET_2010_06/Scripts/GetCompone…" class="nostyle">GetComponents</a> tools. -Distribution: <a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…">https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…</a> +Distribution: <a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…">https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…</a> <a href="https://svn.cactuscode.org/Utilities/branches/ET_2010_06/Scripts/GetCompone…">https://svn.cactuscode.org/Utilities/branches/ET_2010_06/Scripts/GetCompone…</a> </td> - </tr> <tr> @@ -84,7 +84,7 @@ resources. Distribution: CCT SVN -repository: <a href="https://svn.cct.lsu.edu/repos/numrel/simfactory">https://svn.cct.lsu.edu/repos/numrel/simfactory</a></td> +repository: <a href="https://svn.cct.lsu.edu/repos/numrel/simfactory/branches/ET_2010_06">https://svn.cct.lsu.edu/repos/numrel/simfactory/branches/ET_2010_06</a></td> </tr> <tr> @@ -93,9 +93,11 @@ Toolkit. Additional thorns are maintained in external repositories with open access. Distribution: Einstein Toolkit -SVN <a href="https://svn.einsteintoolkit.org/cactus">https://svn.einsteintoolkit.org/cactus/</a> +SVN repository overview +<a href="https://svn.einsteintoolkit.org/cactus">https://svn.einsteintoolkit.org/cactus/</a> </td> </tr> + <tr> <td valign="top">Parameter Files</td> <td>Example Cactus code parameter files for Einstein Toolkit @@ -104,7 +106,7 @@ the documentation for information on how to run these on the QueenBee machine. -Distribution: <a href="https://svn.einsteintoolkit.org/cactus/EinsteinExamples/par">https://svn.einsteintoolkit.org/cactus/EinsteinExamples/par</a> +Distribution: <a href="https://svn.einsteintoolkit.org/cactus/EinsteinExamples/branches/ET_2010_06…">https://svn.einsteintoolkit.org/cactus/EinsteinExamples/branches/ET_2010_06…</a> </td> </tr> </table>

1 0

[svn:einsteintoolkit] www/toolkit/ (Rev. 304)
by schnetter＠cct.lsu.edu 24 Jun '10

24 Jun '10

User: eschnett Date: 2010/06/24 09:33 AM Modified: /toolkit/ index.php Log: Correct url File Changes: Directory: /toolkit/ ==================== File [modified]: index.php Delta lines: +1 -1 =================================================================== --- toolkit/index.php 2010-06-23 15:24:07 UTC (rev 303) +++ toolkit/index.php 2010-06-24 14:33:34 UTC (rev 304) @@ -64,7 +64,7 @@ Retrieval Language and can be checked out using the <a href="https://svn.cactuscode.org/Utilities/branches/ET_2010_06/Scripts/GetCompone…" class="nostyle">GetComponents</a> tools. -Distribution: <a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…">https://svn.einsteintoolkit.org/manifest/einsteintoolkit.th</a> +Distribution: <a href="https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…">https://svn.einsteintoolkit.org/manifest/branches/ET_2010_06/einsteintoolki…</a> <a href="https://svn.cactuscode.org/Utilities/branches/ET_2010_06/Scripts/GetCompone…">https://svn.cactuscode.org/Utilities/branches/ET_2010_06/Scripts/GetCompone…</a> </td>

1 0

[svn:einsteintoolkit] WeylScal4/trunk/test/ (Rev. 57)
by roland.haas＠physics.gatech.edu 23 Jun '10

23 Jun '10

User: rhaas Date: 2010/06/23 10:58 AM Modified: /trunk/test/ teukolskyID.par /trunk/test/teukolskyID/ Psi4i.d.asc, Psi4i.x.asc, Psi4i.y.asc, Psi4i.z.asc, Psi4r.d.asc, Psi4r.x.asc, Psi4r.y.asc, Psi4r.z.asc Log: WeylScal4: fix teukolksyID test suite test suite was non-functional before since I had forgotton to add MoL to the ActiveThorns list File Changes: Directory: /trunk/test/ ======================= File [modified]: teukolskyID.par Delta lines: +1 -1 =================================================================== --- trunk/test/teukolskyID.par 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID.par 2010-06-23 15:58:45 UTC (rev 57) @@ -4,7 +4,7 @@ # $Header$ # Required thorns -ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask LocalReduce IDLinearWaves LoopControl GSL WeylScal4 ReflectionSymmetry" +ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask LocalReduce IDLinearWaves LoopControl GSL WeylScal4 ReflectionSymmetry MoL" # Grid driver::ghost_size = 1 Directory: /trunk/test/teukolskyID/ =================================== File [modified]: Psi4i.d.asc Delta lines: +11 -11 =================================================================== --- trunk/test/teukolskyID/Psi4i.d.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4i.d.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,11 +1,11 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 0 0 0 0 -2.0625 -2.0625 -0.5 -2.8396262443943e+238 -0 0 0 0 0 4 4 4 0 -1.5625 -1.5625 0 2.8396262443943e+238 -0 0 0 0 0 8 8 8 0 -1.0625 -1.0625 0.5 2.8396262443943e+238 -0 0 0 0 0 12 12 12 0 -0.5625 -0.5625 1 2.8396262443943e+238 -0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 2.8396262443943e+238 -0 0 0 0 0 20 20 20 0 0.4375 0.4375 2 2.8396262443943e+238 +0 0 0 0 0 0 0 0 0 -2.0625 -2.0625 -0.5 -0.000947609598144528 +0 0 0 0 0 4 4 4 0 -1.5625 -1.5625 0 -0 +0 0 0 0 0 8 8 8 0 -1.0625 -1.0625 0.5 0.000356282875172737 +0 0 0 0 0 12 12 12 0 -0.5625 -0.5625 1 -0.0263498098998352 +0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 -0.00904684387067402 +0 0 0 0 0 20 20 20 0 0.4375 0.4375 2 -0.00119929445931701 # # # @@ -14,11 +14,11 @@ # # # -0 0 1 0 0 6 6 6 0 -1.3125 -1.3125 0.25 2.8396262443943e+238 -0 0 1 0 0 8 8 8 0 -1.0625 -1.0625 0.5 2.8396262443943e+238 -0 0 1 0 0 10 10 10 0 -0.8125 -0.8125 0.75 2.8396262443943e+238 -0 0 1 0 0 12 12 12 0 -0.5625 -0.5625 1 2.8396262443943e+238 -0 0 1 0 0 14 14 14 0 -0.3125 -0.3125 1.25 2.8396262443943e+238 +0 0 1 0 0 6 6 6 0 -1.3125 -1.3125 0.25 0.000712240897234611 +0 0 1 0 0 8 8 8 0 -1.0625 -1.0625 0.5 0.000356282875172737 +0 0 1 0 0 10 10 10 0 -0.8125 -0.8125 0.75 -0.0112792424499431 +0 0 1 0 0 12 12 12 0 -0.5625 -0.5625 1 -0.0263498098998352 +0 0 1 0 0 14 14 14 0 -0.3125 -0.3125 1.25 -0.0143995291686149 # # # File [modified]: Psi4i.x.asc Delta lines: +31 -31 =================================================================== --- trunk/test/teukolskyID/Psi4i.x.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4i.x.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,39 +1,39 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 0 16 4 0 -2.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 4 16 4 0 -1.5625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 8 16 4 0 -1.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 12 16 4 0 -0.5625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 20 16 4 0 0.4375 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 24 16 4 0 0.9375 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 28 16 4 0 1.4375 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 32 16 4 0 1.9375 -0.0625 0 2.8396262443943e+238 +0 0 0 0 0 0 16 4 0 -2.0625 -0.0625 0 -0 +0 0 0 0 0 4 16 4 0 -1.5625 -0.0625 0 -0 +0 0 0 0 0 8 16 4 0 -1.0625 -0.0625 0 -0 +0 0 0 0 0 12 16 4 0 -0.5625 -0.0625 0 0 +0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 0 +0 0 0 0 0 20 16 4 0 0.4375 -0.0625 0 0 +0 0 0 0 0 24 16 4 0 0.9375 -0.0625 0 0 +0 0 0 0 0 28 16 4 0 1.4375 -0.0625 0 0 +0 0 0 0 0 32 16 4 0 1.9375 -0.0625 0 0 -0 0 1 0 0 6 18 4 0 -1.3125 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 8 18 4 0 -1.0625 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 10 18 4 0 -0.8125 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 12 18 4 0 -0.5625 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 14 18 4 0 -0.3125 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 16 18 4 0 -0.0625 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 20 18 4 0 0.4375 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 22 18 4 0 0.6875 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 24 18 4 0 0.9375 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 26 18 4 0 1.1875 0.1875 0 2.8396262443943e+238 +0 0 1 0 0 6 18 4 0 -1.3125 0.1875 0 0 +0 0 1 0 0 8 18 4 0 -1.0625 0.1875 0 0 +0 0 1 0 0 10 18 4 0 -0.8125 0.1875 0 0 +0 0 1 0 0 12 18 4 0 -0.5625 0.1875 0 0 +0 0 1 0 0 14 18 4 0 -0.3125 0.1875 0 0 +0 0 1 0 0 16 18 4 0 -0.0625 0.1875 0 0 +0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 0 +0 0 1 0 0 20 18 4 0 0.4375 0.1875 0 0 +0 0 1 0 0 22 18 4 0 0.6875 0.1875 0 0 +0 0 1 0 0 24 18 4 0 0.9375 0.1875 0 0 +0 0 1 0 0 26 18 4 0 1.1875 0.1875 0 0 -0 0 2 0 0 11 17 4 0 -0.6875 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 12 17 4 0 -0.5625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 13 17 4 0 -0.4375 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 14 17 4 0 -0.3125 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 15 17 4 0 -0.1875 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 16 17 4 0 -0.0625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 18 17 4 0 0.1875 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 19 17 4 0 0.3125 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 20 17 4 0 0.4375 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 21 17 4 0 0.5625 0.0625 0 2.8396262443943e+238 +0 0 2 0 0 11 17 4 0 -0.6875 0.0625 0 0 +0 0 2 0 0 12 17 4 0 -0.5625 0.0625 0 0 +0 0 2 0 0 13 17 4 0 -0.4375 0.0625 0 0 +0 0 2 0 0 14 17 4 0 -0.3125 0.0625 0 0 +0 0 2 0 0 15 17 4 0 -0.1875 0.0625 0 0 +0 0 2 0 0 16 17 4 0 -0.0625 0.0625 0 0 +0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 0 +0 0 2 0 0 18 17 4 0 0.1875 0.0625 0 0 +0 0 2 0 0 19 17 4 0 0.3125 0.0625 0 0 +0 0 2 0 0 20 17 4 0 0.4375 0.0625 0 0 +0 0 2 0 0 21 17 4 0 0.5625 0.0625 0 0 File [modified]: Psi4i.y.asc Delta lines: +31 -31 =================================================================== --- trunk/test/teukolskyID/Psi4i.y.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4i.y.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,39 +1,39 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 16 0 4 0 -0.0625 -2.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 4 4 0 -0.0625 -1.5625 0 2.8396262443943e+238 -0 0 0 0 0 16 8 4 0 -0.0625 -1.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 12 4 0 -0.0625 -0.5625 0 2.8396262443943e+238 -0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 20 4 0 -0.0625 0.4375 0 2.8396262443943e+238 -0 0 0 0 0 16 24 4 0 -0.0625 0.9375 0 2.8396262443943e+238 -0 0 0 0 0 16 28 4 0 -0.0625 1.4375 0 2.8396262443943e+238 -0 0 0 0 0 16 32 4 0 -0.0625 1.9375 0 2.8396262443943e+238 +0 0 0 0 0 16 0 4 0 -0.0625 -2.0625 0 -0 +0 0 0 0 0 16 4 4 0 -0.0625 -1.5625 0 -0 +0 0 0 0 0 16 8 4 0 -0.0625 -1.0625 0 -0 +0 0 0 0 0 16 12 4 0 -0.0625 -0.5625 0 0 +0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 0 +0 0 0 0 0 16 20 4 0 -0.0625 0.4375 0 0 +0 0 0 0 0 16 24 4 0 -0.0625 0.9375 0 0 +0 0 0 0 0 16 28 4 0 -0.0625 1.4375 0 0 +0 0 0 0 0 16 32 4 0 -0.0625 1.9375 0 0 -0 0 1 0 0 18 6 4 0 0.1875 -1.3125 0 2.8396262443943e+238 -0 0 1 0 0 18 8 4 0 0.1875 -1.0625 0 2.8396262443943e+238 -0 0 1 0 0 18 10 4 0 0.1875 -0.8125 0 2.8396262443943e+238 -0 0 1 0 0 18 12 4 0 0.1875 -0.5625 0 2.8396262443943e+238 -0 0 1 0 0 18 14 4 0 0.1875 -0.3125 0 2.8396262443943e+238 -0 0 1 0 0 18 16 4 0 0.1875 -0.0625 0 2.8396262443943e+238 -0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 18 20 4 0 0.1875 0.4375 0 2.8396262443943e+238 -0 0 1 0 0 18 22 4 0 0.1875 0.6875 0 2.8396262443943e+238 -0 0 1 0 0 18 24 4 0 0.1875 0.9375 0 2.8396262443943e+238 -0 0 1 0 0 18 26 4 0 0.1875 1.1875 0 2.8396262443943e+238 +0 0 1 0 0 18 6 4 0 0.1875 -1.3125 0 0 +0 0 1 0 0 18 8 4 0 0.1875 -1.0625 0 0 +0 0 1 0 0 18 10 4 0 0.1875 -0.8125 0 0 +0 0 1 0 0 18 12 4 0 0.1875 -0.5625 0 0 +0 0 1 0 0 18 14 4 0 0.1875 -0.3125 0 0 +0 0 1 0 0 18 16 4 0 0.1875 -0.0625 0 0 +0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 0 +0 0 1 0 0 18 20 4 0 0.1875 0.4375 0 0 +0 0 1 0 0 18 22 4 0 0.1875 0.6875 0 0 +0 0 1 0 0 18 24 4 0 0.1875 0.9375 0 0 +0 0 1 0 0 18 26 4 0 0.1875 1.1875 0 0 -0 0 2 0 0 17 11 4 0 0.0625 -0.6875 0 2.8396262443943e+238 -0 0 2 0 0 17 12 4 0 0.0625 -0.5625 0 2.8396262443943e+238 -0 0 2 0 0 17 13 4 0 0.0625 -0.4375 0 2.8396262443943e+238 -0 0 2 0 0 17 14 4 0 0.0625 -0.3125 0 2.8396262443943e+238 -0 0 2 0 0 17 15 4 0 0.0625 -0.1875 0 2.8396262443943e+238 -0 0 2 0 0 17 16 4 0 0.0625 -0.0625 0 2.8396262443943e+238 -0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 17 18 4 0 0.0625 0.1875 0 2.8396262443943e+238 -0 0 2 0 0 17 19 4 0 0.0625 0.3125 0 2.8396262443943e+238 -0 0 2 0 0 17 20 4 0 0.0625 0.4375 0 2.8396262443943e+238 -0 0 2 0 0 17 21 4 0 0.0625 0.5625 0 2.8396262443943e+238 +0 0 2 0 0 17 11 4 0 0.0625 -0.6875 0 0 +0 0 2 0 0 17 12 4 0 0.0625 -0.5625 0 0 +0 0 2 0 0 17 13 4 0 0.0625 -0.4375 0 0 +0 0 2 0 0 17 14 4 0 0.0625 -0.3125 0 0 +0 0 2 0 0 17 15 4 0 0.0625 -0.1875 0 0 +0 0 2 0 0 17 16 4 0 0.0625 -0.0625 0 0 +0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 0 +0 0 2 0 0 17 18 4 0 0.0625 0.1875 0 0 +0 0 2 0 0 17 19 4 0 0.0625 0.3125 0 0 +0 0 2 0 0 17 20 4 0 0.0625 0.4375 0 0 +0 0 2 0 0 17 21 4 0 0.0625 0.5625 0 0 File [modified]: Psi4i.z.asc Delta lines: +20 -20 =================================================================== --- trunk/test/teukolskyID/Psi4i.z.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4i.z.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,28 +1,28 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 16 16 0 0 -0.0625 -0.0625 -0.5 -2.8396262443943e+238 -0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 16 8 0 -0.0625 -0.0625 0.5 2.8396262443943e+238 -0 0 0 0 0 16 16 12 0 -0.0625 -0.0625 1 2.8396262443943e+238 -0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 2.8396262443943e+238 -0 0 0 0 0 16 16 20 0 -0.0625 -0.0625 2 2.8396262443943e+238 +0 0 0 0 0 16 16 0 0 -0.0625 -0.0625 -0.5 -0.0204252046237925 +0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 0 +0 0 0 0 0 16 16 8 0 -0.0625 -0.0625 0.5 0.0204252046237925 +0 0 0 0 0 16 16 12 0 -0.0625 -0.0625 1 -0.0617664380239765 +0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 -0.00904684387067402 +0 0 0 0 0 16 16 20 0 -0.0625 -0.0625 2 -0.00904684387067402 -0 0 1 0 0 16 16 2 0 -0.0625 -0.0625 -0.25 -2.8396262443943e+238 -0 0 1 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 1 0 0 16 16 6 0 -0.0625 -0.0625 0.25 2.8396262443943e+238 -0 0 1 0 0 16 16 8 0 -0.0625 -0.0625 0.5 2.8396262443943e+238 -0 0 1 0 0 16 16 10 0 -0.0625 -0.0625 0.75 2.8396262443943e+238 -0 0 1 0 0 16 16 12 0 -0.0625 -0.0625 1 2.8396262443943e+238 -0 0 1 0 0 16 16 14 0 -0.0625 -0.0625 1.25 2.8396262443943e+238 +0 0 1 0 0 16 16 2 0 -0.0625 -0.0625 -0.25 -0.108107083030808 +0 0 1 0 0 16 16 4 0 -0.0625 -0.0625 0 0 +0 0 1 0 0 16 16 6 0 -0.0625 -0.0625 0.25 0.108107083030808 +0 0 1 0 0 16 16 8 0 -0.0625 -0.0625 0.5 0.0204252046237925 +0 0 1 0 0 16 16 10 0 -0.0625 -0.0625 0.75 -0.0472297602104879 +0 0 1 0 0 16 16 12 0 -0.0625 -0.0625 1 -0.0617664380239765 +0 0 1 0 0 16 16 14 0 -0.0625 -0.0625 1.25 -0.0354066409473252 -0 0 2 0 0 16 16 3 0 -0.0625 -0.0625 -0.125 2.8396262443943e+238 -0 0 2 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 2 0 0 16 16 5 0 -0.0625 -0.0625 0.125 2.8396262443943e+238 -0 0 2 0 0 16 16 6 0 -0.0625 -0.0625 0.25 2.8396262443943e+238 -0 0 2 0 0 16 16 7 0 -0.0625 -0.0625 0.375 2.8396262443943e+238 -0 0 2 0 0 16 16 8 0 -0.0625 -0.0625 0.5 2.8396262443943e+238 -0 0 2 0 0 16 16 9 0 -0.0625 -0.0625 0.625 2.8396262443943e+238 +0 0 2 0 0 16 16 3 0 -0.0625 -0.0625 -0.125 -0.119774776965786 +0 0 2 0 0 16 16 4 0 -0.0625 -0.0625 0 0 +0 0 2 0 0 16 16 5 0 -0.0625 -0.0625 0.125 0.119774776965786 +0 0 2 0 0 16 16 6 0 -0.0625 -0.0625 0.25 0.108107083030808 +0 0 2 0 0 16 16 7 0 -0.0625 -0.0625 0.375 0.0676017421667073 +0 0 2 0 0 16 16 8 0 -0.0625 -0.0625 0.5 0.0204252046237925 +0 0 2 0 0 16 16 9 0 -0.0625 -0.0625 0.625 -0.0133627027966291 File [modified]: Psi4r.d.asc Delta lines: +11 -11 =================================================================== --- trunk/test/teukolskyID/Psi4r.d.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4r.d.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,11 +1,11 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 0 0 0 0 -2.0625 -2.0625 -0.5 2.8396262443943e+238 -0 0 0 0 0 4 4 4 0 -1.5625 -1.5625 0 2.8396262443943e+238 -0 0 0 0 0 8 8 8 0 -1.0625 -1.0625 0.5 2.8396262443943e+238 -0 0 0 0 0 12 12 12 0 -0.5625 -0.5625 1 2.8396262443943e+238 -0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 2.8396262443943e+238 -0 0 0 0 0 20 20 20 0 0.4375 0.4375 2 2.8396262443943e+238 +0 0 0 0 0 0 0 0 0 -2.0625 -2.0625 -0.5 -2.37783074190988e-07 +0 0 0 0 0 4 4 4 0 -1.5625 -1.5625 0 -2.29325329647439e-07 +0 0 0 0 0 8 8 8 0 -1.0625 -1.0625 0.5 -9.43217025312457e-06 +0 0 0 0 0 12 12 12 0 -0.5625 -0.5625 1 -0.000134686367417486 +0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 -3.65666746534011e-05 +0 0 0 0 0 20 20 20 0 0.4375 0.4375 2 3.14617153825462e-07 # # # @@ -14,11 +14,11 @@ # # # -0 0 1 0 0 6 6 6 0 -1.3125 -1.3125 0.25 2.8396262443943e+238 -0 0 1 0 0 8 8 8 0 -1.0625 -1.0625 0.5 2.8396262443943e+238 -0 0 1 0 0 10 10 10 0 -0.8125 -0.8125 0.75 2.8396262443943e+238 -0 0 1 0 0 12 12 12 0 -0.5625 -0.5625 1 2.8396262443943e+238 -0 0 1 0 0 14 14 14 0 -0.3125 -0.3125 1.25 2.8396262443943e+238 +0 0 1 0 0 6 6 6 0 -1.3125 -1.3125 0.25 -3.54957385346084e-06 +0 0 1 0 0 8 8 8 0 -1.0625 -1.0625 0.5 -9.43217025312457e-06 +0 0 1 0 0 10 10 10 0 -0.8125 -0.8125 0.75 -3.16134062867657e-05 +0 0 1 0 0 12 12 12 0 -0.5625 -0.5625 1 -0.000134686367417486 +0 0 1 0 0 14 14 14 0 -0.3125 -0.3125 1.25 -8.12440763420316e-05 # # # File [modified]: Psi4r.x.asc Delta lines: +31 -31 =================================================================== --- trunk/test/teukolskyID/Psi4r.x.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4r.x.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,39 +1,39 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 0 16 4 0 -2.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 4 16 4 0 -1.5625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 8 16 4 0 -1.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 12 16 4 0 -0.5625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 20 16 4 0 0.4375 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 24 16 4 0 0.9375 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 28 16 4 0 1.4375 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 32 16 4 0 1.9375 -0.0625 0 2.8396262443943e+238 +0 0 0 0 0 0 16 4 0 -2.0625 -0.0625 0 -0.00223634388343174 +0 0 0 0 0 4 16 4 0 -1.5625 -0.0625 0 -0.00223634388343174 +0 0 0 0 0 8 16 4 0 -1.0625 -0.0625 0 -0.0287052634203998 +0 0 0 0 0 12 16 4 0 -0.5625 -0.0625 0 -0.000402190092767988 +0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 -0.00380414718797543 +0 0 0 0 0 20 16 4 0 0.4375 -0.0625 0 0.0218324172625981 +0 0 0 0 0 24 16 4 0 0.9375 -0.0625 0 -0.0310833376801553 +0 0 0 0 0 28 16 4 0 1.4375 -0.0625 0 -0.00719566764383696 +0 0 0 0 0 32 16 4 0 1.9375 -0.0625 0 -0.00719566764383696 -0 0 1 0 0 6 18 4 0 -1.3125 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 8 18 4 0 -1.0625 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 10 18 4 0 -0.8125 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 12 18 4 0 -0.5625 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 14 18 4 0 -0.3125 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 16 18 4 0 -0.0625 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 20 18 4 0 0.4375 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 22 18 4 0 0.6875 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 24 18 4 0 0.9375 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 26 18 4 0 1.1875 0.1875 0 2.8396262443943e+238 +0 0 1 0 0 6 18 4 0 -1.3125 0.1875 0 -0.012242836303987 +0 0 1 0 0 8 18 4 0 -1.0625 0.1875 0 -0.0266155200058717 +0 0 1 0 0 10 18 4 0 -0.8125 0.1875 0 -0.0258730223353572 +0 0 1 0 0 12 18 4 0 -0.5625 0.1875 0 -0.00413899377540188 +0 0 1 0 0 14 18 4 0 -0.3125 0.1875 0 0.0164710873760525 +0 0 1 0 0 16 18 4 0 -0.0625 0.1875 0 -0.0564077253562977 +0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 -0.00284240990537506 +0 0 1 0 0 20 18 4 0 0.4375 0.1875 0 0.0105843457689293 +0 0 1 0 0 22 18 4 0 0.6875 0.1875 0 -0.0168325588059544 +0 0 1 0 0 24 18 4 0 0.9375 0.1875 0 -0.0289532241479415 +0 0 1 0 0 26 18 4 0 1.1875 0.1875 0 -0.0153415732468004 -0 0 2 0 0 11 17 4 0 -0.6875 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 12 17 4 0 -0.5625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 13 17 4 0 -0.4375 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 14 17 4 0 -0.3125 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 15 17 4 0 -0.1875 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 16 17 4 0 -0.0625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 18 17 4 0 0.1875 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 19 17 4 0 0.3125 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 20 17 4 0 0.4375 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 21 17 4 0 0.5625 0.0625 0 2.8396262443943e+238 +0 0 2 0 0 11 17 4 0 -0.6875 0.0625 0 -0.0140780471998292 +0 0 2 0 0 12 17 4 0 -0.5625 0.0625 0 -0.000402190092767988 +0 0 2 0 0 13 17 4 0 -0.4375 0.0625 0 0.021832417262598 +0 0 2 0 0 14 17 4 0 -0.3125 0.0625 0 0.0427783767090442 +0 0 2 0 0 15 17 4 0 -0.1875 0.0625 0 0.0520239956084546 +0 0 2 0 0 16 17 4 0 -0.0625 0.0625 0 -0.00380414718797543 +0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 -0.00380414718797369 +0 0 2 0 0 18 17 4 0 0.1875 0.0625 0 0.0520239956084533 +0 0 2 0 0 19 17 4 0 0.3125 0.0625 0 0.0427783767090442 +0 0 2 0 0 20 17 4 0 0.4375 0.0625 0 0.0218324172625981 +0 0 2 0 0 21 17 4 0 0.5625 0.0625 0 -0.000572015851518783 File [modified]: Psi4r.y.asc Delta lines: +31 -31 =================================================================== --- trunk/test/teukolskyID/Psi4r.y.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4r.y.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,39 +1,39 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 16 0 4 0 -0.0625 -2.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 4 4 0 -0.0625 -1.5625 0 2.8396262443943e+238 -0 0 0 0 0 16 8 4 0 -0.0625 -1.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 12 4 0 -0.0625 -0.5625 0 2.8396262443943e+238 -0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 20 4 0 -0.0625 0.4375 0 2.8396262443943e+238 -0 0 0 0 0 16 24 4 0 -0.0625 0.9375 0 2.8396262443943e+238 -0 0 0 0 0 16 28 4 0 -0.0625 1.4375 0 2.8396262443943e+238 -0 0 0 0 0 16 32 4 0 -0.0625 1.9375 0 2.8396262443943e+238 +0 0 0 0 0 16 0 4 0 -0.0625 -2.0625 0 0.00225517133772647 +0 0 0 0 0 16 4 4 0 -0.0625 -1.5625 0 0.00225517133772647 +0 0 0 0 0 16 8 4 0 -0.0625 -1.0625 0 0.029169513176083 +0 0 0 0 0 16 12 4 0 -0.0625 -0.5625 0 0.000308978078890513 +0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 -0.00380414718797543 +0 0 0 0 0 16 20 4 0 -0.0625 0.4375 0 -0.0229517575833702 +0 0 0 0 0 16 24 4 0 -0.0625 0.9375 0 0.0317138081444724 +0 0 0 0 0 16 28 4 0 -0.0625 1.4375 0 0.00725546291128947 +0 0 0 0 0 16 32 4 0 -0.0625 1.9375 0 0.00725546291128947 -0 0 1 0 0 18 6 4 0 0.1875 -1.3125 0 2.8396262443943e+238 -0 0 1 0 0 18 8 4 0 0.1875 -1.0625 0 2.8396262443943e+238 -0 0 1 0 0 18 10 4 0 0.1875 -0.8125 0 2.8396262443943e+238 -0 0 1 0 0 18 12 4 0 0.1875 -0.5625 0 2.8396262443943e+238 -0 0 1 0 0 18 14 4 0 0.1875 -0.3125 0 2.8396262443943e+238 -0 0 1 0 0 18 16 4 0 0.1875 -0.0625 0 2.8396262443943e+238 -0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 2.8396262443943e+238 -0 0 1 0 0 18 20 4 0 0.1875 0.4375 0 2.8396262443943e+238 -0 0 1 0 0 18 22 4 0 0.1875 0.6875 0 2.8396262443943e+238 -0 0 1 0 0 18 24 4 0 0.1875 0.9375 0 2.8396262443943e+238 -0 0 1 0 0 18 26 4 0 0.1875 1.1875 0 2.8396262443943e+238 +0 0 1 0 0 18 6 4 0 0.1875 -1.3125 0 0.0124288119092335 +0 0 1 0 0 18 8 4 0 0.1875 -1.0625 0 0.0270366975580574 +0 0 1 0 0 18 10 4 0 0.1875 -0.8125 0 0.0264742495681141 +0 0 1 0 0 18 12 4 0 0.1875 -0.5625 0 0.00382959303505998 +0 0 1 0 0 18 14 4 0 0.1875 -0.3125 0 -0.020087727325365 +0 0 1 0 0 18 16 4 0 0.1875 -0.0625 0 0.0520239956084533 +0 0 1 0 0 18 18 4 0 0.1875 0.1875 0 -0.00284240990537506 +0 0 1 0 0 18 20 4 0 0.1875 0.4375 0 -0.0122164951298324 +0 0 1 0 0 18 22 4 0 0.1875 0.6875 0 0.0171867546851686 +0 0 1 0 0 18 24 4 0 0.1875 0.9375 0 0.0295264979402503 +0 0 1 0 0 18 26 4 0 0.1875 1.1875 0 0.0156133010615632 -0 0 2 0 0 17 11 4 0 0.0625 -0.6875 0 2.8396262443943e+238 -0 0 2 0 0 17 12 4 0 0.0625 -0.5625 0 2.8396262443943e+238 -0 0 2 0 0 17 13 4 0 0.0625 -0.4375 0 2.8396262443943e+238 -0 0 2 0 0 17 14 4 0 0.0625 -0.3125 0 2.8396262443943e+238 -0 0 2 0 0 17 15 4 0 0.0625 -0.1875 0 2.8396262443943e+238 -0 0 2 0 0 17 16 4 0 0.0625 -0.0625 0 2.8396262443943e+238 -0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 2.8396262443943e+238 -0 0 2 0 0 17 18 4 0 0.0625 0.1875 0 2.8396262443943e+238 -0 0 2 0 0 17 19 4 0 0.0625 0.3125 0 2.8396262443943e+238 -0 0 2 0 0 17 20 4 0 0.0625 0.4375 0 2.8396262443943e+238 -0 0 2 0 0 17 21 4 0 0.0625 0.5625 0 2.8396262443943e+238 +0 0 2 0 0 17 11 4 0 0.0625 -0.6875 0 0.0142879718826407 +0 0 2 0 0 17 12 4 0 0.0625 -0.5625 0 0.000308978078890569 +0 0 2 0 0 17 13 4 0 0.0625 -0.4375 0 -0.0229517575833701 +0 0 2 0 0 17 14 4 0 0.0625 -0.3125 0 -0.0452873223411553 +0 0 2 0 0 17 15 4 0 0.0625 -0.1875 0 -0.0564077253562966 +0 0 2 0 0 17 16 4 0 0.0625 -0.0625 0 -0.00380414718797324 +0 0 2 0 0 17 17 4 0 0.0625 0.0625 0 -0.00380414718797369 +0 0 2 0 0 17 18 4 0 0.0625 0.1875 0 -0.0564077253562969 +0 0 2 0 0 17 19 4 0 0.0625 0.3125 0 -0.0452873223411553 +0 0 2 0 0 17 20 4 0 0.0625 0.4375 0 -0.0229517575833701 +0 0 2 0 0 17 21 4 0 0.0625 0.5625 0 0.00011624029177533 File [modified]: Psi4r.z.asc Delta lines: +20 -20 =================================================================== --- trunk/test/teukolskyID/Psi4r.z.asc 2010-06-23 15:52:15 UTC (rev 56) +++ trunk/test/teukolskyID/Psi4r.z.asc 2010-06-23 15:58:45 UTC (rev 57) @@ -1,28 +1,28 @@ # 1D ASCII output created by CarpetIOASCII # -0 0 0 0 0 16 16 0 0 -0.0625 -0.0625 -0.5 2.8396262443943e+238 -0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 0 0 0 16 16 8 0 -0.0625 -0.0625 0.5 2.8396262443943e+238 -0 0 0 0 0 16 16 12 0 -0.0625 -0.0625 1 2.8396262443943e+238 -0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 2.8396262443943e+238 -0 0 0 0 0 16 16 20 0 -0.0625 -0.0625 2 2.8396262443943e+238 +0 0 0 0 0 16 16 0 0 -0.0625 -0.0625 -0.5 -0.000164374418623651 +0 0 0 0 0 16 16 4 0 -0.0625 -0.0625 0 -0.00380414718797543 +0 0 0 0 0 16 16 8 0 -0.0625 -0.0625 0.5 -0.000164374418623651 +0 0 0 0 0 16 16 12 0 -0.0625 -0.0625 1 -0.000543113416260428 +0 0 0 0 0 16 16 16 0 -0.0625 -0.0625 1.5 -3.65666746534011e-05 +0 0 0 0 0 16 16 20 0 -0.0625 -0.0625 2 -3.65666746534011e-05 -0 0 1 0 0 16 16 2 0 -0.0625 -0.0625 -0.25 2.8396262443943e+238 -0 0 1 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 1 0 0 16 16 6 0 -0.0625 -0.0625 0.25 2.8396262443943e+238 -0 0 1 0 0 16 16 8 0 -0.0625 -0.0625 0.5 2.8396262443943e+238 -0 0 1 0 0 16 16 10 0 -0.0625 -0.0625 0.75 2.8396262443943e+238 -0 0 1 0 0 16 16 12 0 -0.0625 -0.0625 1 2.8396262443943e+238 -0 0 1 0 0 16 16 14 0 -0.0625 -0.0625 1.25 2.8396262443943e+238 +0 0 1 0 0 16 16 2 0 -0.0625 -0.0625 -0.25 -0.00217586688787084 +0 0 1 0 0 16 16 4 0 -0.0625 -0.0625 0 -0.00380414718797543 +0 0 1 0 0 16 16 6 0 -0.0625 -0.0625 0.25 -0.00217586688787084 +0 0 1 0 0 16 16 8 0 -0.0625 -0.0625 0.5 -0.000164374418623651 +0 0 1 0 0 16 16 10 0 -0.0625 -0.0625 0.75 -0.000193220264071006 +0 0 1 0 0 16 16 12 0 -0.0625 -0.0625 1 -0.000543113416260428 +0 0 1 0 0 16 16 14 0 -0.0625 -0.0625 1.25 -0.000289840045456915 -0 0 2 0 0 16 16 3 0 -0.0625 -0.0625 -0.125 2.8396262443943e+238 -0 0 2 0 0 16 16 4 0 -0.0625 -0.0625 0 2.8396262443943e+238 -0 0 2 0 0 16 16 5 0 -0.0625 -0.0625 0.125 2.8396262443943e+238 -0 0 2 0 0 16 16 6 0 -0.0625 -0.0625 0.25 2.8396262443943e+238 -0 0 2 0 0 16 16 7 0 -0.0625 -0.0625 0.375 2.8396262443943e+238 -0 0 2 0 0 16 16 8 0 -0.0625 -0.0625 0.5 2.8396262443943e+238 -0 0 2 0 0 16 16 9 0 -0.0625 -0.0625 0.625 2.8396262443943e+238 +0 0 2 0 0 16 16 3 0 -0.0625 -0.0625 -0.125 -0.00333097979395648 +0 0 2 0 0 16 16 4 0 -0.0625 -0.0625 0 -0.00380414718797543 +0 0 2 0 0 16 16 5 0 -0.0625 -0.0625 0.125 -0.00333097979395648 +0 0 2 0 0 16 16 6 0 -0.0625 -0.0625 0.25 -0.00217586688787084 +0 0 2 0 0 16 16 7 0 -0.0625 -0.0625 0.375 -0.000943554132152407 +0 0 2 0 0 16 16 8 0 -0.0625 -0.0625 0.5 -0.000164374418623651 +0 0 2 0 0 16 16 9 0 -0.0625 -0.0625 0.625 -0.000178384667179087

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits June 2010