Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • Aug 05 12:45

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    included exp forms for bessel f… (compare)

  • Aug 05 12:05

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Minor change (compare)

  • Aug 05 11:59

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Refactored mrv_leadterm in grun… (compare)

  • Aug 01 03:22

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Added e.is_negative case for be… (compare)

  • Aug 01 02:43

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    hanled e.is_negative case for b… (compare)

  • Aug 01 02:03

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Replacing exp rewrite with trig… (compare)

  • Jul 31 11:41

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Minor changes (compare)

  • Jul 31 11:30

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    added support when exponent of … (compare)

  • Jul 31 11:00

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Fixed leading term method for b… (compare)

  • Jul 28 04:15

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Added tests for 7224 (compare)

  • Jul 28 03:19

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Added tests for series of besse… (compare)

  • Jul 28 02:51

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Added series test for besseli f… (compare)

  • Jul 28 02:38

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Added tests for leading term me… (compare)

  • Jul 28 02:25

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Added nseries method for bessel… (compare)

  • Jul 27 18:26

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Fixed failing tests (compare)

  • Jul 27 14:49

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    Added and fixed some tests (compare)

  • Jul 27 14:20

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    implemented leading term and ns… (compare)

  • Jul 27 14:16

    anutosh491 on GSoC_Pr4.1_Implementing_few_series_methods_for_bessel_functions

    fixed series for bessely (compare)

  • Jul 26 08:30

    anutosh491 on GSOC_PR3_Implementing_nseries_for_some_inverse_functions

    (compare)

  • Jul 25 16:17

    anutosh491 on GSOC_PR3_Implementing_nseries_for_some_inverse_functions

    Refactored series method in exp… (compare)

Micah
@micahscopes
@jksuom currently I've just been making a list of symbolic multivariate polynomial expressions. The symbols are the components of the two input vectors, E. G. f([u_1, u_2,... , u_n], [v_1, v_2,... , v_n]) would return a list of polynomials of terms from the u and v lists.
In Sage, you can construct algebras over a SymbolicRing, then take the product of two symbolic vectors and convert the result to a list of sympy expressions
(after that I'm using the glsl code printer in sympy)
Micah
@micahscopes
Ah... I found a clue! https://math.stackexchange.com/questions/2889140/can-a-rank-three-tensor-act-as-a-trilinear-bilinear-or-linear-map
Kalevi Suominen
@jksuom
A bilinear map like f:M×M→Mf: M\times M \to M corresponds naturally to a rank three tensor.
A (1, 2) -tensor actually, I think.
Micah
@micahscopes
@jksuom that makes sense... so if I figure out how to represent my algebras as (1, 2)-tensors, I should also be able to find a way, using tensor algebra, to more easily compose them
Kalevi Suominen
@jksuom
I think that is true though I'm not sure what you mean by 'compose'.
Micah
@micahscopes
I'm not sure either, I don't think "compose" is the mathematically correct term, but I'd like to do multiple "base changes"
I'd like to be able to construct some of these division algebras: https://core.ac.uk/download/pdf/25253839.pdf
Kalevi Suominen
@jksuom
Those tensor products of algebras can be represented by products of tensors (taking proper care of the indices).
Ayushman Koul
@ayushmankoul
@jksuom As you pointed if I am not wrong we need to amend subclasses of BaseSeries containing 2D lines series.Please can you tell me how should I proceed further on this issue ?
Kalevi Suominen
@jksuom
You could probably start by adding methods that would create instances of List2DSeries to some geometric classes such as Segment and Triangle, and also Polygon. (Those could then be appended to a suitable Plot object for showing.)
Micah
@micahscopes
@jksuom I think I'm starting to get some intuition for this idea
Kalevi Suominen
@jksuom
The components aijka_{ijk} of the tensor associated with the bilinear product in an algebra can be found as follows. If (ei)(e_i) is a basis of the algebra, then the components are the coefficients of the pairwise products eiej=∑kaijkeke_i e_j = \sum_k a_{ijk} e_k.
Kalevi Suominen
@jksuom
For example, if C\Bbb{C} is considered as an R\Bbb{R}-algebra with basis (1,i)(1, i), then a110=−1a_{110}=-1 is the coefficient of 11 in iiii and a111=0a_{111}=0.
Micah
@micahscopes
@jksuom !thank you so much!
basically then, to construct the tensor of a finite k-algebra, I just need to get the coefficients in k of pairwise products of basis elements of the algebra
Micah
@micahscopes
I think I need to really play with the summation notation
(and the concept of summation/contraction)
Kalevi Suominen
@jksuom
It seems to me that tensor.array might be the best module for handling algebras and their tensor products.
Micah
@micahscopes
yeah, that's my intuition too. I didn't realize how malleable tensor contraction could be... basically it's very important to keep track of the indices
the indices aren't inherently ordered, so it's gonna be up to me to keep them structured
I had been getting caught up on thinking about duality and covariance/contravariance, but for this problem it's maybe not so necessary to worry about that aspect for now
Micah
@micahscopes
by the way... one of my motivations for this is to be able to do automatic differentiation with "dual numbers". have you heard of this?
https://matiasmorant.wordpress.com/2017/12/29/dual-numbers/
Kalevi Suominen
@jksuom
It is a standard construction in algebra. To any ring AA, one can associate the "dual ring" A[δ]A[\delta] where δ2=0\delta^2 = 0.
Micah
@micahscopes
awesome
I didn't realize that
what is that square bracket notation?
is that a "ring extension"?
Kalevi Suominen
@jksuom
Yes, extension of AA by δ\delta.
Its elements are polynomials in δ\delta with coefficients in AA. In this case, only linear polynomials.
Micah
@micahscopes
that makes sense
you can use the dual ring to do automatic differentiation
in my case, I'd like to use it to render fractals made out of clifford algebras and other interesting algebras
with automatic differentiation, you can create a "distance estimator" that can be used to render the surface of a fractal set
(theoretically)
so far I've only been able to render these fractals using a "brute force" approach, which actually goes and tests each point to see if it's in the fractal
this doesn't scale well for large dimensional algebras
http://blog.hvidtfeldts.net/index.php/2011/12/distance-estimated-3d-fractals-vii-dual-numbers/
Abhigyan Khaund
@abhigyank
Hey, I am trying to plot a scatter of points using the Plot module in sympy. From the documention I am unable to figure out how do I plot a scatter of points as we do in matplotlib.
I have 2 arrays; call them A and B. I can do a scatter plot using matplotlib like this:
plt.plot(A,B)
How do I do the same in Sympy?
Aaron Meurer
@asmeurer
Sympy only deals with plotting of sympy expressions. If you want to plot a dataset, you should use matplotlib, just as you have already done
Abhigyan Khaund
@abhigyank
Ok thanks!
Is there a way to plot the dataset in matplot lib and plot the sympy expression using sympy and somehow display the two plots in a single window?
Aaron Meurer
@asmeurer
I think you might be able to access the figure from the sympy plot. I'm not sure.
Another way is to lambdify the sympy expression and plot it manually with np.linspace
udit4
@udit4
Hello all, i am new to the SymPy project.. Can anyone help me where can i find complete documentation of SymPy and some other relevant sources, so that i can get some more insights about project and then start contributing to the same. Thanks in advance.