How should I write this algorithm?

AlgebraicGeometry

asked 2013-04-14 08:16:08 +0100

Neda
65 ●5 ●10 ●20

updated 2013-04-14 08:26:35 +0100

hello I want to write this algorithm in sage, but I don't know how to write it as an algorithm in sage, could you please help me

Algorithm : update basis-update of intermediate basis G with reduct h Given: a finite set G_old ? B[X] and the reduct 0!=h ? B[X]

Find: updates G_new ? B[X] of G_old

begin

1: G_new <-- Ø

2: while G_old !=Ø do

3: selecting from G_old ; G_old <-- G_old {g}

4: if HT(h) not divisible HT(g) then

5: G_new <-- G_new U {g}

6: end

7:end

8: G_new <-- G_new U {h}

9: return G_new

end

thank you so much

edit retag flag offensive close merge delete

add a comment

1 Answer

Sort by » oldest newest most voted

answered 2013-04-15 07:35:34 +0100

slelievre
17684 ●22 ●160 ●349 http://carva.org/samue...

Some of the code you wrote didn't make it. You should indent code so that it is displayed properly.

I don't exactly understand your algorithm, but from what I can grasp it might be something like this.

def update_basis(G,h):
    G_old = G.copy()
    G_new = []
    while G_old:
        g = G_old.pop()
        if ht(h) % ht(g):
            G_new.append(g)
    G_new.append(h)
    return G_new

edit flag offensive delete link

Comments

thank you so so much, I think I should learn python programming ..

Neda ( 2013-04-15 13:37:59 +0100 )edit

add a comment

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

Stats

Asked: 2013-04-14 08:16:08 +0100

Seen: 373 times

Last updated: Apr 15 '13

How should I write this algorithm? edit