User: rhaas Date: 2010/07/19 08:46 AM

Added: /branches/momentum_constraint/test/ bhns_momentum_fd.par /branches/momentum_constraint/test/bhns_momentum_fd/ alp.d.asc, alp.maximum.asc, alp.minimum.asc, alp.norm1.asc, alp.norm2.asc, alp.x.asc, alp.y.asc, alp.z.asc, betax.d.asc, betax.maximum.asc, betax.minimum.asc, betax.norm1.asc, betax.norm2.asc, betax.x.asc, betax.y.asc, betax.z.asc, betay.d.asc, betay.maximum.asc, betay.minimum.asc, betay.norm1.asc, betay.norm2.asc, betay.x.asc, betay.y.asc, betay.z.asc, betaz.d.asc, betaz.maximum.asc, betaz.minimum.asc, betaz.norm1.asc, betaz.norm2.asc, betaz.x.asc, betaz.y.asc, betaz.z.asc, dens.d.asc, dens.maximum.asc, dens.minimum.asc, dens.norm1.asc, dens.norm2.asc, dens.x.asc, dens.y.asc, dens.z.asc, eps.d.asc, eps.maximum.asc, eps.minimum.asc, eps.norm1.asc, eps.norm2.asc, eps.x.asc, eps.y.asc, eps.z.asc, gxx.d.asc, gxx.maximum.asc, gxx.minimum.asc, gxx.norm1.asc, gxx.norm2.asc, gxx.x.asc, gxx.y.asc, gxx.z.asc, gxy.d.asc, gxy.maximum.asc, gxy.minimum.asc, gxy.norm1.asc, gxy.norm2.asc, gxy.x.asc, gxy.y.asc, gxy.z.asc, gxz.d.asc, gxz.maximum.asc, gxz.minimum.asc, gxz.norm1.asc, gxz.norm2.asc, gxz.x.asc, gxz.y.asc, gxz.z.asc, gyy.d.asc, gyy.maximum.asc, gyy.minimum.asc, gyy.norm1.asc, gyy.norm2.asc, gyy.x.asc, gyy.y.asc, gyy.z.asc, gyz.d.asc, gyz.maximum.asc, gyz.minimum.asc, gyz.norm1.asc, gyz.norm2.asc, gyz.x.asc, gyz.y.asc, gyz.z.asc, gzz.d.asc, gzz.maximum.asc, gzz.minimum.asc, gzz.norm1.asc, gzz.norm2.asc, gzz.x.asc, gzz.y.asc, gzz.z.asc, kxx.d.asc, kxx.maximum.asc, kxx.minimum.asc, kxx.norm1.asc, kxx.norm2.asc, kxx.x.asc, kxx.y.asc, kxx.z.asc, kxy.d.asc, kxy.maximum.asc, kxy.minimum.asc, kxy.norm1.asc, kxy.norm2.asc, kxy.x.asc, kxy.y.asc, kxy.z.asc, kxz.d.asc, kxz.maximum.asc, kxz.minimum.asc, kxz.norm1.asc, kxz.norm2.asc, kxz.x.asc, kxz.y.asc, kxz.z.asc, kyy.d.asc, kyy.maximum.asc, kyy.minimum.asc, kyy.norm1.asc, kyy.norm2.asc, kyy.x.asc, kyy.y.asc, kyy.z.asc, kyz.d.asc, kyz.maximum.asc, kyz.minimum.asc, kyz.norm1.asc, kyz.norm2.asc, kyz.x.asc, kyz.y.asc, kyz.z.asc, kzz.d.asc, kzz.maximum.asc, kzz.minimum.asc, kzz.norm1.asc, kzz.norm2.asc, kzz.x.asc, kzz.y.asc, kzz.z.asc, press.d.asc, press.maximum.asc, press.minimum.asc, press.norm1.asc, press.norm2.asc, press.x.asc, press.y.asc, press.z.asc, puncture_u.d.asc, puncture_u.maximum.asc, puncture_u.minimum.asc, puncture_u.norm1.asc, puncture_u.norm2.asc, puncture_u.x.asc, puncture_u.y.asc, puncture_u.z.asc, rho.d.asc, rho.maximum.asc, rho.minimum.asc, rho.norm1.asc, rho.norm2.asc, rho.x.asc, rho.y.asc, rho.z.asc, scon[0].d.asc, scon[0].maximum.asc, scon[0].minimum.asc, scon[0].norm1.asc, scon[0].norm2.asc, scon[0].x.asc, scon[0].y.asc, scon[0].z.asc, scon[1].d.asc, scon[1].maximum.asc, scon[1].minimum.asc, scon[1].norm1.asc, scon[1].norm2.asc, scon[1].x.asc, scon[1].y.asc, scon[1].z.asc, scon[2].d.asc, scon[2].maximum.asc, scon[2].minimum.asc, scon[2].norm1.asc, scon[2].norm2.asc, scon[2].x.asc, scon[2].y.asc, scon[2].z.asc, tau.d.asc, tau.maximum.asc, tau.minimum.asc, tau.norm1.asc, tau.norm2.asc, tau.x.asc, tau.y.asc, tau.z.asc, vel[0].d.asc, vel[0].maximum.asc, vel[0].minimum.asc, vel[0].norm1.asc, vel[0].norm2.asc, vel[0].x.asc, vel[0].y.asc, vel[0].z.asc, vel[1].d.asc, vel[1].maximum.asc, vel[1].minimum.asc, vel[1].norm1.asc, vel[1].norm2.asc, vel[1].x.asc, vel[1].y.asc, vel[1].z.asc, vel[2].d.asc, vel[2].maximum.asc, vel[2].minimum.asc, vel[2].norm1.asc, vel[2].norm2.asc, vel[2].x.asc, vel[2].y.asc, vel[2].z.asc, w_lorentz.d.asc, w_lorentz.maximum.asc, w_lorentz.minimum.asc, w_lorentz.norm1.asc, w_lorentz.norm2.asc, w_lorentz.x.asc, w_lorentz.y.asc, w_lorentz.z.asc

