2019-06-25 11:44:08 -0500 received badge ● Good Answer (source) 2019-06-25 11:44:08 -0500 received badge ● Enlightened (source) 2019-06-18 21:23:27 -0500 received badge ● Popular Question (source) 2019-06-05 19:35:00 -0500 commented answer Bel decomposition Thanks, I understand. 2019-05-22 13:44:30 -0500 commented answer Bel decomposition Thank you very much! I've put BEFORE you code the following (Schwarzschild metric): M = Manifold(4, 'M') print M X. = M.chart(r't r:(0,+oo) th:(0,pi):\theta ph:(0,2*pi):\varphi') X var('m', domain='real') assume(m>=0) g = M.lorentzian_metric('g') g[0,0] = (2*m-r)/r g[1,1] = r/(r-2*m) g[2,2] = r^2 g[3,3] = (r*sin(th))^2 g.display()  I do not understand what is e0 in you code, so I do not know how to proceed. Thanks, Daniel 2019-05-22 03:17:25 -0500 received badge ● Popular Question (source) 2019-05-21 17:01:36 -0500 received badge ● Nice Question (source) 2019-05-18 08:14:37 -0500 asked a question Bel decomposition Hi All, Is there any way to obtain the Bel decomposition (https://en.wikipedia.org/wiki/Bel_dec...) into the electrogravitic, magnetogravitic and topogravitic tensors given the metric? Thanks, Daniel 2019-05-16 06:31:33 -0500 received badge ● Popular Question (source) 2019-04-18 12:49:05 -0500 received badge ● Popular Question (source) 2019-04-14 03:08:34 -0500 commented answer multiprocessing.Pool does not work In Windows it is broken too! 2019-04-11 21:12:18 -0500 marked best answer Differential forms best package I understand there are two different ways to use differential forms in SageMath: 1) in the reference document "manifold.pdf" sage: M = Manifold(4, 'M') sage: a = M.one_form('A') or F = M.diff_form(2, 'F') 2) in the reference document "tensors.pdf" sage: x, y, z = var('x, y, z') sage: U = CoordinatePatch((x, y, z)) sage: F = DifferentialForms(U) sage: form1 = DifferentialForm(F, 0, sin(x*y)); form1 Which is the best? Are they compatible? Can I mix and match from both packages? What are their pros and cons? Thanks. 2019-04-11 21:11:14 -0500 received badge ● Popular Question (source) 2019-03-18 12:49:23 -0500 received badge ● Notable Question (source) 2019-02-27 11:08:30 -0500 commented answer No output when documentation say it should be Thanks, I understand. 2019-02-09 03:16:16 -0500 received badge ● Nice Question (source) 2019-02-09 01:40:47 -0500 commented question No output when documentation say it should be Cont. sage: nab.torsion_form(1)[:] sage: ch_basis = M.automorphism_field() sage: ch_basis[1,1], ch_basis[2,2], ch_basis[3,3] = y, z, x sage: e = M.default_frame().new_frame(ch_basis, 'e') sage: e[:], e[:], e[:] sage: ef = e.coframe() sage: ef[:], ef[:], ef[:] sage: nab.torsion_form(1, e) sage: nab.torsion_form(1, e).comp(e)[:] sage: for i in M.irange(): nab.torsion_form(i, e) == ef[i].exterior_derivative() + \ sum(nab.connection_form(i,j,e).wedge(ef[j]) for j in M.irange())  No output 2019-02-09 01:37:47 -0500 commented question No output when documentation say it should be Sage version is 8.4, Windows Native on a 64bit machine, installed from SageMath-8.4-Installer-v0.4.1.exe I did not compile Sage by myself. I use Jupyter Notebook. I don't use the command line. I did not get any errors. I just got no output at all sage: M = Manifold(3, 'M', start_index=1) sage: c_xyz. = M.chart() sage: nab = M.affine_connection('nabla', r'\nabla') sage: nab[1,1,1], nab[1,1,2], nab[1,1,3] = x*y*z, x^2, -y*z sage: nab[1,2,3], nab[1,3,1], nab[1,3,2] = -x^3, y^2*z, y^2-x^2 sage: nab[2,1,1], nab[2,1,2], nab[2,2,1] = z^2, x*y*z^2, -x^2 sage: nab[2,3,1], nab[2,3,3], nab[3,1,2] = x^2+y^2+z ... 2019-02-08 13:36:30 -0500 asked a question No output when documentation say it should be Hi all, Consider http://doc.sagemath.org/html/en/refer... The very last command of this web page is sage: for i in M.irange(): nab.torsion_form(i, e) == ef[i].exterior_derivative() + \ sum(nab.connection_form(i,j,e).wedge(ef[j]) for j in M.irange())  The output of this command should be: True, True, True, but I get no output at all. In order to duplicate you should include also the previous commands beginning from the torsion_form definition, over all 3 cells. Daniel 2019-01-16 04:04:38 -0500 received badge ● Nice Answer (source) 2019-01-12 17:35:39 -0500 received badge ● Notable Question (source) 2019-01-12 11:00:25 -0500 received badge ● Popular Question (source) 2019-01-12 10:52:38 -0500 answered a question Combine plots with built-in Maxima, trajectory in Sage available? You can combine plots in Maxima with the command "xfun" inside plotdf e.g. plotdf([x+y,x/2+y], [xfun,"-0.7071*x;0.7071*x"])$ In you case: plotdf(4*y/x+x*sqrt(y),[xfun,"x^4*(log(x)+2)^2/4"], [trajectory_at,1,1],[x,0,4],[y,0,10])$  you can use several functions inside the quotes, separated by semi colons. Daniel 2018-12-24 08:37:16 -0500 received badge ● Famous Question (source) 2018-12-16 01:56:13 -0500 commented question Velocity of time I am using Windows 10 on a 64bit machine with SageMath native. Just to clarify, I can't see the 3D plots in the Jupyter nbviewer rendition of the file. If I download the ipynb file to my machine, change the 3D viewer to 'tachyon' then I can see the 3D plots. I'm using Microsoft Edge to open Jupyter notebooks with no problems. As I reported a while ago, the only 3D viewer that works on Windows 10 is 'tachyon'. Daniel 2018-12-16 01:38:32 -0500 commented question Velocity of time Microsoft Edge. 2018-12-16 01:36:54 -0500 marked best answer Wrong Pointer Hi All, When you open: http://nbviewer.jupyter.org/github/sa... and try to download the corresponding ipynb file using the link pointed at by here, it points to SM_Schwarzschild.ipynb not to SM_basic_Schwarzschild.ipynb. Please fix that. Thanks, Daniel 2018-12-14 14:27:23 -0500 commented answer Velocity of time Cont. The big question is Dt_0. Is this something the sensor system should measure, imposed on my by the outside conditions? I understand from your explanation that there is another constraint: the fact that the norm of the tangent vector should be constant to preserve a time-like trajectory, and from this constraint I can calculate Dt_0 using the other measured values. I have several Maxima files in which I calculate numerically the trajectory of a mass in different geometries using rkf45. I always wondered what value should I select for Dt_0. I usually select 1. Thanks, Daniel 2018-12-14 14:21:35 -0500 commented answer Velocity of time Cont. In order to do that, the navigation system needs the initial conditions. I understand t_0=0, this is arbitrary, I understand (r_0, th_0, ph_0), my current radial distance, its polar angle, and the azimuth angle, those are values that should be measured by some kind of on board sensor system, based on information from the outside. I don't know how it is measured but I guess it can be done. I understand (Dr_0, Dth_0, Dph_0), my current radial velocity, and polar and azimuthal angular velocities, I also guess those can be measured too. So (r_0, th_0, ph_0) and (Dr_0, Dth_0, Dph_0) are measured, that is, imposed on me, based on information from the outside and referred to the selected coordinate system. 2018-12-14 14:18:58 -0500 commented answer Velocity of time Hi Karim, Thank you very much for taking the time to answer this question so thoroughly, I really appreciate it. If it is not much of an imposition, I'll recommend including this explanation in the file itself. Sometimes I imagine being in a spaceship, the screen showing a big Black Hole and me wandering if I am going to crash into the Black Hole or I will be able to acquire a stable orbit. I guess the on board navigation system is able to choose the easiest Boyer-Lindquist coordinate system (it is a Schwarzschild Black Hole, no rotation and no charge) and calculate my trajectory. 2018-12-11 08:00:27 -0500 asked a question Velocity of time Hi All, I was looking at http://nbviewer.jupyter.org/github/sa... at the paragraph Computing numerical solutions/Timelike geodesics/Bounded orbit, at the following command: params_values_bounded = {m:1, s_0:0, s_max:1500, t_0:0, r_0:8, th_0:pi/2, ph_0:1e-12, Dt_0:sqrt(80.81)/(4*sqrt(3)), Dr_0:0, Dth_0:0, Dph_0:4.1/64}  I understand there are 4-coordinates (t,r,th,ph), and for each one there is a 2nd order Geodesic equation, so in total should be 8 initial conditions. How did you choose the value of Dt_0? How did you calculate this value? Is this the velocity of time at time zero. What is the physical significance of this value? P.D. Did you notice there are no 3D figures in this Notebook? Thanks, Daniel 2018-12-11 03:19:04 -0500 marked best answer set_name causes an error on external_derivative display I have the following code: R3 = Manifold(3, 'R3', start_index=1, latex_name=r'\mathbb{R}^3') cartesian3d. = R3.chart() omega = R3.diff_form(1, 'omega', latex_name=r'\omega') omega[:] = (x^2+y*z+e^y, y^3+x*z+x*e^y+z^2, e^z+x*y+2*y*z) print omega.category() omega.display() domega = omega.exterior_derivative() domega.set_name(domega, r'\mathrm{d}\omega') print domega.category() domega.display()  I get a TypeError: unsupported operand type(s) for +: 'DiffFormFreeModule_with_category.element_class' and 'str' If I delete the set_name command the code works fine. Daniel