conversion of function from Matlab to python

Question

I have a MATLAB function :

Bits=30
NBits= ceil(fzero(@(x)2^(x) - x -1 - Bits, max(log2(Bits),1)))

I want to convert it to python, I wrote something like this so far:

from numpy import log, log2
from scipy.optimize import root_scalar

def func(x,Bits):
    return ((x)2^(x)-x-1-Bits, max(log2(Bits)))

However it says that it need to be (x)*2^ Does anybody know first, if the conversion from Matlab to python is correct? and second if * has to be added?

Upon suggestion I wrote this lambda function:

lambda x: (2^(x) -x -1 -Bits) , max(log2(Bits))

but I get this error: TypeError: 'numpy.float64' object is not iterable

The python syntax for "to the power of" is ** rather than ^. So to raise something to the power of 2 you would do x**2. — scotty3785
– scotty3785, Commented Dec 16, 2020 at 14:21
Note that the @(x) in your MATLAB code indicates an anonymous function. In Python this can be created using lambda — scotty3785
– scotty3785, Commented Dec 16, 2020 at 14:24
oh right, so I just need to use the lambda so I apply it to the second argument in this case — Paolo Lorenzini
– Paolo Lorenzini, Commented Dec 16, 2020 at 14:26
I think there was some mistake in the original code, there should be a space — Paolo Lorenzini
– Paolo Lorenzini, Commented Dec 16, 2020 at 14:27
ceil(fzero(@(x) 2^(x) - x -1 - numDataBits, max(log2(numDataBits),1))) — Paolo Lorenzini
– Paolo Lorenzini, Commented Dec 16, 2020 at 14:27

scotty3785 · Accepted Answer · 2020-12-16 14:44:37Z

2

I don't have numpy or scipy on this computer so here is my best attempt at an answer.

def YourFunc(Bits):
    return math.ceil(root_scalar(lambda x: (2**x)-x-1-Bits, x0 = max(log2(Bits),1)))

Bits = 30
NBits = YourFunc(30)

print(NBits)

I used this function for log2 rather than the one from numpy. Try it

def log2(x):
    return math.log(x,2)

edited Dec 16, 2020 at 14:44

answered Dec 16, 2020 at 14:39

scotty3785

7,0711 gold badge29 silver badges38 bronze badges

Sign up to request clarification or add additional context in comments.

9 Comments

Paolo Lorenzini Over a year ago

Traceback (most recent call last): File "<stdin>", line 1, in <module> File "<stdin>", line 2, in YourFunc File "/Library/Frameworks/Python.framework/Versions/3.9/lib/python3.9/site-packages/scipy/optimize/_root_scalar.py", line 230, in root_scalar raise ValueError('Unable to select a solver as neither bracket ' ValueError: Unable to select a solver as neither bracket nor starting point provided.

Paolo Lorenzini Over a year ago

I think we are on the right path but I got this error now after importing root_scalar

Paolo Lorenzini Over a year ago

Traceback (most recent call last): File "<stdin>", line 1, in <module> File "<stdin>", line 2, in YourFunc File "/Library/Frameworks/Python.framework/Versions/3.9/lib/python3.9/site-packages/scipy/optimize/_root_scalar.py", line 251, in root_scalar raise ValueError('x1 must not be None for %s' % method) ValueError: x1 must not be None for secant

Paolo Lorenzini Over a year ago

sorry this now after I used the edited function with only one argument, without x, yeah no need to specify x

scotty3785 Over a year ago

@PaoloLorenzini I don't know how the root_scalar function works I'm afraid. I'd try either adding x1=1 as a third argument to root_scalar or changing the solver method method='bisect' or one of the other methods list in the documentation.

|

Collectives™ on Stack Overflow

conversion of function from Matlab to python

1 Answer 1

9 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

9 Comments

Your Answer

Sign up or log in

Post as a guest

Related