Reading file: Data/LS.2D.16.r.dat n = 16 m + m_ = 104 (n - 1)*d - d = 28 n_active = 16 m_active = 40 The average, max and min relative gaps are: 1.2241055595586683 0.9577418144371066 0 The number of perfect (to within 1E-3) contacts is 6 Now at compression attempt 1 LOQO 6.03: iterlim=100 convex verbose=1 It's a QP. 1 0.000000e+00 9.9e+02 4.278411e+04 2.0e+03 2 -4.427967e+01 4.9e+02 3.666305e+04 3.8e-13 DF 3 -2.721532e+00 3.0e+01 2.745226e+04 1.1e-12 DF 4 -3.181665e-01 3.5e+00 4.025890e+03 9.0e-14 DF 5 -3.075589e-02 3.4e-01 2.824948e+02 1.0e-12 DF 6 -3.355424e-05 2.2e-03 3.252854e+01 8.0e-13 DF 7 1.833565e-04 3.1e-13 1.683780e+00 1.6e-13 PF DF 8 2.038897e-04 4.1e-13 1.668114e-01 1.2e-14 1 PF DF 9 8.680816e-04 3.9e-13 8.601906e-02 6.1e-15 1 PF DF 10 9.935039e-03 3.2e-13 7.093008e-02 1.5e-14 1 PF DF 11 1.149042e-02 3.2e-13 1.837092e-02 3.7e-15 2 PF DF 12 1.252813e-02 1.9e-12 1.362484e-02 3.9e-16 3 PF DF 13 1.269866e-02 4.5e-12 1.294672e-02 2.4e-16 4 PF DF 14 1.273834e-02 3.9e-11 1.279816e-02 5.7e-16 4 PF DF 15 1.274542e-02 5.0e-11 1.276344e-02 5.7e-15 5 PF DF 16 1.274747e-02 1.3e-11 1.275040e-02 2.4e-14 6 PF DF 17 1.274797e-02 7.5e-11 1.274813e-02 4.4e-15 7 PF DF 18 1.274800e-02 7.9e-11 1.274801e-02 1.1e-14 8 PF DF D, dD, dD/D for id= 1 = 0.339924028545668 0.012748000352406502 0.03750249844633692 Strain is: -1.7398416717271639 0.00021302361728656104 0.00021302361728656104 1.739841671566896 Average relative nodal displacement: 6.238378140280849 Scaling for unjamming motion is: 0.042684011293378024 Updating frame # 1 Average: 0.2658534362252299 and maximal: 0.4404697216194792 relative total displacement Now at compression attempt 2 LOQO 6.03: iterlim=100 convex verbose=1 It's a QP. 1 1.274800e-02 1.9e+02 4.285767e+04 2.5e+03 2 -4.491318e+01 9.5e+01 3.665959e+04 6.8e-13 DF 3 -3.277841e+00 6.9e+00 2.753491e+04 1.4e-12 DF 4 -1.025907e+00 2.2e+00 4.677574e+03 2.0e-12 DF 5 -3.201546e-01 6.6e-01 7.681864e+02 1.1e-12 DF 6 -6.504544e-02 1.4e-01 1.456783e+02 5.5e-13 DF 7 -5.397505e-03 1.3e-02 1.504523e+01 2.6e-13 DF 8 5.932419e-04 1.6e-15 1.857900e+00 5.5e-13 PF DF 9 6.047400e-04 1.5e-15 1.103215e-01 1.2e-14 1 PF DF 10 8.565113e-04 1.4e-15 2.287896e-02 1.9e-15 2 PF DF 11 1.139768e-02 1.4e-15 1.584255e-02 1.2e-15 2 PF DF 12 1.304129e-02 1.5e-15 1.330231e-02 8.6e-17 4 PF DF 13 1.312947e-02 1.3e-15 1.316750e-02 1.5e-16 4 PF DF 14 1.313898e-02 1.4e-15 1.314954e-02 1.3e-16 5 PF DF 15 1.314336e-02 1.3e-15 1.314606e-02 6.2e-16 6 PF DF 16 1.314435e-02 1.2e-15 1.314494e-02 1.1e-15 6 PF DF 17 1.314455e-02 1.5e-15 1.314469e-02 4.5e-15 7 PF DF 18 1.314459e-02 1.1e-15 1.314462e-02 2.5e-14 8 PF DF D, dD, dD/D for id= 1 = 0.3404681643366781 0.013144592080478494 0.03860740432541654 Strain is: 0.0025043284592766815 -2.805540081197848e-05 -2.805540081197848e-05 -0.002504333023558513 Average relative nodal displacement: 0.0005539728917005788 Scaling for unjamming motion is: 1.0000322278102465 Updating frame # 2 Average: 0.24493111209492074 and maximal: 0.40243334004066145 relative total displacement Now at compression attempt 3 LOQO 6.03: iterlim=100 convex verbose=1 It's a QP. 1 1.314459e-02 1.0e+03 4.277702e+04 2.5e+03 2 -4.491539e+01 5.1e+02 3.665862e+04 1.2e-13 DF 3 -2.628650e+00 3.0e+01 2.753715e+04 6.8e-13 DF 4 -2.459191e-01 2.8e+00 3.871443e+03 1.2e-12 DF 5 -3.101504e-02 3.5e-01 2.381773e+02 9.4e-13 DF 6 -1.766359e-03 2.0e-02 1.222626e+01 4.6e-13 DF 7 -6.423297e-05 7.6e-04 6.313791e-01 1.2e-12 DF 8 2.530967e-06 6.0e-15 4.270919e-02 8.0e-13 1 PF DF 9 2.553003e-06 8.5e-15 2.313337e-03 3.0e-14 3 PF DF 10 3.012518e-06 7.2e-15 2.924226e-04 6.3e-15 4 PF DF 11 1.853381e-05 9.3e-15 1.716854e-04 3.1e-15 4 PF DF 12 5.160917e-05 6.6e-15 6.622302e-05 4.2e-16 5 PF DF 13 5.508660e-05 7.1e-15 5.584737e-05 2.5e-16 6 PF DF 14 5.526239e-05 7.0e-15 5.531817e-05 1.0e-16 7 PF DF 15 5.527758e-05 5.2e-15 5.528746e-05 2.7e-16 8 PF DF D, dD, dD/D for id= 1 = 0.35361318003857595 5.5277580118746375e-05 0.0001563221713419056 Strain is: -2.453558813523727e-05 1.9319750447926897e-05 1.9319750447926897e-05 2.4531996510516868e-05 Average relative nodal displacement: 0.00019982845954991034 Scaling for unjamming motion is: 1.0000399647973486 Updating frame # 3 Average: 0.24481640345956146 and maximal: 0.4021102525985475 relative total displacement Now at compression attempt 4 LOQO 6.03: iterlim=100 convex verbose=1 It's a QP. 1 5.527758e-05 1.0e+03 4.277659e+04 2.5e+03 2 -4.492194e+01 5.2e+02 3.665845e+04 1.1e-12 DF 3 -2.629202e+00 3.0e+01 2.753704e+04 6.5e-13 DF 4 -2.461607e-01 2.9e+00 3.871449e+03 7.7e-13 DF 5 -3.091083e-02 3.6e-01 2.382414e+02 9.4e-13 DF 6 -1.786160e-03 2.1e-02 1.220488e+01 1.4e-13 DF 7 -8.933934e-05 1.0e-03 6.102644e-01 2.8e-13 DF 8 -4.458679e-06 5.2e-05 3.052030e-02 7.1e-13 2 DF 9 -2.146408e-07 2.5e-06 1.533071e-03 8.7e-14 3 DF 10 -2.438844e-09 3.8e-08 8.322292e-05 6.3e-13 4 PF DF 11 8.719720e-10 7.9e-15 4.740122e-06 6.6e-14 5 PF DF 12 9.071177e-10 1.0e-14 3.258879e-07 8.2e-15 6 PF DF 13 1.742098e-09 7.4e-15 1.021899e-07 2.8e-15 7 PF DF 14 1.558175e-08 5.3e-15 4.531589e-08 2.4e-16 8 PF DF D, dD, dD/D for id= 1 = 0.353668459827852 1.5581748921818448e-08 4.4057502128979384e-08 Strain is: -2.5009107173546877e-09 -5.767570529695549e-09 -5.767570529695549e-09 1.7904082588577555e-09 Average relative nodal displacement: 3.807337303939402e-08 Scaling for unjamming motion is: 1.0538780070967726 Updating frame # 4 Average: 0.24481641558189898 and maximal: 0.402110306088592 relative total displacement Now at compression attempt 5 LOQO 6.03: iterlim=100 convex verbose=1 It's a QP. 1 1.558175e-08 1.0e+03 4.277659e+04 2.5e+03 2 -4.492197e+01 5.2e+02 3.665845e+04 1.4e-12 DF 3 -2.629205e+00 3.0e+01 2.753703e+04 1.3e-12 DF 4 -2.461617e-01 2.9e+00 3.871450e+03 8.8e-13 DF 5 -3.090953e-02 3.6e-01 2.382420e+02 3.7e-13 DF 6 -1.786098e-03 2.1e-02 1.220488e+01 5.4e-13 DF 7 -8.934346e-05 1.0e-03 6.102579e-01 7.7e-13 DF 8 -4.466110e-06 5.2e-05 3.051381e-02 5.4e-13 2 DF 9 -2.222376e-07 2.6e-06 1.526602e-03 2.8e-12 3 DF 10 -1.004400e-08 1.2e-07 7.723341e-05 4.6e-13 4 PF DF 11 1.119006e-10 9.3e-15 4.608993e-06 5.4e-13 5 PF DF 12 8.407778e-11 3.5e-09 2.934340e-07 6.8e-14 7 PF DF 13 1.185423e-10 7.7e-15 2.614194e-08 4.0e-15 8 PF DF D, dD, dD/D for id= 1 = 0.3536684762491145 1.1854225775234314e-10 3.351790326623433e-10 Strain is: -1.4644472091777107e-09 1.2832957166768318e-10 1.2832957166768318e-10 -4.5755508526498663e-10 Average relative nodal displacement: 2.9511066154096034e-09 Scaling for unjamming motion is: 1.5334159720330134 Updating frame # 5 Average: 0.24481641379986024 and maximal: 0.4021102974013991 relative total displacement n = 16 m + m_ = 104 (n - 1)*d - d = 28 l_min = 0.353668 l_max = 0.707337 The average, max and min relative gaps are: 1.1320508257092963 1.000000013525978 -6.278326164826323e-16 The number of perfect (to within 1E-3) contacts is 48 Writing compressed packing to: Packings/LS.2D.16.D.dat Writing VRML animation to: Scenes/LS.2D.16.compression.wrl with 5 frames!