What is the cleanest (most "pythonic") way of printing every unique combination of three items in the following list? strats = ["1","2","3","4"] The solution needs to avoid duplicates: 1 and 1 and 1 No (not unique) 1 and 1 and 2 No (not unique) 1 and 1 and 3 No (not unique) … 1 and… Read More Python print every unique combination of list items
I would like to transfer some of variables calculated in Matlab (R2022a) to Python environment (3.10). I have figured out that the SciPy package has function loadmat() and I may save vars in mat file and then I should read the file with Python. However, I’m not able to use the function. import scipy scipy.loadmat(‘myVars.mat’)… Read More How to transfer variables from Matlab to Python
I am new to Python, and am trying to solve this differential equation using scipy.interpolate.odeint. However, I keep getting a TypeError. I have looked, and cannot find how to fix this issue to get the odeint module to work. Below is my code: import pandas as pd import matplotlib.pyplot as plt import numpy as np… Read More Scipy.integrate.odeint TypeError: Float object is not callable
I have two functions defined as integrals where the lower and upper limits are (0-1) and (0-10), respectively. I found out the scipy library and was trying to use it for this purpose. Below I am showing my functions. f(x) = (p*x) dx g(x) = alpha*(ln(f(x))) dx I am trying to calculate the value of… Read More how to integrate a composite function using scipy
I would like to know where I made a mistake in the following code. I suspect it is some basic python mistake that has nothing to do with scipy. I am trying to pass a function to optimize to scipy.optimize.differential_evolution. def func_to_opt(x, TRANS_MIN_BV=TRANS_MIN_BV, SUBS_VAL=100, model=model): """Returns Rsp if BV is above TRANS_MIN_BV, SUBS_VAL if not.… Read More scipy.optimize.differential_evolution – unable to pass a function to optimize into it
I want to calculate complex eigenvectors (Psi functions for shrodinger equation) for 20k*20k matrices, but always gets smth like this: [ 2.99009782e-09 -1.12381299e-08 -2.82346868e-08 …, 6.20967928e-34 -4.80675528e-34 -3.84848719e-35] [ 4.07337553e-08 -1.45976681e-07 -3.47961439e-07 …, 7.43558322e-34 -5.74815572e-34 -4.61317607e-35] [ 5.51921102e-07 -1.88491289e-06 -4.26000711e-06 …, 9.29535407e-34 -7.15304083e-34 -5.78846547e-35] As I understand it just replaces part of data with dots,… Read More eigenvectors created by numpy.linalg.eig have dots instead of data
I have a problem with my code. It looks different everytime i run it. Any ideas? I don’t see any problem. I am looking at this code since 2h and I can’t find the problem… import numpy as np import matplotlib.pyplot as plt from scipy import stats x = np.arange(0,24, 1) y = stats.poisson.pmf(x, mu=13)… Read More Plot looks different everytime i run the code
I’m currently working on computing various similarity metrics between vectors such as cosine similarity, euclidean distance, mahalanobis distance, etc. As I’m working with vectors that can be very large, I need compute time to be minimal. I’m struggling to understand how to work with vectors of different dimensions (they, do, however, share one dimension) and… Read More Understanding broadcasting and arithmetic operations on different dimension tensors
I was comparing the result of my following python calculation with Mathematica: https://www.wolframalpha.com/input?i=sum+%28500+choose+r+%29%28-1%29%5Er+%2F%28r%21%29+%2C+r%3D0+to+500 import numpy as np from decimal import * import scipy.special from scipy.special import factorial getcontext().prec = 30 i = 500 sum(np.array([scipy.special.comb(Decimal(i), (r), exact=True)*pow(-1, r)/Decimal(factorial(r, exact=False)) for r in range(i+1)])) Both calculations are giving almost same value until i = 400 but failing… Read More How to perform binomial-coefficient and factorial calculation with more precision?
How can I write import foo.bar in __init__.py so it will load the system-wide version of foo/bar.py when run from most places, but will load the local version of bar.py when run from within foo‘s source directory? # foo/__init__.py from foo.bar import baz baz() # foo/bar.py def baz: print(‘Hello") This will always load the site-wide… Read More Python imports in deployed/local packages