Ask Your Question

How to create a loop with two variables?

asked 2015-01-03 05:55:13 -0500

Kevin Smith gravatar image

updated 2015-01-03 06:03:40 -0500

tmonteil gravatar image

I have defined the following functions:

def I(x):
    return 1

def Moebius(x):
    return sum(moebius(i) * I(x*i^(-1)) for i in range(1,x+1))

The output is of course Mertens function. How do I turn Moebius(x) into an operator, so that I can define Moebius^2(x) as the Moebius transform of the output?

edit retag flag offensive close merge delete

1 answer

Sort by ยป oldest newest most voted

answered 2015-01-03 06:29:08 -0500

tmonteil gravatar image

updated 2015-01-03 06:33:06 -0500

First, the * I(x*i^(-1)) part in your Moebius function will only add some useless computation since it is a multiplication by 1.

Then, you can define the composition of your Moebius function by itself as follows:

M2 = lambda x : Moebius(Moebius(x))

If you want the ^2 notation to work, you have to define a Python class since ^2 calls the __pow__ method of your object, this will require some more work around the following lines:

class MyOperator():
    def __init__(self, func):
        self.func = func
    def __pow__(self, exponent):
        def iterfunc(x):
            result = x
            for i in range(exponent):
                result = self.func(result)
            return result
        return MyOperator(iterfunc)
    def __call__(self, x):
        return self.func(x)

Then you can do:

sage: f = MyOperator(Moebius)
sage: f(10)
sage: g = f^2
sage: g(10)
sage: f(f(10))

Or similarly:

sage: inc = lambda x : x+1
sage: f = MyOperator(inc)
sage: f(14)
sage: g = f^12
sage: g(14)
sage: h = f^0
sage: h(100)
edit flag offensive delete link more


Thank you very much for such a prompt and clear reply. The I function was there because I'd hoped to have it as non-trivial input to the function.

Kevin Smith gravatar imageKevin Smith ( 2015-01-03 06:51:31 -0500 )edit

Without the presence of a function such as I(x), I don't see how Moebius even has an input ( other than x). Can you explain how your solution had another input please? Thank you

Kevin Smith gravatar imageKevin Smith ( 2015-01-03 08:37:26 -0500 )edit

You could just write

def Moebius(x):
    return sum(moebius(i) for i in range(1,x+1))

Which other input than x would you like ?

tmonteil gravatar imagetmonteil ( 2015-01-06 01:07:12 -0500 )edit

Excuse my absolute ignorance of programming- though I am now certain this is complete nonsense from sage's point of view, in an algebraic sense, I want the function I(x) to be something in the domain of the operator that I name Moebius. The reason I set I(x) to be constant initially was merely to demonstrate that I get the right output in this case. What I really want to do is take the output as an input to further applications of Moebius, so that I get the higher order Mertens functions, that is, the Mellin transforms of the negative powers of the Riemann zeta function. Of course there ought to be another way to do this involving extracting the number of times each primes divides a given integer, but that looked rather more difficult at first. It's equally likely that I am mistaken about that too though.

Kevin Smith gravatar imageKevin Smith ( 2015-01-18 04:56:44 -0500 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower


Asked: 2015-01-03 05:55:13 -0500

Seen: 187 times

Last updated: Jan 03 '15