|\^/| Maple V Release 4 (University of Bath) ._|\| |/|_. Copyright (c) 1981-1996 by Waterloo Maple Inc. All rights \ MAPLE / reserved. Maple and Maple V are registered trademarks of <____ ____> Waterloo Maple Inc. | Type ? for help. # g2direct-sample # Maple code example of a direct computation of the hyperKaehler # potential for the next-to-minimal nilpotent orbit in G2 # by # Piotr Kobak and Andrew Swann # # This code is generated from a noweb source file g2.nw # See that for further description and comments. # RCS info from g2.nw: # $Id: mapleoutput-g2direct-sample,v 1.1 2000/01/05 14:11:32 swann Exp $ > restart; > MetricNormalisation:=k^2; 2 MetricNormalisation := k > read `g2nmin-direct`; Warning, new definition for norm Warning, new definition for trace Our base point is [ 1/2 ] [0 s 0 t 2 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] X := [ 1/2] [0 0 0 0 0 0 -t 2 ] [ ] [0 -t 0 0 0 0 0 ] [ ] [0 0 t 0 0 0 -s ] [ ] [0 0 0 0 0 0 0 ] Defining J bytes used=1000068, alloc=786288, time=0.78 Computing change of variables 2 2 2 eta1 := 2 k (s + 3 t ) 2 4 2 2 4 eta2 := 4 k (s + 6 s t + 3 t ) bytes used=2000484, alloc=1179432, time=1.76 bytes used=3000676, alloc=1441528, time=2.69 > A := matrix(3,3,0): > A[2,3] := 1/s: > print(A); [0 0 0 ] [ ] [0 0 1/s] [ ] [0 0 0 ] > xA := so7Lb(X,sl3(A)); [0 0 1 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] xA := [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 -1] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] > A1 := J2stP(xA,1); bytes used=4000988, alloc=1965720, time=4.02 bytes used=5001208, alloc=1965720, time=5.47 [ 2 ] [ rhos t rhot rhos ] [0 , 0 , - 1/4 ----- - 1/4 ----------- + 1 , 0 , 0 , 0 , 0] [ 4 4 ] [ k k s ] [ ] [0 , 0 , 0 , 0 , 0 , 0 , 0] [ ] [0 , 0 , 0 , 0 , 0 , 0 , 0] [ ] A1 := [0 , 0 , 0 , 0 , 0 , 0 , 0] [ ] [ 2 ] [ rhos t rhot rhos ] [0 , 0 , 0 , 0 , 0 , 0 , 1/4 ----- + 1/4 ----------- - 1] [ 4 4 ] [ k k s ] [ ] [0 , 0 , 0 , 0 , 0 , 0 , 0] [ ] [0 , 0 , 0 , 0 , 0 , 0 , 0] > print(`First equation`); First equation > e1 := -numer(A1[1,3]); 2 4 e1 := rhos s + t rhot rhos - 4 k s > A := matrix(3,3,0): > A[1,1]:=1: > A[3,3]:=-1: > print(A); [1 0 0] [ ] [0 0 0] [ ] [0 0 -1] > xA := so7Lb(X,sl3(A)); [ 1/2 ] [0 -s 0 -t 2 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] xA := [ 1/2] [0 0 0 0 0 0 t 2 ] [ ] [0 t 0 0 0 0 0 ] [ ] [0 0 -t 0 0 0 s ] [ ] [0 0 0 0 0 0 0 ] > A2 := J2stP(xA,1); bytes used=6001516, alloc=1965720, time=7.01 bytes used=7001792, alloc=1965720, time=8.44 bytes used=8002112, alloc=1965720, time=9.93 bytes used=9110140, alloc=2162292, time=14.47 A2 := [ 2 2 [ s rhos t s rhot rhoss t rhot rhost t rhot rhos [0 , 1/4 ------- + 1/8 -------------- + 1/8 ------------- + 1/4 ----------- [ 4 4 4 4 [ k k k k 2 2 s rhos rhoss t s rhos rhost t rhos rhott + 1/4 ------------- + 3/8 -------------- - s + 1/8 ------------- , 0 , 4 4 4 k k k 2 1/2 2 1/2 1/2 1/2 t rhot 2 t rhos 2 t s rhos 2 rhoss -t 2 + 1/36 ------------ + 1/4 ------------ + 1/4 ------------------- 4 4 4 k k k 2 1/2 1/2 t rhos 2 rhost t s rhot 2 rhost + 1/4 ------------------ + 1/36 ------------------- 4 4 k k 2 1/2 2 2 t rhot 2 rhott s rhot t rhos rhoss + 1/36 ------------------ , 0 , 1/144 ------- - 3/16 ------------- 4 4 4 k k k 2 t s rhot rhoss t rhot rhost t rhot rhos - 1/48 -------------- - 1/48 ------------- - 1/24 ----------- 4 4 4 k k k 2 3 t s rhos rhost s rhot rhost t rhos rhost - 1/48 -------------- + 1/144 ------------- - 3/16 ------------- 4 4 4 k k k s 2 2 2 ] t rhos t s rhot rhott t rhos rhott ] - 3/16 -------- + 1/144 -------------- - 1/48 ------------- , 0] 4 4 4 ] k s k k ] [ 2 2 [ s rhot t rhos rhoss [0 , 0 , 0 , 0 , 0 , 0 , - 1/144 ------- + 3/16 ------------- [ 4 4 [ k k 2 t s rhot rhoss t rhot rhost t rhot rhos + 1/48 -------------- + 1/48 ------------- + 1/24 ----------- 4 4 4 k k k 2 3 t s rhos rhost s rhot rhost t rhos rhost + 1/48 -------------- - 1/144 ------------- + 3/16 ------------- 4 4 4 k k k s 2 2 2 ] t rhos t s rhot rhott t rhos rhott] + 3/16 -------- - 1/144 -------------- + 1/48 -------------] 4 4 4 ] k s k k ] [0 , 0 , 0 , 0 , 0 , 0 , 0] [ 2 1/2 1/2 [ t rhot 2 rhost t s rhos 2 rhost [0 , 0 , 1/48 ------------------ + 1/48 ------------------- [ 4 4 [ k k 1/2 1/2 t rhot 2 rhos t s rhot 2 rhoss + 1/24 ---------------- + 1/48 ------------------- 4 4 k k 2 1/2 2 1/2 s rhot 2 rhost t rhos 2 rhoss - 1/144 ------------------ + 3/16 ------------------ 4 4 k k 3 1/2 2 2 1/2 1/2 t rhos 2 rhost t rhos 2 t s rhot 2 rhott + 3/16 ------------------ + 3/16 ------------- - 1/144 ------------------- 4 4 4 k s k s k 2 1/2 2 1/2 t rhos 2 rhott s rhot 2 1/2 + 1/48 ------------------ - 1/144 ------------ , 0 , 0 , 0 , t 2 4 4 k k 2 1/2 2 1/2 1/2 t rhot 2 t rhos 2 t s rhos 2 rhoss - 1/36 ------------ - 1/4 ------------ - 1/4 ------------------- 4 4 4 k k k 2 1/2 1/2 t rhos 2 rhost t s rhot 2 rhost - 1/4 ------------------ - 1/36 ------------------- 4 4 k k 2 1/2 ] t rhot 2 rhott] - 1/36 ------------------] 4 ] k ] [ 2 2 [ t rhot t s rhos rhoss t rhos [0 , t - 1/36 ------- - 1/4 -------------- - 1/4 ------- [ 4 4 4 [ k k k 2 2 t s rhot rhost t rhos rhost t rhot rhott - 1/36 -------------- - 1/4 ------------- - 1/36 ------------- , 0 , 4 4 4 k k k 2 1/2 1/2 t rhot 2 rhost t s rhos 2 rhost - 1/48 ------------------ - 1/48 ------------------- 4 4 k k 1/2 1/2 t rhot 2 rhos t s rhot 2 rhoss - 1/24 ---------------- - 1/48 ------------------- 4 4 k k 2 1/2 2 1/2 s rhot 2 rhost t rhos 2 rhoss + 1/144 ------------------ - 3/16 ------------------ 4 4 k k 3 1/2 2 2 1/2 1/2 t rhos 2 rhost t rhos 2 t s rhot 2 rhott - 3/16 ------------------ - 3/16 ------------- + 1/144 ------------------- 4 4 4 k s k s k 2 1/2 2 1/2 ] t rhos 2 rhott s rhot 2 ] - 1/48 ------------------ + 1/144 ------------ , 0 , 0 , 0] 4 4 ] k k ] [ 2 2 [ t rhot t s rhos rhoss t rhos [0 , 0 , -t + 1/36 ------- + 1/4 -------------- + 1/4 ------- [ 4 4 4 [ k k k 2 2 t s rhot rhost t rhos rhost t rhot rhott + 1/36 -------------- + 1/4 ------------- + 1/36 ------------- , 0 , 0 , 0 4 4 4 k k k 2 2 s rhos t s rhot rhoss t rhot rhost t rhot rhos , - 1/4 ------- - 1/8 -------------- - 1/8 ------------- - 1/4 ----------- 4 4 4 4 k k k k 2 2 ] s rhos rhoss t s rhos rhost t rhos rhott] - 1/4 ------------- - 3/8 -------------- + s - 1/8 -------------] 4 4 4 ] k k k ] [0 , 0 , 0 , 0 , 0 , 0 , 0] > A2e := [numer(A2[1,2]), > numer(A2[1,4]/t/sqrt(2)), > numer(A2[2,7])]; 2 2 A2e := [2 rhos s + t s rhot rhoss + t rhot rhost + 2 t rhot rhos 2 4 2 4 + 2 s rhos rhoss + 3 t s rhos rhost - 8 k s + t rhos rhott, -2 t (36 k 2 2 - rhot - 9 rhos - 9 s rhos rhoss - 9 t rhos rhost - s rhot rhost 2 2 2 2 - t rhot rhott), -s rhot + 27 t rhos rhoss s + 3 t s rhot rhoss 2 2 3 + 3 t rhot rhost s + 6 t rhot rhos s + 3 t s rhos rhost - s rhot rhost 3 2 2 2 2 + 27 t rhos rhost + 27 t rhos - t s rhot rhott + 3 t rhos rhott s] > A2c := map(collect,A2e,{rhoss,rhost,rhott},distribute,factor); 2 A2c := [t (3 s rhos + t rhot) rhost + t rhos rhott 2 4 + s (2 s rhos + t rhot) rhoss + 2 rhos s - 8 k s + 2 t rhot rhos, 2 2 t (9 t rhos + rhot s) rhost + 2 t rhot rhott + 18 t s rhos rhoss 4 2 2 - 2 t (36 k - rhot - 9 rhos ), 2 3 2 3 (3 t rhot s + 27 rhos t + 3 s t rhos - s rhot) rhost + t s (-rhot s + 3 t rhos) rhott + 3 t s (9 t rhos + rhot s) rhoss + (9 t rhos - rhot s) (3 t rhos + rhot s)] > A2e2 := [ coeff(A2c[2],rhoss)*A2c[1] - coeff(A2c[1],rhoss)*A2c[2], > coeff(A2c[3],rhoss)*A2c[1] - coeff(A2c[1],rhoss)*A2c[3] ]; 2 A2e2 := [18 t s rhos (t (3 s rhos + t rhot) rhost + t rhos rhott + s %1 rhoss 2 4 + 2 rhos s - 8 k s + 2 t rhot rhos) - s %1 ( 2 2 t (9 t rhos + rhot s) rhost + 2 t rhot rhott + 18 t s rhos rhoss 4 2 2 - 2 t (36 k - rhot - 9 rhos )), 3 t s (9 t rhos + rhot s) ( 2 2 t (3 s rhos + t rhot) rhost + t rhos rhott + s %1 rhoss + 2 rhos s 4 - 8 k s + 2 t rhot rhos) - s %1 ( 2 3 2 3 (3 t rhot s + 27 rhos t + 3 s t rhos - s rhot) rhost + t s (-rhot s + 3 t rhos) rhott + 3 t s (9 t rhos + rhot s) rhoss + (9 t rhos - rhot s) (3 t rhos + rhot s))] %1 := 2 s rhos + t rhot > A2c2 := map(collect,A2e2,{rhoss,rhost,rhott},distribute,factor); 2 2 2 A2c2 := [2 t s (-2 rhot rhos s + 9 t rhos - t rhot ) rhost 2 2 2 + 2 t s (-2 rhot rhos s + 9 t rhos - t rhot ) rhott 2 2 4 + 2 t rhot s (-2 rhot rhos s - t rhot + 9 t rhos + 36 k t), 2 3 2 2 2 3 2 2 s (27 t rhos - 6 t rhos s + 2 s rhos rhot + t s rhot ) rhost 3 2 2 2 3 2 2 + s t (27 t rhos - 6 t rhos s + 2 s rhos rhot + t s rhot ) rhott - s 2 4 3 2 2 2 2 4 (216 s t rhos k - 27 t rhot rhos + 6 s rhos t rhot + 24 s t rhot k 3 2 3 2 - 2 s rhos rhot - t rhot s )] > A2e3 := expand(coeff(A2c2[2],rhost)*A2c2[1] > - coeff(A2c2[1],rhost)*A2c2[2]); 4 4 3 4 4 4 2 4 5 3 2 4 A2e3 := 3888 s t rhos k - 432 s t rhot rhos k - 864 s t rhos rhot k 6 2 2 4 5 3 3 4 3 5 2 4 + 48 s rhos rhot t k + 24 s t rhot k + 1944 s t rhos rhot k > print(`Second equation`); Second equation > e2 := factor(A2e3); 4 2 3 2 e2 := 24 k t s (2 s rhos + t rhot) (9 t rhos - rhot s) > print(`Solutions to e2`); Solutions to e2 > sollist := [ -2*s*rhos/t, 9*t*rhos/s ]; s rhos t rhos sollist := [-2 ------, 9 ------] t s > print(`These two expressions should give zero`); These two expressions should give zero bytes used=10112140, alloc=2227816, time=15.89 > subs(rhot=sollist[1],e2); 0 > subs(rhot=sollist[2],e2); 0 > print(`Substitute first solution in e1`); Substitute first solution in e1 > factor(subs(rhot=op(1,sollist),e1)); 2 4 -s (rhos + 4 k ) > print(`Substitute second solution in e1`); Substitute second solution in e1 > s2 := simplify(subs(rhot=op(2,sollist),e1)); 2 2 2 2 4 2 -rhos s - 9 t rhos + 4 k s s2 := - -------------------------------- s > s3 := collect(numer(s2),rhos); 2 2 2 4 2 s3 := (s + 9 t ) rhos - 4 k s > print(`Solutions for rhos and rhot are`); Solutions for rhos and rhot are > solrhos := 2*k^2*s/sqrt(s^2+9*t^2); 2 k s solrhos := 2 -------------- 2 2 1/2 (s + 9 t ) > solrhot := 18*k^2*t/sqrt(s^2+9*t^2); 2 k t solrhot := 18 -------------- 2 2 1/2 (s + 9 t ) > print(`The following three expressions should be zero`); The following three expressions should be zero > subs(rhos=solrhos,s3); 0 > simplify(subs([rhos=solrhos,rhot=solrhot],e1)); 0 > subs([rhos=solrhos,rhot=solrhot],e2); 0 > int(solrhos,s)+fun1(t); 2 2 1/2 2 2 (s + 9 t ) k + fun1(t) > int(solrhot,t)+fun2(s); 2 2 1/2 2 2 (s + 9 t ) k + fun2(s) > solrho := int(solrhos,s); 2 2 1/2 2 solrho := 2 (s + 9 t ) k > print(`The following should be zero`); The following should be zero > radsimp(k*sqrt(2*(eta1+sqrt(6*(eta1^2-k^2*eta2))))) > -solrho; 0 > rhos:=diff(solrho,s); 2 k s rhos := 2 -------------- 2 2 1/2 (s + 9 t ) > rhot:=diff(solrho,t); 2 k t rhot := 18 -------------- 2 2 1/2 (s + 9 t ) > rhoss:=diff(rhos,s); 2 2 2 k k s rhoss := 2 -------------- - 2 -------------- 2 2 1/2 2 2 3/2 (s + 9 t ) (s + 9 t ) > rhost:=diff(rhos,t); 2 k s t rhost := -18 -------------- 2 2 3/2 (s + 9 t ) > rhott:=diff(rhot,t); 2 2 2 k k t rhott := 18 -------------- - 162 -------------- 2 2 1/2 2 2 3/2 (s + 9 t ) (s + 9 t ) > print(`All the remaining matrices should be zero`); All the remaining matrices should be zero > for i from 1 to 3 do > for j from 1 to 3 do > A := matrix(3,3,0); > A[i,j]:=1; > if i=j then A[3,3]:=-A[i,i] fi; > xA := so7Lb(X,sl3(A)); > print(map(simplify,J2stP(xA,1))); > od; > od; bytes used=11112456, alloc=2227816, time=17.43 bytes used=12112684, alloc=2227816, time=19.23 bytes used=13115316, alloc=2227816, time=21.65 bytes used=14125600, alloc=2293340, time=26.52 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=15125800, alloc=2293340, time=27.92 bytes used=16126132, alloc=2293340, time=29.34 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=17126340, alloc=2358864, time=30.77 bytes used=18126556, alloc=2358864, time=32.20 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=19126764, alloc=2358864, time=33.70 bytes used=20127012, alloc=2358864, time=35.44 bytes used=21127412, alloc=2358864, time=37.12 bytes used=22127876, alloc=2358864, time=38.39 bytes used=23128584, alloc=2358864, time=39.65 bytes used=24128804, alloc=2358864, time=40.88 bytes used=25129248, alloc=2358864, time=42.58 bytes used=26129648, alloc=2358864, time=43.89 bytes used=27130100, alloc=2358864, time=45.17 bytes used=28279020, alloc=2358864, time=47.03 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=29279808, alloc=2555436, time=48.66 bytes used=30280036, alloc=2555436, time=50.20 bytes used=31455428, alloc=2555436, time=53.62 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=32455716, alloc=2620960, time=55.37 bytes used=33455932, alloc=2620960, time=56.90 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=34456196, alloc=2620960, time=58.33 bytes used=35456540, alloc=2620960, time=59.89 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=36456760, alloc=2620960, time=61.56 bytes used=37457132, alloc=2620960, time=63.32 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] Warning: matrix passed to sl3 is not trace-free bytes used=38457428, alloc=2620960, time=64.85 bytes used=39457736, alloc=2620960, time=66.36 bytes used=40767288, alloc=2817532, time=69.89 [ 3 1/2 2 ] [ t 2 t s ] [0 , 0 , 0 , - --------- , 0 , - 2/3 --------- , 0] [ 2 2 2 2 ] [ s + 9 t s + 9 t ] [ ] [ 2 ] [ t s ] [0 , 0 , 0 , 0 , 0 , 0 , 2/3 ---------] [ 2 2] [ s + 9 t ] [ ] [0 , 0 , 0 , 0 , 0 , 0 , 0] [ ] [ 2 1/2 3 1/2 ] [ t s 2 t 2 ] [0 , 0 , - 1/3 --------- , 0 , 0 , 0 , ---------] [ 2 2 2 2] [ s + 9 t s + 9 t ] [ ] [ 3 2 1/2 ] [ t t s 2 ] [0 , -2 --------- , 0 , 1/3 --------- , 0 , 0 , 0] [ 2 2 2 2 ] [ s + 9 t s + 9 t ] [ ] [ 3 ] [ t ] [0 , 0 , 2 --------- , 0 , 0 , 0 , 0] [ 2 2 ] [ s + 9 t ] [ ] [0 , 0 , 0 , 0 , 0 , 0 , 0] > for i from 1 to 3 do > v:=[0,0,0]; > v[i]:=1; > xA := so7Lb(X,V10(v)); > print(map(simplify,J2stP(xA,1))); > xA := so7Lb(X,V01(v)); > print(map(simplify,J2stP(xA,1))); > od; v := [0, 0, 0] v[1] := 1 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] xA := [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] bytes used=41767624, alloc=2817532, time=71.76 bytes used=42767836, alloc=2817532, time=73.28 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ 1/2 1/2 ] [t 2 , 0 , 0 , 0 , - 1/2 s 2 , 0 , 0] [ ] [ 1/2 ] [0 , - 1/2 t 2 , 0 , 0 , 0 , 0 , 0] [ ] [ 1/2 1/2] [0 , 0 , - 1/2 t 2 , 0 , 0 , 0 , 1/2 s 2 ] [ ] xA := [0 , -s , 0 , 0 , 0 , 0 , 0] [ ] [ 1/2 ] [0 , 0 , 0 , 0 , 1/2 t 2 , 0 , 0] [ ] [ 1/2 ] [0 , 0 , 0 , s , 0 , 1/2 t 2 , 0] [ ] [ 1/2] [0 , 0 , 0 , 0 , 0 , 0 , -t 2 ] bytes used=43768040, alloc=2817532, time=74.76 bytes used=44768288, alloc=2817532, time=76.41 bytes used=45768500, alloc=2817532, time=78.37 bytes used=46768864, alloc=2817532, time=79.58 bytes used=47769256, alloc=2817532, time=80.83 bytes used=48769516, alloc=2817532, time=82.05 bytes used=49769924, alloc=2817532, time=83.30 bytes used=50770316, alloc=2817532, time=84.65 bytes used=51770660, alloc=2817532, time=85.87 bytes used=52771260, alloc=2817532, time=87.12 bytes used=53771468, alloc=2817532, time=88.33 bytes used=54771720, alloc=2817532, time=89.58 bytes used=55771988, alloc=2817532, time=90.84 bytes used=56772272, alloc=2817532, time=93.18 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] v := [0, 0, 0] v[2] := 1 [ 1/2 ] [0 0 0 s 0 -t 2 0 ] [ ] [ 1/2] [0 0 0 0 0 0 t 2 ] [ ] [0 0 0 0 0 0 0 ] [ ] xA := [0 0 2 t 0 0 0 -s ] [ ] [ 1/2 ] [0 - 1/2 s 2 0 -2 t 0 0 0 ] [ ] [ 1/2 ] [0 0 1/2 s 2 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] bytes used=57772644, alloc=2817532, time=95.06 bytes used=58773524, alloc=2817532, time=96.85 bytes used=60089240, alloc=2817532, time=101.53 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ 1/2 ] [0 3/2 t 2 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] xA := [0 0 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] [ 1/2] [0 0 0 0 0 0 - 3/2 t 2 ] [ ] [0 0 0 0 0 0 0 ] bytes used=61089588, alloc=2817532, time=103.42 bytes used=62089912, alloc=2817532, time=104.80 bytes used=63429108, alloc=2817532, time=108.52 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] v := [0, 0, 0] v[3] := 1 [ 1/2 ] [0 0 0 0 -t 2 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] [ 1/2] [0 0 0 0 0 0 t 2 ] xA := [ ] [0 -2 t 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] [0 0 0 2 t 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] bytes used=64429320, alloc=2817532, time=110.34 bytes used=65429564, alloc=2817532, time=111.88 bytes used=66430064, alloc=2817532, time=113.44 bytes used=67430428, alloc=2817532, time=114.67 bytes used=68430832, alloc=2817532, time=115.98 bytes used=69431136, alloc=2817532, time=117.24 bytes used=70460092, alloc=2817532, time=118.77 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ 1/2 ] [0 0 3/2 t 2 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] [ ] xA := [0 0 0 0 0 0 0 ] [ ] [ 1/2] [0 0 0 0 0 0 - 3/2 t 2 ] [ ] [0 0 0 0 0 0 0 ] [ ] [0 0 0 0 0 0 0 ] bytes used=71460332, alloc=2817532, time=120.37 bytes used=72460560, alloc=2817532, time=121.95 [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] [ ] [0 0 0 0 0 0 0] > quit; bytes used=72775452, alloc=2817532, time=122.39