Changes between Initial Version and Version 2 of Ticket #30946


Ignore:
Timestamp:
11/22/20 20:34:11 (10 months ago)
Author:
jipilab
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • Ticket #30946

    • Property Status changed from new to needs_review
  • Ticket #30946 – Description

    initial v2  
    44
    55{{{
    6 sage: A=[[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1],[1/4,1/4,1/4,1/4]]                                                                                                                                               
    7 sage: n=len(A)                                                                                                                                                                                                     
    8 sage: A=[vector(v) for v in A]                                                                                                                                                                                     
    9 sage: AP = Polyhedron(vertices=A)                                                                                                                                                                                 
    10 sage: M,b=AP.affine_hull_projection(orthonormal=True,extend=True,as_affine_map=True)                                                                                                                               
     6sage: A=[[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1],[1/4,1/4,1/4,1/4]]
     7sage: n=len(A)
     8sage: A=[vector(v) for v in A]
     9sage: AP = Polyhedron(vertices=A)
     10sage: M,b=AP.affine_hull_projection(orthonormal=True,extend=True,as_affine_map=True)
    1111sage: V=[]
    12 ....: for i in range(n):
    13 ....:     for j in range(i+1,n):
    14 ....:             V.append((A[i]-A[j])/2)
    15 ....:                                                                                                                                                                                                             
    16 sage: Z=polytopes.zonotope(V)                                                                                                                                                                                     
    17 sage: T = M.matrix().transpose()                                                                                                                                                                                   
    18 sage: timeit('T*Z')                                                                                                                                                                                               
     12....: for i in range(n):
     13....:     for j in range(i+1,n):
     14....:             V.append((A[i]-A[j])/2)
     15sage: Z=polytopes.zonotope(V)
     16sage: T = M.matrix().transpose()
     17sage: timeit('T*Z')
    19185 loops, best of 3: 78.5 ms per loop
    2019}}}
     
    2322
    2423{{{
    25 sage: A=[[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1],[1/4,1/4,1/4,1/4]]                                                                                                                                               
    26 sage: n=len(A)                                                                                                                                                                                                     
    27 sage: A=[vector(v) for v in A]                                                                                                                                                                                     
    28 sage: AP = Polyhedron(vertices=A)                                                                                                                                                                                 
    29 sage: M,b=AP.affine_hull_projection(orthonormal=True,extend=True,as_affine_map=True,minimal=True)                                                                                                                 
    30 sage: V=[] 
     24sage: A=[[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1],[1/4,1/4,1/4,1/4]]
     25sage: n=len(A)
     26sage: A=[vector(v) for v in A]
     27sage: AP = Polyhedron(vertices=A)
     28sage: M,b=AP.affine_hull_projection(orthonormal=True,extend=True,as_affine_map=True,minimal=True)
     29sage: V=[]
    3130....: for i in range(n):
    3231....:     for j in range(i+1,n):
    3332....:             V.append((A[i]-A[j])/2)
    34 ....:                                                                                                                                                                                                             
    35 sage: Z=polytopes.zonotope(V)                                                                                                                                                                                     
    36 sage: T = M.matrix().transpose()                                                                                                                                                                                   
    37 sage: timeit('T*Z')                                                                                                                                                                                               
     33sage: Z=polytopes.zonotope(V)
     34sage: T = M.matrix().transpose()
     35sage: timeit('T*Z')
    383625 loops, best of 3: 18 ms per loop
    3937}}}