All Questions
10
questions
1
vote
1
answer
261
views
Implementing Preconditioned conjugate gradient
I have implemented the Preconditioned Conjugate Gradient (PCG) method for solving a system of linear equations in Python and I would appreciate it if someone could verify its correctness since I am ...
2
votes
1
answer
286
views
Finite difference estimation of partial derivatives for arbitrary grids in python
I am currently working on testing some codes on python to solve differential equations by finite differences, but some problems and equations i am dealing with are:
Too big to keep manually writing ...
1
vote
2
answers
460
views
Get the best rational approximations for any real number in builtin Python
Best rational approximations is a technical term and not subjective. The best rational approximations are the convergents of a value plus it's semi convergents that are closer to the original value ...
2
votes
1
answer
141
views
Arbitrary precision Euler-Mascheroni constant via a Brent-McMillan algorithm with no math module
Utilizing the below relation, I am able to compute the Euler constant to great precision on a single thread quickly and simply.
My process is to compute the natural log via the AGM(utilizing Pi an ...
1
vote
1
answer
84
views
Euler-Mascheroni Single Thread Speed Improvements
The below code was written to generate γ, for educational purposes.
My general methodology is as follows:
Compute Gamma via the accepted answer's algorithm here.
In order to do this I need to ...
4
votes
1
answer
667
views
Calculate the Euler-Mascheroni constant without the math module
The below code was written to generate γ, for educational purposes.
Single threaded, no functional zeroes required, no binary splitting(which can all be used to compute competitively like y-cruncher, ...
4
votes
2
answers
650
views
Python numerical integration
Could the time complexity of this definite integral algorithm be improved?
...
4
votes
1
answer
4k
views
Compute Gini Coefficient
Recently, I was given a math assignment to calculate Gini Indexes for a table of percent distributions of aggregate income.
The table takes the form of:
...
6
votes
2
answers
10k
views
Newton's method to solve cubic equations
I have used the Newton-Raphson method to solve Cubic equations of the form $$ax^3+bx^2+cx+d=0$$ by first iteratively finding one solution, and then reducing the polynomial to a quadratic $$a1*x^2+b1*x+...
8
votes
1
answer
5k
views
Python class that implements the Newton method
Here is a python function I wrote to implement the Newton method for optimization for the case where you are trying to optimize a function that takes a vector input and gives a scalar output. I use ...