Macaulay2 » Documentation
Packages » NumericalSchubertCalculus :: solutionsToAffineCoords
next | previous | forward | backward | up | index | toc

solutionsToAffineCoords -- writes solutions in global coords to affine coordinates.

Description

Takes a list of solutions $s\in Gr(k,n)$ in global coordinates and writes them as solutions in coordinates of the affine patch that has an identity matrix in the last $k$ rows.

i1 : Pblm = {
         ({2,1}, random(RR^6,RR^6)),
         ({2,1}, random(RR^6,RR^6)),
         ({2,1}, random(RR^6,RR^6))
         }

o1 = {({2, 1}, | .892712  .89189  .0741835 .741046 .830833 .191734  |), ({2,
               | .673395  .231053 .808694  .108386 .538554 .403215  |       
               | .29398   .461944 .362835  .348931 .873665 .615911  |       
               | .632944  .775187 .706096  .562428 .415912 .0147867 |       
               | .0258884 .909047 .127435  .246268 .606588 .223028  |       
               | .714827  .314897 .254482  .153346 .848005 .388829  |       
     ------------------------------------------------------------------------
     1}, | .557119 .96518  .0647412 .174853 .444183 .184779 |), ({2, 1}, |
         | .873708 .681683 .877846  .626892 .644366 .370833 |            |
         | .7037   .914199 .0340514 .350611 .194945 .339222 |            |
         | .681869 .887381 .507989  .379495 .518585 .062212 |            |
         | .276259 .169813 .150294  .237252 .987173 .465736 |            |
         | .605659 .965004 .656391  .116721 .568273 .40273  |            |
     ------------------------------------------------------------------------
     .164647 .562839  .0645275 .501243 .205375  .0821679 |)}
     .713493 .629991  .283709  .154289 .276652  .10394   |
     .909537 .479826  .977573  .467203 .0958269 .280679  |
     .566034 .815167  .212436  .765564 .605398  .475179  |
     .305423 .97723   .592747  .305946 .883549  .130004  |
     .732358 .0595849 .831802  .53632  .942865  .522285  |

o1 : List
i2 : S = solveSchubertProblem(Pblm, 3,6)

o2 = {| .362674-1.92986e-16ii  5.87881-7.96063e-14ii .067722+1.82633e-16ii |,
      | -.168148-1.45574e-16ii 10.2268-1.45803e-13ii .359604+2.40887e-16ii | 
      | .287667-6.35526e-17ii  5.41301-7.40938e-14ii .178413+2.13998e-16ii | 
      | .399967-1.3683e-16ii   10.4365-1.44272e-13ii .172408+2.07443e-16ii | 
      | .893699                2.90938-4.09052e-14ii .127632+7.36638e-17ii | 
      | -.108865-1.54531e-16ii 5.44614-7.8727e-14ii  .099289+1.93934e-16ii | 
     ------------------------------------------------------------------------
     | .675164+8.15266e-16ii  .770431+1.31064e-15ii -.0909566-3.82716e-16ii 
     | .0675706+6.14976e-16ii .870501+2.4005e-15ii  .150312-5.04789e-16ii   
     | .390574+2.68476e-16ii  .658364+1.21988e-15ii -.00751668-4.48442e-16ii
     | .621526+5.78034e-16ii  1.17845+2.3753e-15ii  -.0078269-4.34707e-16ii 
     | .902762                .284467+6.73463e-16ii .0636297-1.54366e-16ii  
     | .141356+6.52813e-16ii  .394186+1.29616e-15ii -.0692081-4.06397e-16ii 
     ------------------------------------------------------------------------
     |}
     |
     |
     |
     |
     |

o2 : List
i3 : solutionsToAffineCoords S

o3 = {| 1.17594  -.247384 -1.04186 |, | -.0717886 .515973  1.79675  |}
      | -2.59821 1.69933  5.94898  |  | 1.97727   -.796664 -3.12795 |
      | -.743029 .890962  1.94183  |  | .504284   .067644  .113771  |
      | 1        0        0        |  | 1         0        0        |
      | 0        1        0        |  | 0         1        0        |
      | 0        0        1        |  | 0         0        1        |

o3 : List

Caveat

This function may fail if the solutions are not in general position (if they cannot fit the specific local coordinates) One way to avoid this is by applying a random linear transformation to the solutions before calling this function

Ways to use solutionsToAffineCoords:

  • solutionsToAffineCoords(List)

For the programmer

The object solutionsToAffineCoords is a method function.


The source of this document is in NumericalSchubertCalculus/doc.m2:474:0.