Modified: /branches/momentum_constraint/ param.ccl /branches/momentum_constraint/src/ Source.c, TwoPunctures.c, TwoPunctures.h

Log: TwoPunctures: add option to compute Aij using finite differencing

Adds an option use_finite_differencing_for_Aij (defaults to no) which computes the extrinsic curvature via finite differencing when putting it onto the Cactus grid. Gives smoother results than using the spectral values directly. Test suite (and example parameter file) are included in this commit.

Currently terribly inefficient since it evaluates the spectral representation nine times for each Cactus grid point.

File Changes:

Directory: /branches/momentum_constraint/test/bhns_momentum_fd/ ===============================================================

File [added]: alp.d.asc Delta lines: <Unavailable>

File [added]: alp.maximum.asc Delta lines: <Unavailable>

File [added]: alp.minimum.asc Delta lines: <Unavailable>

File [added]: alp.norm1.asc Delta lines: <Unavailable>

File [added]: alp.norm2.asc Delta lines: <Unavailable>

File [added]: alp.x.asc Delta lines: <Unavailable>

File [added]: alp.y.asc Delta lines: <Unavailable>

File [added]: alp.z.asc Delta lines: <Unavailable>

File [added]: betax.d.asc Delta lines: <Unavailable>

File [added]: betax.maximum.asc Delta lines: <Unavailable>

File [added]: betax.minimum.asc Delta lines: <Unavailable>

File [added]: betax.norm1.asc Delta lines: <Unavailable>

File [added]: betax.norm2.asc Delta lines: <Unavailable>

File [added]: betax.x.asc Delta lines: <Unavailable>

File [added]: betax.y.asc Delta lines: <Unavailable>

File [added]: betax.z.asc Delta lines: <Unavailable>

File [added]: betay.d.asc Delta lines: <Unavailable>

