[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-10-31 Thread Imran Ali
I forgot to mention in the SO post (not that I have mentioned plenty already there) that sometimes I get the following error : Traceback (most recent call last): File "tensor.py", line 274, in RC = R.find_Christoffel_tensor() File "tensor.py", line 210, in find_Christoffel_tensor R =

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-10-31 Thread Imran Ali
On Saturday, October 31, 2015 at 6:04:14 PM UTC+1, Imran Ali wrote: > > I just asked the following question at stackoverflow : > > > http://stackoverflow.com/questions/33453941/determining-christoffel-symbols-and-curvature-tensors-for-a-flat-sphere-using-sy > > In essence, the question is two fol

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-05 Thread Imran Ali
If I simply hard code each metric, every thing works properly. For example, take egg carton surface, from sympy.diffgeom import Manifold, Patch, CoordSystem, TensorProduct from sympy import sin,cos dim = 2 m = Manifold("M",dim) patch = Patch("P",m) system = CoordSystem('eg

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-06 Thread Kalevi Suominen
On Thursday, November 5, 2015 at 3:56:09 PM UTC+2, Imran Ali wrote: > > If I simply hard code each metric, every thing works properly. For > example, take egg carton surface, > > from sympy.diffgeom import Manifold, Patch, CoordSystem, TensorProduct > from sympy import sin,cos > dim

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-08 Thread Imran Ali
Hi Kalevi! I got the following output from srepr : Symbol('u') Differential(BaseScalarField(CoordSystem(Symbol('egg_carton'), Patch(Symbol('P'), Manifold(Symbol('M'), Integer(2))), Tuple(Symbol('u'), Symbol('v'))), Integer(0))) I see now what the error is! My code used u as a sympy expression,

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-08 Thread Imran Ali
I was a bit quick there. srepr(u) gives the following output : BaseScalarField(CoordSystem(Symbol('egg_carton'), Patch(Symbol('P'), Manifold(Symbol('M'), Integer(2))), Tuple(Symbol('u'), Symbol('v'))), Integer(0)) On Sunday, November 8, 2015 at 4:40:45 PM UTC+1, Imran Ali wrote: > > Hi Kalevi!

[sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-08 Thread Kalevi Suominen
On Sunday, November 8, 2015 at 5:45:44 PM UTC+2, Imran Ali wrote: > > I was a bit quick there. srepr(u) gives the following output : > > BaseScalarField(CoordSystem(Symbol('egg_carton'), Patch(Symbol('P'), > Manifold(Symbol('M'), Integer(2))), Tuple(Symbol('u'), Symbol('v'))), > Integer(0)) > >

Re: [sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-06 Thread Aaron Meurer
What you should do is post the complete code to stackoverflow, as you have done here. Aaron Meurer On Fri, Nov 6, 2015 at 10:51 AM, Kalevi Suominen wrote: > > > On Thursday, November 5, 2015 at 3:56:09 PM UTC+2, Imran Ali wrote: >> >> If I simply hard code each metric, every thing works properly

Re: [sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-06 Thread Francesco Bonazzi
I answered to you on StackOverflow. You should be more clear by posting questions, this looks like a problem in your calculation logic, not in SymPy. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving em

Re: [sympy] Re: Determining Christoffel Symbols and Curvature Tensors for a Flat Sphere using SymPy.Diffgeom

2015-11-08 Thread Imran Ali
Hey Aaron! Thanks for the tip, I posted the code on SO. But, I think that Kalevi solved the problem. I found the issue. On Friday, November 6, 2015 at 5:54:36 PM UTC+1, Aaron Meurer wrote: > > What you should do is post the complete code to stackoverflow, as you > have done here. > > Aaron Meu