From b0ae812addbd9b8ae4c6815484e5ff366bd05016 Mon Sep 17 00:00:00 2001 From: Thibaut Paumard Date: Tue, 22 Oct 2019 12:52:57 +0200 Subject: [PATCH] Preparing to release 1.4.0 --- ChangeLog | 760 ++++++++++++++++++++++++++++++++++++++++++++++++++- NEWS | 36 ++- configure | 24 +- configure.ac | 6 +- 4 files changed, 806 insertions(+), 20 deletions(-) diff --git a/ChangeLog b/ChangeLog index 995bffe8..b9785712 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,6 +1,36 @@ commit unknown Author: Thibaut Paumard -Date: Fri Sep 27 16:03:01 2019 +0200 +Date: Tue Oct 22 12:52:37 2019 +0200 + + Releasing 1.4.0 + +commit eb1da8b77098bbb8ac8035312c3634079aea9e9c +Author: Thibaut Paumard +Date: Tue Oct 22 12:50:35 2019 +0200 + + Fix documentation in Register and Scenery + +commit 1e30036d8d4b216ea25cb2a35fe1367d8ea3f7c9 +Author: Thibaut Paumard +Date: Sat Aug 3 09:13:51 2019 +0200 + + Releasing 1.3.5 + +commit 31195c8f583085d36911ae0b567019a90d83c108 +Author: Thibaut Paumard +Date: Fri Aug 2 22:37:38 2019 +0200 + + Releasing 1.3.4 (fix merge in master) + +commit d75a030dec028dd4f6f44e478db56d666d9d502f +Author: Thibaut Paumard +Date: Fri Aug 2 20:53:30 2019 +0200 + + Releasing 1.3.3 (fix import in master) + +commit 904c849c86fd9552837d1bd55efad9fb5b88cd06 +Author: Thibaut Paumard +Date: Fri Sep 27 16:04:19 2019 +0200 Releasing 1.3.6 @@ -31,6 +61,40 @@ Date: Fri Sep 27 14:37:09 2019 +0200 Fix example.py after changing Worldline.getCoord +commit ae6b2131a1bbe7ebdbfc757070985c395aa29787 +Author: Frederic Vincent +Date: Fri Sep 27 08:49:35 2019 +0200 + + ThickDisk: commiting proper versions of StdPlug and lib/Makefile.in + +commit c2ec8e22fed7319d3599fb0401cd598fec03a96e +Author: Frederic Vincent +Date: Fri Sep 27 08:20:25 2019 +0200 + + New Astrobj ThickDisk + This is a geometrically thick, optically thin disk, described + by its opening angle (between BH spin axis and outer disk surface), + and its inner radius. It emits thermal synchrotron parametrized + by an inner temperature, inner density, and a magnetization. + Its velocity is Keplerian outside ISCO, and can be tuned + inside ISCO (see VelocityBelowIsco) + +commit d74acb8469a716a77a9b56b2b9fa4af9aaffbc98 +Author: Frederic Vincent +Date: Fri Sep 27 08:15:18 2019 +0200 + + Jet: add a Thermal Synchrotron spectrum and change jet velocity to (Vphi,Vr) + Velocity stuff: previous implementation was assuming pure radial + velocity in jet. Now a suerimposition of azimuthal and radial + is allowed with Vphi/Vr chosen in xml + +commit f230d7fb226911fe5679fc166fef3f82536d08c8 +Author: Frederic Vincent +Date: Fri Sep 27 08:11:31 2019 +0200 + + in Torus: add a Thermal Synchrotron spectrum and a radiativeQ() + recipe to allow realistic astrophysics + commit f277f6ebf3cc3810c7b56bd0b779540723272545 Author: Thibaut Paumard Date: Sun Aug 25 19:32:33 2019 +0200 @@ -125,6 +189,137 @@ Date: Fri Aug 9 09:37:16 2019 +0200 python: install util.rayTrace as Scenery.rayTrace +commit a71e2ef9f529a9c0e49a662d90b8fedf969d5a07 +Author: Thibaut Paumard +Date: Tue Sep 10 10:25:36 2019 +0200 + + GyotoDefs.h: give derived constants to machine precision + + Some constants like GYOTO_G_OVER_C_SQUARED where only given to the same + precision as the constant they derive from (GYOTO_G in this case). + + While it is correct that e.g. G is known only to 1e5, mixiin guses of GYOTO_G + and GYOTO_G_OVER_C_SQUARED both specified only to that precision leads to + additional errors. GYOTO_G_OVER_C_SQUARE, GYOTO_G_OVER_C_SQUARE_CGS, + GYOTO_PLANCK_OVER_C_SQUARE and GYOTO_PLANCK_OVER_BOLTZMANN are now given with + many more places so that GYOTO_G/GYOTO_G_OVER_C_SQUARE actually yields GYOTO_C^2. + +commit 933632b26e4138b1f9c6ac11dde3ea53eb89aaf9 +Author: Thibaut Paumard +Date: Thu Sep 5 10:41:17 2019 +0200 + + doc/user_guide: correct example about setting the observer tetrad + +commit aa0c11e9953f89c687be267548f5152757e13e76 +Author: Thibaut Paumard +Date: Sun Aug 25 19:36:06 2019 +0200 + + Worldline: fix bug introduced with proper + + time_ was not properly upated after xExpand (or xFill, which calls xExpand) + +commit 770de5b9e28c7490e6e3297fec4337e7609fe4c6 +Author: Thibaut Paumard +Date: Sun Aug 25 19:33:48 2019 +0200 + + GyotoDefs.h: introduce GYOTO_DEBUG_THIS and GYOTO_DEBUG_THIS_EXPR + + Those versions output unconditionnally and are meant for temporary use during debugging sessions. + +commit 59b68193fcb2947bdb80c47ddc9debd3403448dc +Author: Thibaut Paumard +Date: Sun Aug 25 19:32:33 2019 +0200 + + gyoto.util: make getCoord(double t, stat_t coord) work again + +commit 9f7d8453019f49efba9fb466b0c359eef2d62a5a +Author: Thibaut Paumard +Date: Fri Aug 16 09:32:15 2019 +0200 + + Python: make sure Worldline.getCoord works both with NumPy and C arrays + +commit 930d25d714354771e4b5ff2380905991fd22f52a +Author: Thibaut Paumard +Date: Thu Aug 15 22:35:51 2019 +0200 + + Python: make sure Worldline.getCartesian works both with NumPy and C arrays + +commit c77a243dbe042703afe59eb51b4a5ffda65b9baa +Author: Thibaut Paumard +Date: Thu Aug 15 09:01:36 2019 +0200 + + gyoto.util.rayTrace: return arrays or the right dimension + + if j and/or i are scalars, decrease dimension of output arrays + +commit 13ab900e9371bfe34dc389d69e53bda3da12d089 +Author: Thibaut Paumard +Date: Wed Aug 14 21:22:25 2019 +0200 + + gyoto.util.rayTrace: reorder parameters + + The order j, i is more Pythonic than i.j. This ways, + sc[j, i]['Quantity'] == sc[:,:]['Quantity'][j,i]. + + Having height and width last allows writing: + sc[j, i, core.Bucket, None] + skipping height and width to remove progress output. + +commit ec2095da7da28c171f639d4a25e93b8b97b8eb25 +Author: Thibaut Paumard +Date: Wed Aug 14 14:31:42 2019 +0200 + + Gyoto::Screen::Coord1dSet: copy beffur in Indices and Angles + + Much safer in Python + +commit 7f685adf5b748ef908709f170957846ecf0984b3 +Author: Thibaut Paumard +Date: Sun Aug 11 14:07:43 2019 +0200 + + gyoto.util: improve doc comments + +commit b010c42ef536885d3096d54206ba797b1b2c3d93 +Author: Thibaut Paumard +Date: Sun Aug 11 14:00:19 2019 +0200 + + gyoto.util.rayTrace/Scenery_getitem: improve dimensionality + + - Scenery_getitem accepts all rayTrace parameters + - return 1D array when the Coord2dSet is 1D + +commit 6da372f854ce2146a9093f122b58300df1644abf +Author: Thibaut Paumard +Date: Sun Aug 11 07:53:12 2019 +0200 + + gyoto.util: support angles in rayTrace + +commit 0fd2ad270afc7616408211439e5bc8d6b23ad216 +Author: Thibaut Paumard +Date: Sun Aug 11 07:45:42 2019 +0200 + + gyoto.core.Scenery: add __getitem__ method + + shortcut for rayTrace + +commit 7094bfe1607f0b9ad8be4d49c9d3487718204c80 +Author: Thibaut Paumard +Date: Sun Aug 11 07:21:59 2019 +0200 + + gyoto.util.rayTrace: support indices, ranges and list of indices + +commit 194b3c919b1d7a050b2c69bbd6a11541656458ca +Author: Thibaut Paumard +Date: Sun Aug 11 06:21:56 2019 +0200 + + gyoto.util.rayTrace: accept paremeters i and j + +commit 15d6d9ccaa0050267b8b823d62f9ba685eb16d36 +Author: Thibaut Paumard +Date: Fri Aug 9 09:37:16 2019 +0200 + + python: install util.rayTrace as Scenery.rayTrace + commit bb7ce60adc3c9f476dafb0429b35df57ae99fe0a Author: Thibaut Paumard Date: Sat Aug 3 09:13:51 2019 +0200 @@ -253,6 +448,393 @@ Date: Thu Jan 25 11:31:02 2018 +0100 Fixup: specify encoding in python/tests/core.py +commit 1998eb3b8ef1e69f5f047fbd7b100ec2639e5e73 +Author: Thibaut Paumard +Date: Fri Aug 2 14:54:31 2019 +0200 + + Fix printing issues during ray-tracing, mostly visible in Python + + - Make sure cout is flushed after Scenery::rayTrace finishes + - Screen::Grid now keeps a copy of the prefix string + - util.rayTrace fmt defaults to "\r j = " + +commit 6ab168f94799ef8c668c665d452dfa3564a128ab +Author: Thibaut Paumard +Date: Fri Aug 2 14:24:12 2019 +0200 + + Screen.C: apply PALN also for other observerKinds + +commit 93e49c0fef6ccdd5dc6753b00d34d33ec095c3bb +Author: Thibaut Paumard +Date: Thu Aug 1 13:56:29 2019 +0200 + + gyoto.animate: add static_screen video kind + +commit 964d581ded5c73c46dc11d6d5cbe3de429dae25a +Author: Thibaut Paumard +Date: Wed Jul 31 15:32:34 2019 +0200 + + doc/user_guide: document how to set observer tetrad + +commit 073d516079a50966cdf712e8e075111cc71a6768 +Author: Thibaut Paumard +Date: Tue Jul 30 16:54:55 2019 +0200 + + doc: update documentation on the Python modules + +commit e1a2779bd7e4bdc7b611cac0f5a04fe906495ba3 +Author: Thibaut Paumard +Date: Tue Jul 30 16:22:38 2019 +0200 + + doc: warn in the user guide that the Yorick plug-in should not be used + + The Yorick plug-in will not see new features. It's maintained for gyotoy only. + +commit ee82c057d2208306b33a403156a6179ff388c7cf +Author: Thibaut Paumard +Date: Tue Jul 30 15:50:11 2019 +0200 + + gyoto.animate: add NullVideoWriter, roll angle in _forward + +commit a46ff58d3da3f6569fc3437974e1d03681a378c9 +Author: Thibaut Paumard +Date: Tue Jul 30 11:54:31 2019 +0200 + + mk_video.C: actually include it in git + +commit 68024b0dabfd75c61e541ca9717ff9fdaf087fc6 +Author: Thibaut Paumard +Date: Tue Jul 30 11:42:48 2019 +0200 + + gyoto.animate: improve projection in orbiting_screen_forward + +commit 84a70b3fb343f69a72f83a189ee337e7bb9d8c9e +Author: Thibaut Paumard +Date: Mon Jul 29 17:17:25 2019 +0200 + + gyoto.animate: add orbiting_screen_forward + +commit 4b046b5bc964ff3823818555e3a62c54cccc9459 +Author: Thibaut Paumard +Date: Mon Jul 29 17:15:38 2019 +0200 + + Metric: fix bug in Generic::cartesianVelocity + + Was only correct in the equatorial plane + +commit 3ffdef6805430eef2b6f60185b51d6bf0773464a +Author: Thibaut Paumard +Date: Mon Jul 29 14:44:35 2019 +0200 + + Metric: fix projectFourVect again + + Previous fix was worse than original bug. + + Also introduce test case in python/tests + +commit 1e80b1de6c001d637105effd050fe51664d772df +Author: Thibaut Paumard +Date: Sun Jul 28 22:21:52 2019 +0200 + + Metric.C: remove debugging cerr + +commit 2337fd3ac575dddaa1271f63f5c9f5de78e581d3 +Author: Thibaut Paumard +Date: Sun Jul 28 19:59:18 2019 +0200 + + Metric: fix zamoVelocity and projectFourVect + + zamoVelocity should now be correct in all cases. We now take + (1, 0, 0, 0) and project it onto the hyperplan orthogonal to + ephi. + + projectFourVect had a bug and pnly worked when u2 had norm a + norm of + or -1. Fixed. + +commit 512e484ae84dc78597791ccdcdeeee5c3f862271 +Author: Thibaut Paumard +Date: Fri Jul 26 22:52:28 2019 +0200 + + Metric: make observerTetrad work for any metric + + + Metric: new helper methods norm, multiplyFourVect, addFourVect, + projectFourVect, GramSchmidt; renormalize fourvel in observerTetrad. + + + KerrBL: renormalize 4-velocity in observerTetrad + +commit 6dac2074a07b8d2fb47b2d674b08ef0f0ae087f4 +Author: Thibaut Paumard +Date: Fri Jul 26 22:51:13 2019 +0200 + + ThinDisk::sphericalPhi: return value between 0 and 2 pi + +commit 5fe2b21be1957a426a29a7fee8d09765c1a72e4b +Author: Thibaut Paumard +Date: Fri Jul 26 22:50:19 2019 +0200 + + Python: debug some methods in Worldline + +commit 250d256006d30bea6b6db5e74a8760093e18c4a8 +Author: Thibaut Paumard +Date: Fri Jul 26 22:49:35 2019 +0200 + + gyoto.util: add readScenery and writeObject + +commit 7316ec87d22a57d4dda85db79c86bda2b34cd0b3 +Author: Thibaut Paumard +Date: Fri Jul 26 15:37:05 2019 +0200 + + python: gyoto.animate.rayTraceFrame: width/height override resolution + +commit 214ca5466b1830cbe3f11c9f31f719fa28b91fb1 +Author: Thibaut Paumard +Date: Fri Jul 26 15:35:19 2019 +0200 + + python: new module gyoto.util containing rayTrace() + + rayTrace is a utility function to hide the complexity of calling Scenery.rayTrace + +commit 0f935bd7712502886f5999620703c8eddad0b294 +Author: Thibaut Paumard +Date: Thu Jul 25 14:34:11 2019 +0200 + + gyoto.C: allow making movies directly from the gyoto command line + + gyoto mk-video [options] will try to load the Python plug-in + and use it to run code in the gyoto.animate module. + +commit 7f3649ebea23be6d38a40eca656732a39714d667 +Author: Thibaut Paumard +Date: Wed Jul 24 15:03:08 2019 +0200 + + animate.py: add PyAVVideoWriter, add cmap option + +commit c34397fce9d01c2bb5bacde0373152aa319bcfee +Author: Thibaut Paumard +Date: Wed Jul 24 13:48:40 2019 +0200 + + python: call setup.py build instead of build_ext + +commit 83495797cbb52e50617cb34ca2263db39ff46bf9 +Author: Thibaut Paumard +Date: Wed Jul 24 12:34:28 2019 +0200 + + Screen: observerkind is now stored internally as an unsigned int + +commit e4c077a0eeeca6064e33ef527ab1a3dcc5c2af0a +Author: Thibaut Paumard +Date: Wed Jul 24 12:33:52 2019 +0200 + + ThinDisk: add VelocityKind property to choose between Keplerian and ZAMO + +commit 7202aced515b5e4eaa505a55c7bdcb58ef075e4a +Author: Thibaut Paumard +Date: Wed Jul 24 10:47:53 2019 +0200 + + Metric: add normalizeFourVel + + Also fix a bug with normalization in Screen + +commit 6741faf27ea4e4d02c4768e2c11ffb6e754213ac +Author: Thibaut Paumard +Date: Tue Jul 23 20:35:32 2019 +0200 + + animate.py: pass nframes to func + +commit e7a5e281b402f71a96bf46ee438355834ad508a2 +Author: Thibaut Paumard +Date: Tue Jul 23 16:06:35 2019 +0200 + + Astrobj: allow changing deltaMax inside rMax + +commit facb69c9d1d16b018b5a8f84d686e17c5f5e1899 +Author: Thibaut Paumard +Date: Mon Jul 22 21:49:17 2019 +0200 + + animate.py: add growing mass, change rin to 0 + +commit 6a5d56304ea35d42abf816fec4263ae8bdbf1c0c +Author: Thibaut Paumard +Date: Fri Jul 19 17:16:46 2019 +0200 + + Keep track of proper time in integration + + * Worldline: + + New member tau_ with accessor get_tau; + + xStore take additional argument tau; + + Methods that take time as argument take new argument 'proper': + various versions of getCoord, xFill. + + * Worldline::IntegState: + + nextStep takes additionnal parameter tau and updates it + (not implemented for Legacy integrator). + + * Metric::Generic: + + renormalize fourvel in observerTetrad as it may now come from an + integration to later be able to check orthonormality of the tetrad + + turn orthonormality check from error to warning. + + * Photon: + + integrate and store affine parameter in hit. + +commit d8bd05b2c85b2a1b7b5be6860500ec57f8fcfa03 +Author: Thibaut Paumard +Date: Tue Jul 16 19:37:41 2019 +0200 + + python: new module gyoto.animate for creating Gyoto-based movies + +commit f822f4dd72176daf5507b0ea1cf47c1ef87181be +Author: Thibaut Paumard +Date: Tue Jul 16 19:35:12 2019 +0200 + + bugfix: fix norm of ephi in KerrBL::observerTetrad + +commit d5100e96e353b1ef7f0f2cd426b37b0829fa3e66 +Author: Thibaut Paumard +Date: Sun Jul 14 12:08:22 2019 +0200 + + Support computing Screen tetrad from arbitrary velocity + + * Metric: + + new API zamoVelocity(), default assumes null 3-velocity. + + new method dualOneForm(). + + new API observerTetrad(in pos, in fourvel, out screen1,2,3), + default throws an Error. + * KerrBL: + + implement zamoVelocity(); + + remove observerTetrad(obskind...); + + implement observerTetrad(pos, fourvel, ...) for any fourvel. + * Screen: supports two new obskinds: VelocitySpecified and + FullySpecified. + +commit 032dcc8851fe40d62b0090c978e8bca302420c9f +Author: Thibaut Paumard +Date: Sat Jul 13 00:31:35 2019 +0200 + + Screen.C: change definition of dangle(1|2) for anglekind != spherical + + Older definition only works for spherical angles (at most?) + + This definition actually rotates the camera. + +commit 6f60d0318d522dbd119b803d575e27bb5afae2b3 +Author: Thibaut Paumard +Date: Fri Jul 12 10:06:34 2019 +0200 + + python: fix segfault in gyoto.lorene + + numpy was not properly initialized in gyoto_lorene.i + +commit 71501663d3d5d7baa6208ae437de82df9f3726aa +Author: Thibaut Paumard +Date: Thu Jul 11 21:10:46 2019 +0200 + + python: fix segfault in gyoto.std + + numpy was not properly initialized in gyoto_std.i + +commit 8d3a7b2c308581ff4406f519337c09c0bf317006 +Author: Thibaut Paumard +Date: Thu Jul 11 15:26:54 2019 +0200 + + python: add example-patterndisk.py + +commit 0f16e0cb2c8ab9bd5b96e3053b0868325403c1c7 +Author: Frederic Vincent +Date: Thu May 16 15:44:25 2019 +0200 + + Blob.C: update signature of radiativeQ + +commit ab742924fcb73e0234cceee38324a926d1825f89 +Author: Thibaut Paumard +Date: Mon May 6 11:51:19 2019 +0200 + + Update INSTALL.Gyoto.md with some information specific to Debian derivatives + +commit 652932021efef69d6266873bb03db3c0c8acb47e +Author: Frederic Vincent +Date: Tue Apr 16 15:57:44 2019 +0200 + + Website update for Picture Gallery + +commit 281c1619d50d538854e0e33b3cecf88d66765c06 +Author: Frederic Vincent +Date: Tue Apr 16 15:55:04 2019 +0200 + + Again website update for Picture Gallery + +commit feeb8b475b91a8feb3b2fd1293a7c5200f3d1c6b +Author: Frederic Vincent +Date: Tue Apr 16 15:51:32 2019 +0200 + + Website update for Picture Gallery + +commit a55baa2b5b2e315b466801166bfa6a59411c0965 +Author: Frederic Vincent +Date: Tue Apr 16 15:37:54 2019 +0200 + + Update website with picture gallery + +commit e38618af5d599b12f74f905aa40a373e12016594 +Author: Frederic Vincent +Date: Tue Mar 26 14:40:35 2019 +0100 + + FlaredDiskSynchro: adding z variation of density and temperature + +commit 8ea563747715c34734f2c01cb45aec9b79bb97a1 +Author: Frederic Vincent +Date: Fri Mar 1 16:19:15 2019 +0100 + + FlaredDisk: correct copy constructor for multithreading + +commit c8b09e7f9dc1d7fef4404f409196d1eae7b07516 +Author: Frederic Vincent +Date: Thu Feb 28 14:16:07 2019 +0100 + + GridData2D: bilinear interpolation outside time domain + +commit 965a467aa329d07f495302013bdf4e351b9f3e05 +Author: Frederic Vincent +Date: Thu Feb 28 12:09:37 2019 +0100 + + Adding Makefiles and configure + +commit d75c68b8bd521adb2af5c7aa4d2b3bc73c6a2c1e +Author: Frederic Vincent +Date: Thu Feb 28 11:54:10 2019 +0100 + + FlaredDisk: adding time translation and temperature variation + +commit 8b369a20483639b283060c5ef22206254318592b +Author: Frederic Vincent +Date: Wed Feb 27 16:40:08 2019 +0100 + + GridData2D: adding interpolating function + +commit 68a39f6e653e350f976b6c2a3126563e64304bd4 +Author: Frederic Vincent +Date: Fri Feb 22 09:59:24 2019 +0100 + + Correct Makefiles in lib + +commit bde487c5159375d76fb9518c36a84127bbc623f5 +Author: Frederic Vincent +Date: Thu Feb 21 18:11:06 2019 +0100 + + New class GridData2D and subclass FlaredDiskSynchrotron + + GridData2D: allows to deal with 2D (r,phi) data in a grid + with time dependence. The class can read and write from and + to FITS files with a python interface. + The data are supposed to be on a regular grid in all dimensions. + + FlaredDiskSynchrotron: very preliminary submission of a class + using GridData2D to read density and velocity from a 2D grid + in the equatorial plane. The density will be analytically + prescribed in the z dimension, with a specified aspect ratio z/r. + The flared disk will emit optically thin synchrotron radiation. + Not ready yet for doing physics. + commit 3534961b0a31f2e2975e3786391dc78407cd956f Author: Thibaut Paumard Date: Tue Jan 22 14:03:23 2019 +0100 @@ -288,6 +870,42 @@ Date: Mon Jan 21 17:33:02 2019 +0100 This macro prepends the error message with __FILE__:__LINE__ in __PRETTY_FUNCTION__ +commit f31bec703a0170eca81ec69bd427612d92707df3 +Author: Thibaut Paumard +Date: Mon Jan 21 17:33:02 2019 +0100 + + New macro GYOTO_ERROR around throwError + + This macro prepends the error message with __FILE__:__LINE__ in __PRETTY_FUNCTION__ + +commit 0e9850d9c08f94b8cceae99562c1d57825fc225e +Author: Thibaut Paumard +Date: Mon Jan 21 15:05:01 2019 +0100 + + UniformSphere: deprecate Alpha + + This was a tweak to implement the same thing as integrateEmission more + efficiently in a specific situation. + The more generic implementation still works. + +commit 31caeab224aa42eb5ed1ff5f3e689aa04d2e034f +Author: Thibaut Paumard +Date: Sat Jan 19 13:54:10 2019 +0100 + + Astrobj: let emission/transmission/radiativeQ communicate + + The default implementations of these methods call each-other + and use a private member to cache whether each method is the + default implementation. + + transmission() has now the same signature as emission(). + +commit 987ebfe0b233dff3814d512831b6bb77dcdf2d62 +Author: Thibaut Paumard +Date: Thu Jan 17 14:07:39 2019 +0100 + + Fixup: make sure 'polar' compiles after merging 'master' + commit 7de5ee601f1ac18e95d4f2d11414976e0937d73e Author: Thibaut Paumard Date: Wed Jan 16 11:21:46 2019 +0100 @@ -921,6 +1539,22 @@ Date: Mon Feb 5 20:15:03 2018 +0100 New Metric: Hayward +commit 85b51d6cfa24db506025066fc8529563dbcfd23d +Author: Thibaut Paumard +Date: Wed Jan 31 17:59:25 2018 +0100 + + Test suite: test conservation of base vector orthogonality and norm + + In Python-based test suite. + +commit 85975855b5bad310807deb94266d09dbe177a381 +Author: Thibaut Paumard +Date: Wed Jan 31 17:59:12 2018 +0100 + + ugfix: add time-did-not-evolve check in Worldline::xFill() + + Same test as in Photon::hit(). + commit 8d11cd22c20668155c5a6e9370599438b12a073a Author: Thibaut Paumard Date: Wed Jan 31 11:14:16 2018 +0100 @@ -954,6 +1588,130 @@ Date: Tue Jan 30 00:12:24 2018 +0100 Bugfix: wrong allocated size for Inu in PolishDoughnut::integrateEmission +commit 74bda170afb5721aa9caf1d73773abaf39120308 +Author: Thibaut Paumard +Date: Tue Jan 30 11:15:55 2018 +0100 + + Draft polarized radiative transfer in Astrobj::Generic and Photon + + Provide new interface for polarized radiative transfer (overload of + radiativeQ()). The default implementation calls the unpolarized implementation + to provide consistent unpolarized emission and polarization-neutral + transmission. + + Provide new interface Photon::transfer() to both perform polarized transfer + from the volume element of the emitter to the observer and accumulate the + transfer function of this element into the light-ray transfer function. + This draft implementation actually only treats the unpolarized case. + + Depending on ph -> ParallelTransport, Astrobj::Generic::processHitQuantitive() + calls either this new polarized Astrobj::Generic::radiativeQ() and + Photon::transfer() or use the older unpolarized implementation + (calling unpolarized radiativeQ and Photon::transmit()). + +commit 3a15113c3fb48ae248b3554be94670ceb7da165e +Author: Thibaut Paumard +Date: Tue Jan 30 10:51:14 2018 +0100 + + Bugfix: make PolishDoughnut::emission() and radiativeQ() equivalent + + Previously emission() worked only for komissarov_==false and + radiativeQ only for komissarov_==true. + +commit 4974282751a8aa3900de8c79812115d4ffc60300 +Author: Thibaut Paumard +Date: Tue Jan 30 09:58:33 2018 +0100 + + fixup: correct indices in Worldline::getCoord(double t, state_t &coord) + +commit 05820f1131f7e96a8393585158ad187b893ba929 +Author: Thibaut Paumard +Date: Tue Jan 30 00:12:24 2018 +0100 + + Bugfix: wrong allocated size for Inu in PolishDoughnut::integrateEmission + +commit ea910af29e3752e63a09bccd17ecb4f2b760a2ce +Author: Thibaut Paumard +Date: Mon Jan 29 13:52:37 2018 +0100 + + Prepare Astrobj interface for polarization + + * Change prototype of processHitQuantities; + * Change type of photon coordinate from double[8] to + state_t in all emission/transmission methods; + * Make more parameters const; + * Change typedefs to use a single Gyoto::state_t typedef. + +commit 1966012b9827462be5462edb535fd318ce288928 +Author: Thibaut Paumard +Date: Thu Jan 25 22:30:40 2018 +0100 + + Implement initialization of parallel transport + + Raytracing one of the STOKES quantities should perform parallel transport. + +commit 2f2196f18d0487dd49c2903ad71a8e8408ffe1f8 +Author: Thibaut Paumard +Date: Thu Jan 25 12:27:47 2018 +0100 + + Implement parallel transport of additional vectors + +commit 1d626ad19c22941210b01002d178e5bc70c64586 +Author: Thibaut Paumard +Date: Thu Jan 25 11:33:07 2018 +0100 + + Fixup: change sign of precision in assertAlmostEqual + +commit 634ecab4e2629861a2d17b43da872b9089f965b7 +Author: Thibaut Paumard +Date: Thu Jan 25 11:31:02 2018 +0100 + + Fixup: specify encoding in python/tests/core.py + +commit 69442966bac040b6d905f0099c10da268147249f +Author: Frederic Vincent +Date: Wed Jan 24 18:07:19 2018 +0100 + + Implement the polarization basis (Ephi, Etheta) in Screen + + The polarization basis (Ephi, Etheta, k) will be + parallel-transported along the ray. + Ephi and Etheta are glued to the observer's screen orientation. + Also added a corresponding test in python/test/core.py + to check the orthonormality of the polarization basis. + +commit 73e2580761ac70a772463a6ec22b8c4de6fb0193 +Author: Thibaut Paumard +Date: Wed Jan 24 11:33:12 2018 +0100 + + Implement memory allocation for parallel transport of base vectors + + * add ep0_ ... et3_ tables in Worldline; + * implement logic to allocate/deallocate them depending on parallel_transport_; + * implement parallel transport property, acessibel through e.g. XML. + +commit b942096fac9dac492a6f0ca212ef81143f7e2b97 +Author: Thibaut Paumard +Date: Wed Jan 24 10:54:55 2018 +0100 + + Fixup: add parantheses around the definition of compound quantities GYOTO_QUANTITY_SPECTRAL and GYOTO_QUANTITY_STOKES. + +commit 628026ad297aac81df74d2437452370634e28086 +Author: Thibaut Paumard +Date: Tue Jan 23 17:40:02 2018 +0100 + + Add quantities for Stokes parameters + + Add new constants GYOTO_QUANTITY_SPECTRUM_STOKES_Q, U and V, support them in Scenery and gyoto. + + Still not implemeted in processHitQuantities. + +commit f9ba59e19cf5d8dc581ed65714d5f1ba63878bb4 +Author: Thibaut Paumard +Date: Mon Jan 22 14:38:39 2018 +0100 + + Start implementing parallel transport of local tetrad + commit 02ead454357ccdbc25c6cda8597eaeb139deaf8d Author: Frederic Vincent Date: Wed Dec 20 15:01:15 2017 +0100 diff --git a/NEWS b/NEWS index bca782e9..e674dbae 100644 --- a/NEWS +++ b/NEWS @@ -1,10 +1,24 @@ -NEXT ABI +1.4.0 2019/10/22 ABI 8:0:0 + * Ongoing work to support polarization: + * Support parallel transport of local tetrad + + in many methods, double* coord_ph becomes state_t &coord_ph + + Worldline: new APIs parallelTransport(), eAllocate(), + eDeallocate(), eExpand() + + Scenery: new API parallelTransport() + + Screen: new API getRayTriad() + * Introduce polarized quantities (Inu and alphanu for each + component of the Stokes vector) + + Several methods take additional arguments + + Photon: new method transfer() + + Scenery: new API getSpectralQuantitiesCount() * Support computing Screen tetrad from arbitrary velocity * Metric: + new API zamoVelocity(), default assumes null 3-velocity. - + new method dualOneForm(). + new API observerTetrad(in pos, in fourvel, out screen1,2,3), - default throws an Error. + default throws an Error. + + new methods to help vector computations: dualOneForm(), + normalizeFourVel(), norm(), multiplyFourVect(), + addFourVect(), projectFourVect(), GramSchmidt() * KerrBL: + implement zamoVelocity(); + remove observerTetrad(obskind...); @@ -12,8 +26,22 @@ NEXT ABI * Screen: supports two new obskinds: VelocitySpecified and FullySpecified. obskind is now stored as and int and needs to be specified as e.g. ZAMO in XML. - * Screen: support changing camera orientation (with dangle1 and dangle2) + * Astrobj: + * Generic: new property DeltaMaxInsideRMax + * New kinds FlaredDiskSynchrotron, ThickDisk + * Jet: + + support thermal distribution + + new property JetVphiOverVr + * ThinDisk: new property VelocityKind + * Torus: add a Thermal Synchrotron spectrum + * UniformSphere: deprecate property Alpha + * Screen: + * support changing camera orientation (with dangle1 and dangle2) + * change type of observerkind_ to obskind_t * Python: new module gyoto.animate to help writing Gyoto-based movies. + * Command line: new task mk_video to help writing Gyoto-based movies. + * Register: loadPlugin() now returns a handle o the dlopen'ed plug-in. + * Star/StarTrace: allow expressing dates in proper or coordinate time. 1.3.6 2019/09/27 FTR * Python: improve interface in gyoto.util diff --git a/configure b/configure index b5e24f0a..3172b6df 100755 --- a/configure +++ b/configure @@ -1,6 +1,6 @@ #! /bin/sh # Guess values for system-dependent variables and create Makefiles. -# Generated by GNU Autoconf 2.69 for Gyoto 1.3.6. +# Generated by GNU Autoconf 2.69 for Gyoto 1.4.0. # # Report bugs to . # @@ -590,8 +590,8 @@ MAKEFLAGS= # Identity of this package. PACKAGE_NAME='Gyoto' PACKAGE_TARNAME='gyoto' -PACKAGE_VERSION='1.3.6' -PACKAGE_STRING='Gyoto 1.3.6' +PACKAGE_VERSION='1.4.0' +PACKAGE_STRING='Gyoto 1.4.0' PACKAGE_BUGREPORT='gyoto@sympa.obspm.fr' PACKAGE_URL='' @@ -1488,7 +1488,7 @@ if test "$ac_init_help" = "long"; then # Omit some internal or obsolete options to make the list less imposing. # This message is too long to be a string in the A/UX 3.1 sh. cat <<_ACEOF -\`configure' configures Gyoto 1.3.6 to adapt to many kinds of systems. +\`configure' configures Gyoto 1.4.0 to adapt to many kinds of systems. Usage: $0 [OPTION]... [VAR=VALUE]... @@ -1560,7 +1560,7 @@ fi if test -n "$ac_init_help"; then case $ac_init_help in - short | recursive ) echo "Configuration of Gyoto 1.3.6:";; + short | recursive ) echo "Configuration of Gyoto 1.4.0:";; esac cat <<\_ACEOF @@ -1809,7 +1809,7 @@ fi test -n "$ac_init_help" && exit $ac_status if $ac_init_version; then cat <<\_ACEOF -Gyoto configure 1.3.6 +Gyoto configure 1.4.0 generated by GNU Autoconf 2.69 Copyright (C) 2012 Free Software Foundation, Inc. @@ -2508,7 +2508,7 @@ cat >config.log <<_ACEOF This file contains any messages produced by compilers while running configure, to aid debugging if configure makes a mistake. -It was created by Gyoto $as_me 1.3.6, which was +It was created by Gyoto $as_me 1.4.0, which was generated by GNU Autoconf 2.69. Invocation command line was $ $0 $@ @@ -3484,7 +3484,7 @@ fi # Define the identity of the package. PACKAGE='gyoto' - VERSION='1.3.6' + VERSION='1.4.0' cat >>confdefs.h <<_ACEOF @@ -15741,8 +15741,8 @@ ac_config_headers="$ac_config_headers config.h include/GyotoConfig.h" # http://www.gnu.org/software/libtool/manual/html_node/Updating-version-info.html # C: current; R: revision; A: age # Linux library name is: libgyoto.so.(C-A).A.R -VERSINFO_C=7 -VERSINFO_R=1 +VERSINFO_C=8 +VERSINFO_R=0 VERSINFO_A=0 VERSINFO=${VERSINFO_C}:${VERSINFO_R}:${VERSINFO_A} @@ -21955,7 +21955,7 @@ cat >>$CONFIG_STATUS <<\_ACEOF || ac_write_fail=1 # report actual input values of CONFIG_FILES etc. instead of their # values after options handling. ac_log=" -This file was extended by Gyoto $as_me 1.3.6, which was +This file was extended by Gyoto $as_me 1.4.0, which was generated by GNU Autoconf 2.69. Invocation command line was CONFIG_FILES = $CONFIG_FILES @@ -22021,7 +22021,7 @@ _ACEOF cat >>$CONFIG_STATUS <<_ACEOF || ac_write_fail=1 ac_cs_config="`$as_echo "$ac_configure_args" | sed 's/^ //; s/[\\""\`\$]/\\\\&/g'`" ac_cs_version="\\ -Gyoto config.status 1.3.6 +Gyoto config.status 1.4.0 configured by $0, generated by GNU Autoconf 2.69, with options \\"\$ac_cs_config\\" diff --git a/configure.ac b/configure.ac index 63895f05..359a5f28 100644 --- a/configure.ac +++ b/configure.ac @@ -10,7 +10,7 @@ dnl without any warranty. AC_PREREQ([2.69]) -AC_INIT([Gyoto], [1.3.6], [gyoto@sympa.obspm.fr]) +AC_INIT([Gyoto], [1.4.0], [gyoto@sympa.obspm.fr]) gyoto_test_CXXFLAGS=${CXXFLAGS+set} gyoto_test_CPPFLAGS=${CPPFLAGS+set} AC_CANONICAL_SYSTEM @@ -26,8 +26,8 @@ AC_SUBST([VERBATIM], []) # http://www.gnu.org/software/libtool/manual/html_node/Updating-version-info.html # C: current; R: revision; A: age # Linux library name is: libgyoto.so.(C-A).A.R -VERSINFO_C=7 -VERSINFO_R=1 +VERSINFO_C=8 +VERSINFO_R=0 VERSINFO_A=0 VERSINFO=${VERSINFO_C}:${VERSINFO_R}:${VERSINFO_A} AC_SUBST([VERSINFO])