File [added]: betay.maximum.asc Delta lines: <Unavailable>

File [added]: betay.minimum.asc Delta lines: <Unavailable>

File [added]: betay.norm1.asc Delta lines: <Unavailable>

File [added]: betay.norm2.asc Delta lines: <Unavailable>

File [added]: betay.x.asc Delta lines: <Unavailable>

File [added]: betay.y.asc Delta lines: <Unavailable>

File [added]: betay.z.asc Delta lines: <Unavailable>

File [added]: betaz.d.asc Delta lines: <Unavailable>

File [added]: betaz.maximum.asc Delta lines: <Unavailable>

File [added]: betaz.minimum.asc Delta lines: <Unavailable>

File [added]: betaz.norm1.asc Delta lines: <Unavailable>

File [added]: betaz.norm2.asc Delta lines: <Unavailable>

File [added]: betaz.x.asc Delta lines: <Unavailable>

File [added]: betaz.y.asc Delta lines: <Unavailable>

File [added]: betaz.z.asc Delta lines: <Unavailable>

File [added]: dens.d.asc Delta lines: <Unavailable>

File [added]: dens.maximum.asc Delta lines: <Unavailable>

File [added]: dens.minimum.asc Delta lines: <Unavailable>

File [added]: dens.norm1.asc Delta lines: <Unavailable>

File [added]: dens.norm2.asc Delta lines: <Unavailable>

File [added]: dens.x.asc Delta lines: <Unavailable>

File [added]: dens.y.asc Delta lines: <Unavailable>

File [added]: dens.z.asc Delta lines: <Unavailable>

File [added]: eps.d.asc Delta lines: <Unavailable>

File [added]: eps.maximum.asc Delta lines: <Unavailable>

File [added]: eps.minimum.asc Delta lines: <Unavailable>

File [added]: eps.norm1.asc Delta lines: <Unavailable>

File [added]: eps.norm2.asc Delta lines: <Unavailable>

File [added]: eps.x.asc Delta lines: <Unavailable>

File [added]: eps.y.asc Delta lines: <Unavailable>

File [added]: eps.z.asc Delta lines: <Unavailable>

File [added]: gxx.d.asc Delta lines: <Unavailable>

File [added]: gxx.maximum.asc Delta lines: <Unavailable>

File [added]: gxx.minimum.asc Delta lines: <Unavailable>

File [added]: gxx.norm1.asc Delta lines: <Unavailable>

File [added]: gxx.norm2.asc Delta lines: <Unavailable>

File [added]: gxx.x.asc Delta lines: <Unavailable>

File [added]: gxx.y.asc Delta lines: <Unavailable>

File [added]: gxx.z.asc Delta lines: <Unavailable>

File [added]: gxy.d.asc Delta lines: <Unavailable>

File [added]: gxy.maximum.asc Delta lines: <Unavailable>

File [added]: gxy.minimum.asc Delta lines: <Unavailable>

File [added]: gxy.norm1.asc Delta lines: <Unavailable>

File [added]: gxy.norm2.asc Delta lines: <Unavailable>

File [added]: gxy.x.asc Delta lines: <Unavailable>

File [added]: gxy.y.asc Delta lines: <Unavailable>

File [added]: gxy.z.asc Delta lines: <Unavailable>

File [added]: gxz.d.asc Delta lines: <Unavailable>

File [added]: gxz.maximum.asc Delta lines: <Unavailable>

File [added]: gxz.minimum.asc Delta lines: <Unavailable>

File [added]: gxz.norm1.asc Delta lines: <Unavailable>

File [added]: gxz.norm2.asc Delta lines: <Unavailable>

File [added]: gxz.x.asc Delta lines: <Unavailable>

File [added]: gxz.y.asc Delta lines: <Unavailable>

File [added]: gxz.z.asc Delta lines: <Unavailable>

File [added]: gyy.d.asc Delta lines: <Unavailable>

File [added]: gyy.maximum.asc Delta lines: <Unavailable>

