24-Dec-2018 09:01:51 CVT_1D_LUMPING MATLAB version Use lumping to implement a version of Lloyd's algorithm in the 1D interval [-1,1]. Number of generators is 33 Number of iterations is 1000 Number of samples is 4000000 Initial generators (Uniform, staggered): 1 -0.941176 2 -0.882353 3 -0.823529 4 -0.764706 5 -0.705882 6 -0.647059 7 -0.588235 8 -0.529412 9 -0.470588 10 -0.411765 11 -0.352941 12 -0.294118 13 -0.235294 14 -0.176471 15 -0.117647 16 -0.058824 17 0.000000 18 0.058824 19 0.117647 20 0.176471 21 0.235294 22 0.294118 23 0.352941 24 0.411765 25 0.470588 26 0.529412 27 0.588235 28 0.647059 29 0.705882 30 0.764706 31 0.823529 32 0.882353 33 0.941176 Saved graphics file "density.png" Saved graphics file "energy_033_u.png" Saved graphics file "motion_033_u.png" Saved graphics file "evolution_033_u.png" CVT data: K G H rho H^3 mu H 1 -1.0000 0.004597 0.1541 0.5362 2 -0.9908 0.0136 0.001023 0.1007 3 -0.9728 0.02249 0.0009125 0.09699 4 -0.9459 0.03116 0.0008855 0.09603 5 -0.9104 0.03956 0.0008747 0.09564 6 -0.8668 0.0476 0.0008694 0.09544 7 -0.8152 0.05521 0.0008663 0.09533 8 -0.7563 0.06232 0.0008644 0.09526 9 -0.6906 0.06886 0.0008632 0.09521 10 -0.6186 0.07479 0.0008624 0.09518 11 -0.5410 0.08003 0.0008618 0.09516 12 -0.4585 0.08455 0.0008613 0.09515 13 -0.3719 0.08831 0.000861 0.09514 14 -0.2819 0.09127 0.0008608 0.09513 15 -0.1894 0.0934 0.0008607 0.09512 16 -0.0951 0.09469 0.0008606 0.09512 17 0.0000 0.09509 0.0008599 0.09509 18 0.0951 0.09464 0.0008592 0.09507 19 0.1893 0.09335 0.0008594 0.09507 20 0.2818 0.09123 0.0008596 0.09508 21 0.3717 0.08828 0.0008598 0.09509 22 0.4583 0.08453 0.0008602 0.0951 23 0.5408 0.08001 0.0008607 0.09512 24 0.6184 0.07478 0.0008613 0.09515 25 0.6903 0.06886 0.0008623 0.09518 26 0.7561 0.06233 0.0008636 0.09523 27 0.8150 0.05522 0.0008655 0.0953 28 0.8665 0.04762 0.0008687 0.09542 29 0.9102 0.03959 0.0008741 0.09561 30 0.9457 0.0312 0.0008849 0.09601 31 0.9726 0.02253 0.000912 0.09698 32 0.9908 0.01365 0.001021 0.1007 33 0.9999 0.004644 0.07649 0.4245 Chebyshev Zero Node data for comparison K G H rho H^3 mu H 1 -0.9989 0.005656 0.001679 0.1189 2 -0.9898 0.01353 0.0008589 0.09506 3 -0.9718 0.02241 0.0008589 0.09506 4 -0.9450 0.03109 0.0008589 0.09506 5 -0.9096 0.03949 0.0008589 0.09506 6 -0.8660 0.04753 0.0008589 0.09506 7 -0.8146 0.05514 0.0008589 0.09506 8 -0.7557 0.06225 0.0008589 0.09506 9 -0.6901 0.0688 0.0008589 0.09506 10 -0.6182 0.07472 0.0008589 0.09506 11 -0.5406 0.07997 0.0008589 0.09506 12 -0.4582 0.08449 0.0008589 0.09506 13 -0.3717 0.08825 0.0008589 0.09506 14 -0.2817 0.09121 0.0008589 0.09506 15 -0.1893 0.09334 0.0008589 0.09506 16 -0.0951 0.09463 0.0008589 0.09506 17 0.0000 0.09506 0.0008589 0.09506 18 0.0951 0.09463 0.0008589 0.09506 19 0.1893 0.09334 0.0008589 0.09506 20 0.2817 0.09121 0.0008589 0.09506 21 0.3717 0.08825 0.0008589 0.09506 22 0.4582 0.08449 0.0008589 0.09506 23 0.5406 0.07997 0.0008589 0.09506 24 0.6182 0.07472 0.0008589 0.09506 25 0.6901 0.0688 0.0008589 0.09506 26 0.7557 0.06225 0.0008589 0.09506 27 0.8146 0.05514 0.0008589 0.09506 28 0.8660 0.04753 0.0008589 0.09506 29 0.9096 0.03949 0.0008589 0.09506 30 0.9450 0.03109 0.0008589 0.09506 31 0.9718 0.02241 0.0008589 0.09506 32 0.9898 0.01353 0.0008589 0.09506 33 0.9989 0.005656 0.001679 0.1189 ||G-ChebyshevZero|| = 0.00330913 Saved CVT coordinates in file "cvt_033_u.txt" Saved Chebyshev Zero node coordinates in file "chebyzero_033.txt" Iteration required 26.3245 seconds. CVT_1D_LUMPING Normal end of execution. 24-Dec-2018 09:02:20 Output saved in "output_033_u.txt"