File [added]: gyy.minimum.asc Delta lines: <Unavailable>

File [added]: gyy.norm1.asc Delta lines: <Unavailable>

File [added]: gyy.norm2.asc Delta lines: <Unavailable>

File [added]: gyy.x.asc Delta lines: <Unavailable>

File [added]: gyy.y.asc Delta lines: <Unavailable>

File [added]: gyy.z.asc Delta lines: <Unavailable>

File [added]: gyz.d.asc Delta lines: <Unavailable>

File [added]: gyz.maximum.asc Delta lines: <Unavailable>

File [added]: gyz.minimum.asc Delta lines: <Unavailable>

File [added]: gyz.norm1.asc Delta lines: <Unavailable>

File [added]: gyz.norm2.asc Delta lines: <Unavailable>

File [added]: gyz.x.asc Delta lines: <Unavailable>

File [added]: gyz.y.asc Delta lines: <Unavailable>

File [added]: gyz.z.asc Delta lines: <Unavailable>

File [added]: gzz.d.asc Delta lines: <Unavailable>

File [added]: gzz.maximum.asc Delta lines: <Unavailable>

File [added]: gzz.minimum.asc Delta lines: <Unavailable>

File [added]: gzz.norm1.asc Delta lines: <Unavailable>

File [added]: gzz.norm2.asc Delta lines: <Unavailable>

File [added]: gzz.x.asc Delta lines: <Unavailable>

File [added]: gzz.y.asc Delta lines: <Unavailable>

File [added]: gzz.z.asc Delta lines: <Unavailable>

File [added]: kxx.d.asc Delta lines: <Unavailable>

File [added]: kxx.maximum.asc Delta lines: <Unavailable>

File [added]: kxx.minimum.asc Delta lines: <Unavailable>

File [added]: kxx.norm1.asc Delta lines: <Unavailable>

File [added]: kxx.norm2.asc Delta lines: <Unavailable>

File [added]: kxx.x.asc Delta lines: <Unavailable>

File [added]: kxx.y.asc Delta lines: <Unavailable>

File [added]: kxx.z.asc Delta lines: <Unavailable>

File [added]: kxy.d.asc Delta lines: <Unavailable>

File [added]: kxy.maximum.asc Delta lines: <Unavailable>

File [added]: kxy.minimum.asc Delta lines: <Unavailable>

File [added]: kxy.norm1.asc Delta lines: <Unavailable>

File [added]: kxy.norm2.asc Delta lines: <Unavailable>

File [added]: kxy.x.asc Delta lines: <Unavailable>

File [added]: kxy.y.asc Delta lines: <Unavailable>

File [added]: kxy.z.asc Delta lines: <Unavailable>

File [added]: kxz.d.asc Delta lines: <Unavailable>

File [added]: kxz.maximum.asc Delta lines: <Unavailable>

File [added]: kxz.minimum.asc Delta lines: <Unavailable>

File [added]: kxz.norm1.asc Delta lines: <Unavailable>

File [added]: kxz.norm2.asc Delta lines: <Unavailable>

File [added]: kxz.x.asc Delta lines: <Unavailable>

File [added]: kxz.y.asc Delta lines: <Unavailable>

File [added]: kxz.z.asc Delta lines: <Unavailable>

File [added]: kyy.d.asc Delta lines: <Unavailable>

File [added]: kyy.maximum.asc Delta lines: <Unavailable>

File [added]: kyy.minimum.asc Delta lines: <Unavailable>

File [added]: kyy.norm1.asc Delta lines: <Unavailable>

File [added]: kyy.norm2.asc Delta lines: <Unavailable>

File [added]: kyy.x.asc Delta lines: <Unavailable>

File [added]: kyy.y.asc Delta lines: <Unavailable>

File [added]: kyy.z.asc Delta lines: <Unavailable>

File [added]: kyz.d.asc Delta lines: <Unavailable>

File [added]: kyz.maximum.asc Delta lines: <Unavailable>

File [added]: kyz.minimum.asc Delta lines: <Unavailable>

File [added]: kyz.norm1.asc Delta lines: <Unavailable>

File [added]: kyz.norm2.asc Delta lines: <Unavailable>

File [added]: kyz.x.asc Delta lines: <Unavailable>

File [added]: kyz.y.asc Delta lines: <Unavailable>

File [added]: kyz.z.asc Delta lines: <Unavailable>

File [added]: kzz.d.asc Delta lines: <Unavailable>

File [added]: kzz.maximum.asc Delta lines: <Unavailable>

File [added]: kzz.minimum.asc Delta lines: <Unavailable>

File [added]: kzz.norm1.asc Delta lines: <Unavailable>

File [added]: kzz.norm2.asc Delta lines: <Unavailable>

File [added]: kzz.x.asc Delta lines: <Unavailable>

File [added]: kzz.y.asc Delta lines: <Unavailable>

File [added]: kzz.z.asc Delta lines: <Unavailable>

File [added]: press.d.asc Delta lines: <Unavailable>

File [added]: press.maximum.asc Delta lines: <Unavailable>

File [added]: press.minimum.asc Delta lines: <Unavailable>

File [added]: press.norm1.asc Delta lines: <Unavailable>

File [added]: press.norm2.asc Delta lines: <Unavailable>

File [added]: press.x.asc Delta lines: <Unavailable>

File [added]: press.y.asc Delta lines: <Unavailable>

File [added]: press.z.asc Delta lines: <Unavailable>

File [added]: puncture_u.d.asc Delta lines: <Unavailable>

File [added]: puncture_u.maximum.asc Delta lines: <Unavailable>

File [added]: puncture_u.minimum.asc Delta lines: <Unavailable>

File [added]: puncture_u.norm1.asc Delta lines: <Unavailable>

File [added]: puncture_u.norm2.asc Delta lines: <Unavailable>

File [added]: puncture_u.x.asc Delta lines: <Unavailable>

File [added]: puncture_u.y.asc Delta lines: <Unavailable>

File [added]: puncture_u.z.asc Delta lines: <Unavailable>

File [added]: rho.d.asc Delta lines: <Unavailable>

File [added]: rho.maximum.asc Delta lines: <Unavailable>

File [added]: rho.minimum.asc Delta lines: <Unavailable>

File [added]: rho.norm1.asc Delta lines: <Unavailable>

File [added]: rho.norm2.asc Delta lines: <Unavailable>

File [added]: rho.x.asc Delta lines: <Unavailable>

File [added]: rho.y.asc Delta lines: <Unavailable>

File [added]: rho.z.asc Delta lines: <Unavailable>

File [added]: scon[0].d.asc Delta lines: <Unavailable>

File [added]: scon[0].maximum.asc Delta lines: <Unavailable>

File [added]: scon[0].minimum.asc Delta lines: <Unavailable>

File [added]: scon[0].norm1.asc Delta lines: <Unavailable>

File [added]: scon[0].norm2.asc Delta lines: <Unavailable>

File [added]: scon[0].x.asc Delta lines: <Unavailable>

File [added]: scon[0].y.asc Delta lines: <Unavailable>

File [added]: scon[0].z.asc Delta lines: <Unavailable>

File [added]: scon[1].d.asc Delta lines: <Unavailable>

File [added]: scon[1].maximum.asc Delta lines: <Unavailable>

File [added]: scon[1].minimum.asc Delta lines: <Unavailable>

File [added]: scon[1].norm1.asc Delta lines: <Unavailable>

File [added]: scon[1].norm2.asc Delta lines: <Unavailable>

File [added]: scon[1].x.asc Delta lines: <Unavailable>

File [added]: scon[1].y.asc Delta lines: <Unavailable>

File [added]: scon[1].z.asc Delta lines: <Unavailable>

File [added]: scon[2].d.asc Delta lines: <Unavailable>

File [added]: scon[2].maximum.asc Delta lines: <Unavailable>

File [added]: scon[2].minimum.asc Delta lines: <Unavailable>

File [added]: scon[2].norm1.asc Delta lines: <Unavailable>

File [added]: scon[2].norm2.asc Delta lines: <Unavailable>

File [added]: scon[2].x.asc Delta lines: <Unavailable>

File [added]: scon[2].y.asc Delta lines: <Unavailable>

File [added]: scon[2].z.asc Delta lines: <Unavailable>

File [added]: tau.d.asc Delta lines: <Unavailable>

File [added]: tau.maximum.asc Delta lines: <Unavailable>

File [added]: tau.minimum.asc Delta lines: <Unavailable>

File [added]: tau.norm1.asc Delta lines: <Unavailable>

File [added]: tau.norm2.asc Delta lines: <Unavailable>

File [added]: tau.x.asc Delta lines: <Unavailable>

File [added]: tau.y.asc Delta lines: <Unavailable>

File [added]: tau.z.asc Delta lines: <Unavailable>

File [added]: vel[0].d.asc Delta lines: <Unavailable>

File [added]: vel[0].maximum.asc Delta lines: <Unavailable>

File [added]: vel[0].minimum.asc Delta lines: <Unavailable>

File [added]: vel[0].norm1.asc Delta lines: <Unavailable>

File [added]: vel[0].norm2.asc Delta lines: <Unavailable>

File [added]: vel[0].x.asc Delta lines: <Unavailable>

File [added]: vel[0].y.asc Delta lines: <Unavailable>

File [added]: vel[0].z.asc Delta lines: <Unavailable>

File [added]: vel[1].d.asc Delta lines: <Unavailable>

File [added]: vel[1].maximum.asc Delta lines: <Unavailable>

File [added]: vel[1].minimum.asc Delta lines: <Unavailable>

File [added]: vel[1].norm1.asc Delta lines: <Unavailable>

File [added]: vel[1].norm2.asc Delta lines: <Unavailable>

File [added]: vel[1].x.asc Delta lines: <Unavailable>

File [added]: vel[1].y.asc Delta lines: <Unavailable>

File [added]: vel[1].z.asc Delta lines: <Unavailable>

File [added]: vel[2].d.asc Delta lines: <Unavailable>

File [added]: vel[2].maximum.asc Delta lines: <Unavailable>

File [added]: vel[2].minimum.asc Delta lines: <Unavailable>

File [added]: vel[2].norm1.asc Delta lines: <Unavailable>

File [added]: vel[2].norm2.asc Delta lines: <Unavailable>

File [added]: vel[2].x.asc Delta lines: <Unavailable>

File [added]: vel[2].y.asc Delta lines: <Unavailable>

File [added]: vel[2].z.asc Delta lines: <Unavailable>

File [added]: w_lorentz.d.asc Delta lines: <Unavailable>

File [added]: w_lorentz.maximum.asc Delta lines: <Unavailable>

File [added]: w_lorentz.minimum.asc Delta lines: <Unavailable>

File [added]: w_lorentz.norm1.asc Delta lines: <Unavailable>

File [added]: w_lorentz.norm2.asc Delta lines: <Unavailable>

File [added]: w_lorentz.x.asc Delta lines: <Unavailable>

File [added]: w_lorentz.y.asc Delta lines: <Unavailable>

File [added]: w_lorentz.z.asc Delta lines: <Unavailable>

Directory: /branches/momentum_constraint/ =========================================

File [modified]: param.ccl Delta lines: <Unavailable>

Directory: /branches/momentum_constraint/src/ =============================================

File [modified]: Source.c Delta lines: <Unavailable>

File [modified]: TwoPunctures.c Delta lines: <Unavailable>

File [modified]: TwoPunctures.h Delta lines: <Unavailable>

Directory: /branches/momentum_constraint/test/ ==============================================

File [added]: bhns_momentum_fd.par Delta lines: <Unavailable>