diff --git a/Benchmarks/ConfigFiles/m1931_modif_210Img.txt b/Benchmarks/ConfigFiles/m1931_modif_210Img.txt
index 535aeac..6a43880 100644
--- a/Benchmarks/ConfigFiles/m1931_modif_210Img.txt
+++ b/Benchmarks/ConfigFiles/m1931_modif_210Img.txt
@@ -1,208 +1,228 @@
#REFERENCE 0 0.0000000 0.0000000
-1 2.6786527 2.9037603 34.846982 0.044513 62.43534 1.5000 7.31
-1 -11.8484881 13.3774383 41.530182 0.058299 177.79656 1.5000 6.83
-1 9.9036907 -10.4750539 34.634955 0.075692 176.98494 1.5000 6.74
-2 -2.6672309 18.3345021 41.997432 0.110967 173.97052 1.5000 6.12
-2 10.9447338 9.1407289 21.765016 0.201954 53.10618 1.5000 6.18
-2 15.5033588 1.7617761 21.466066 0.359257 2.07257 1.5000 5.57
-3 2.5286856 2.1814895 8.491896 0.162661 47.99557 1.5000 7.44
-3 -12.8724137 13.6578294 6.819269 0.276049 22.52509 1.5000 7.10
-3 9.0081475 -10.0961433 6.817751 0.353202 31.25535 1.5000 6.84
-4 -1.3358012 -2.8790666 23.057318 0.057944 27.85238 1.5000 7.48
-4 6.8407717 -8.1121706 20.574407 0.084857 91.24205 1.5000 7.19
-4 -14.8580323 14.8551029 11.821152 0.067028 86.10019 1.5000 8.04
-5 3.0331801 -14.6505707 44.887170 0.021107 16.70145 1.5000 7.85
-5 -16.7625653 10.1664425 28.244176 0.021254 39.18539 1.5000 8.35
-5 -8.4052010 -5.9872161 23.055401 0.027855 44.51135 1.5000 8.27
-6 -11.5124223 15.9970478 40.696439 0.013649 36.64909 1.5000 8.43
-6 10.1123445 -6.2419290 56.155712 0.018888 63.79062 1.5000 7.73
-6 6.8681360 4.5905742 38.939242 0.020371 11.59938 1.5000 8.05
-7 24.4507222 35.4518592 12.030922 0.300848 1.46343 1.5000 6.39
-8 0.4294637 0.3368837 17.875325 0.032617 41.79915 1.5000 8.38
-8 1.6975484 1.4120697 17.112294 0.036880 38.61008 1.5000 8.30
-8 -2.6231902 -2.2017254 12.024044 0.043511 37.79877 1.5000 8.50
-8 8.4786219 -10.5709866 7.510960 0.077451 57.12046 1.5000 8.38
-8 -13.4278497 13.3064510 6.248381 0.071115 48.69072 1.5000 8.68
-9 -12.2729271 23.8607406 20.470970 0.136267 2.27813 1.5000 6.68
-10 8.3811315 7.9105283 16.010105 0.133434 131.96659 1.5000 6.96
-10 13.6029638 -3.8834001 17.642846 0.159210 104.80730 1.5000 6.66
-10 -6.8153537 16.9665578 9.442873 0.187264 133.56479 1.5000 7.17
-11 -5.1498385 15.3208303 173.626665 0.031134 14.17090 1.5000 5.95
-11 14.4490163 -7.0793369 140.119512 0.032200 49.70954 1.5000 6.15
-11 5.5394577 10.3171621 119.187491 0.048467 84.58532 1.5000 5.88
-13 7.1884598 14.6438521 8.810685 0.197648 149.64275 1.5000 7.19
-13 3.4752244 16.8592447 5.944581 0.243989 146.03268 1.5000 7.39
-13 17.4301938 2.1957167 2.643589 0.166403 100.27012 1.5000 8.68
-15 15.5982332 -6.8318229 25.200524 0.009008 63.02755 1.5000 9.40
-15 -2.4642730 15.0794757 31.326845 0.012990 15.38250 1.5000 8.77
-15 3.9724585 12.3620713 21.562202 0.021541 103.31489 1.5000 8.63
-18 1.0479069 -1.1044685 10.768333 0.006847 22.52448 1.5000 10.62
-18 -15.2217938 15.3201427 5.277769 0.004422 164.44813 1.5000 11.87
-18 6.2890313 -7.2140096 8.200618 0.007841 155.73619 1.5000 10.77
-18 -0.7098203 -2.9315479 6.970110 0.009213 81.11879 1.5000 10.77
-19 3.6379121 -17.6226603 7.147027 0.110295 19.29235 1.5000 8.04
-19 -16.0320537 5.4916741 5.593126 0.221276 35.01584 1.5000 7.55
-19 -11.8452915 -3.3468708 4.615314 0.203069 43.22132 1.5000 7.85
-21 -20.5880303 8.9329582 0.336281 0.012562 64.84756 1.5000 9.16
-22 17.3720134 2.1882537 0.401604 0.028619 88.59146 1.5000 7.71
-22 7.5287582 14.3644280 0.621351 0.062721 157.01493 1.5000 6.38
-22 3.0852536 17.0240154 0.263595 0.118443 129.67801 1.5000 6.62
-23 -15.6866276 15.6845567 0.055011 0.007139 161.67654 1.5000 13.42
-23 5.5055808 -6.4106321 0.103302 0.014270 154.36511 1.5000 11.99
-23 -0.1804981 -3.2185423 0.079898 0.021746 92.26818 1.5000 11.81
-24 16.5899706 -0.7754639 0.712518 0.072540 87.02951 1.5000 8.20
-24 7.3400718 12.9906877 0.760430 0.122910 154.27397 1.5000 7.56
-24 0.0386740 17.0817512 0.360582 0.204579 144.39873 1.5000 7.82
-25 -0.5582720 18.1081624 0.501990 0.040558 170.87681 1.5000 12.51
-25 10.3472270 11.1997347 0.595960 0.072387 123.22072 1.5000 11.69
-25 16.4070668 2.0675156 0.195808 0.113744 29.24929 1.5000 12.41
-26 -3.6566538 -0.9732215 0.613580 0.136538 20.52408 1.5000 8.39
-26 9.7897135 -12.7521977 0.387439 0.165176 60.03832 1.5000 8.68
-26 -11.8832764 11.3910034 0.388550 0.197227 54.76396 1.5000 8.49
-28 -8.8301584 -10.2296004 6.278205 0.234335 139.37575 1.5000 6.45
-28 -4.0456315 -14.0621632 4.996734 0.265758 137.81187 1.5000 6.56
-28 -20.3074952 8.7186627 1.783389 0.229330 104.42894 1.5000 7.84
-29 -15.4001065 14.6578526 0.554882 0.080878 168.20715 1.5000 7.98
-29 6.1406896 -8.3811360 0.850465 0.127398 160.62235 1.5000 7.02
-29 -1.9455502 -3.6483147 0.607322 0.104519 75.02974 1.5000 7.60
-30 16.9039033 -0.7407439 0.359729 0.006926 79.75539 1.5000 10.59
-30 1.7452545 16.5076313 0.192199 0.032516 12.43573 1.5000 9.59
-30 5.9668771 14.1843117 0.142788 0.047680 18.69737 1.5000 9.50
-31 16.5139312 -3.8451362 0.953037 0.000768 78.65420 1.5000 13.92
-31 1.3123810 15.4732942 0.736126 0.005786 8.96831 1.5000 12.01
-31 3.6765900 14.3002919 0.396902 0.009968 50.91878 1.5000 12.09
-32 -19.3763552 6.4541388 0.128487 0.031001 54.30277 1.5000 10.09
-32 -1.0478628 -16.8624882 0.128406 0.034213 18.77814 1.5000 9.98
-32 -11.6234609 -7.1394199 0.093792 0.032315 27.56403 1.5000 10.38
-33 4.4835691 5.8872015 0.940979 0.014051 53.91017 1.5000 12.37
-33 11.5215383 -9.5055949 1.026644 0.019262 61.39061 1.5000 11.94
-33 -9.7699119 14.0765786 0.908140 0.022506 44.63471 1.5000 11.90
-34 -10.3351451 10.6553738 0.528612 0.034140 157.42627 1.5000 13.51
-34 10.9153487 -13.3958584 0.355411 0.030735 151.89381 1.5000 14.05
-34 -0.1692590 2.6923973 0.460052 0.085688 92.53567 1.5000 12.66
-35 4.9563344 -16.6556871 1.930383 0.327434 22.06552 1.5000 7.96
-35 -15.3680550 7.0269874 1.355196 0.533507 46.15432 1.5000 7.81
-35 -10.6479548 -3.3933905 1.012174 0.499609 26.32416 1.5000 8.20
-37 -2.1269883 17.7618777 0.839149 0.146064 4.33839 1.5000 10.16
-37 15.9185406 -0.2597843 0.821566 0.184455 67.47268 1.5000 9.93
-37 9.4646689 10.7235580 0.648071 0.368732 111.52255 1.5000 9.44
-38 -6.2210215 -15.8704320 0.937844 0.024682 150.95757 1.5000 12.26
-38 -13.2038524 -10.2148699 0.700936 0.031004 127.61932 1.5000 12.33
-38 -22.1953724 5.4720839 0.242843 0.040746 173.16403 1.5000 13.18
-39 3.0312580 -15.8087360 1.252698 0.234004 3.25013 1.5000 8.95
-39 -9.7784925 -5.3082813 0.738278 0.286702 83.26595 1.5000 9.30
-39 -16.7500612 8.6158344 0.771764 0.333547 8.36658 1.5000 9.09
-40 14.4943275 -4.2786351 0.899064 0.075604 70.33597 1.5000 7.07
-40 -5.2517624 16.7019063 0.654067 0.082554 26.17356 1.5000 7.32
-40 7.7887901 9.4448365 0.461530 0.148344 29.43926 1.5000 7.06
-41 -21.3995271 7.4375134 0.805568 0.051520 45.86022 1.5000 10.84
-41 -6.4802887 -13.9904447 1.839369 0.135518 163.28206 1.5000 8.89
-41 -9.5515612 -11.4941910 1.284929 0.250311 107.80091 1.5000 8.61
-42 -13.1218821 -2.4426522 2.740368 0.076499 108.90595 1.5000 7.35
-42 -15.9520465 3.4705114 2.675136 0.095900 126.65988 1.5000 7.13
-42 3.2709225 -18.5152424 1.415411 0.073096 157.24543 1.5000 8.12
-43 2.4040917 -15.7224455 1.788478 0.078973 3.45380 1.5000 11.06
-43 -17.1344846 8.7164655 1.119390 0.106268 14.08818 1.5000 11.25
-43 -9.7154364 -5.7390332 1.035433 0.100593 75.58388 1.5000 11.39
-44 6.0647302 12.9685543 3.176132 0.076161 149.54889 1.5000 9.26
-44 16.2815642 -2.9996360 1.610419 0.078817 101.50112 1.5000 9.96
-44 -0.7772337 16.5018612 2.063908 0.102766 149.23518 1.5000 9.40
-45 9.4626405 -6.2544580 1.954653 0.014537 64.97761 1.5000 11.32
-45 -12.2639022 15.9438784 1.317880 0.010645 39.84226 1.5000 12.09
-45 6.5910853 3.9203102 1.399579 0.015605 4.33053 1.5000 11.61
-46 -3.0687036 -17.7236144 0.356505 0.042735 171.67775 1.5000 9.74
-46 -13.9649508 -8.2788025 0.235761 0.056968 21.99631 1.5000 9.88
-46 -20.6025554 4.0245322 0.274445 0.072585 38.35050 1.5000 9.45
-47 7.8688411 -16.4944952 2.140231 0.052530 36.56283 1.5000 5.67
-47 -13.2317046 6.5044642 2.135389 0.089221 58.89322 1.5000 5.10
-47 -9.9045632 -1.3440293 2.118868 0.101374 171.47875 1.5000 4.97
-48 -11.8767244 -5.4322162 1.277910 0.000346 122.06425 1.5000 6.19
-48 -17.8993462 5.9518549 1.256308 0.000450 122.67098 1.5000 5.92
-48 1.1979541 -17.4473601 1.104821 0.000408 143.89249 1.5000 6.17
-49 15.1040611 2.3571116 20.405094 0.181247 112.06935 1.5000 5.63
-49 -3.2500200 18.5620044 9.660994 0.185627 143.94883 1.5000 6.41
-49 11.0148057 8.7776669 7.289802 0.142045 128.95353 1.5000 7.01
-50 4.0302890 1.6646021 1.472455 0.043042 52.02903 1.5000 11.25
-50 -14.3131270 15.3193885 0.715348 0.038570 19.59704 1.5000 12.15
-50 7.4497300 -7.3163268 0.854058 0.072140 27.68554 1.5000 11.28
-51 -3.3306534 15.3708033 1.049641 0.198236 3.90524 1.5000 8.52
-51 4.9104924 11.6912630 0.876250 0.255629 120.78366 1.5000 8.44
-51 15.2979744 -6.5626779 0.655289 0.224586 52.14670 1.5000 8.90
-53 12.4238247 -9.5226491 2.782814 0.158592 59.54957 1.5000 5.89
-53 4.4819744 7.4432173 4.128941 0.264958 70.80977 1.5000 4.90
-53 -8.4071970 13.9983475 2.582500 0.205369 39.71261 1.5000 5.69
-55 13.7061891 1.7502961 5.328571 0.056648 109.47795 1.5000 6.36
-55 10.8979905 7.0403978 3.840162 0.048798 128.49451 1.5000 6.88
-55 -5.6997319 18.4884138 1.410669 0.058820 140.27996 1.5000 7.76
-56 6.0935140 -16.6090186 1.581300 0.092701 33.50078 1.5000 7.39
-56 -14.6184750 6.8227033 1.472476 0.122279 63.75621 1.5000 7.17
-56 -10.5002861 -2.6646404 1.105215 0.117458 175.23476 1.5000 7.52
-57 16.5121525 -3.5196819 0.091809 0.001033 78.06127 1.5000 12.52
-57 0.8381572 15.7962753 0.075629 0.004756 7.95312 1.5000 11.07
-57 4.3548701 14.0324012 0.037846 0.008976 53.77684 1.5000 11.13
-58 6.1939611 7.6041748 3.148253 0.781497 148.17302 1.5000 5.26
-58 12.4841128 -7.6487509 3.288863 0.829458 100.17621 1.5000 5.15
-58 -8.4859528 15.2787147 2.403395 1.040799 125.67280 1.5000 5.24
-59 1.2664786 2.7720698 0.338377 0.058125 35.35947 1.5000 13.41
-59 10.2287300 -12.1803663 0.215495 0.074500 76.11091 1.5000 13.63
-59 -11.3726996 11.9228981 0.200816 0.095615 73.42239 1.5000 13.44
-60 14.2660661 -10.8918824 1.461158 0.020111 48.44266 1.5000 6.97
-60 -4.2386057 12.1078222 1.845690 0.043224 19.30371 1.5000 5.89
-62 1.6682920 3.3763413 0.732986 0.064828 53.95586 1.5000 8.12
-62 -10.8262629 12.0695095 0.581357 0.084547 19.10402 1.5000 8.08
-62 10.6503513 -12.0025208 0.473102 0.082742 29.42353 1.5000 8.33
-63 1.1545882 -1.0972397 0.741401 0.013959 30.83257 1.5000 11.72
-63 -0.6725456 -3.0000851 0.466250 0.018268 64.73455 1.5000 11.93
-63 -15.2962124 15.3637616 0.236516 0.012636 0.67847 1.5000 13.07
-63 6.1748702 -7.1198104 0.347677 0.024568 175.97714 1.5000 11.93
-64 -2.6973417 0.2940260 0.482623 0.019284 42.39270 1.5000 9.32
-64 -10.6855249 10.9131743 0.241891 0.018527 161.38476 1.5000 10.12
-64 10.6884502 -13.1698494 0.170452 0.017441 157.84230 1.5000 10.56
-65 2.3094589 17.6823759 1.660124 0.394168 164.03633 1.5000 11.59
-65 9.2854380 13.3080509 2.051221 0.492479 140.77914 1.5000 11.12
-65 17.2189593 3.0433873 1.046768 0.369445 86.13514 1.5000 12.16
-66 16.7834979 -1.7923174 2.029341 0.021093 80.57855 1.5000 11.19
-66 1.6352172 16.1829209 1.262533 0.108044 10.85484 1.5000 9.93
-66 5.2103675 14.2777820 0.777049 0.179898 35.13864 1.5000 9.90
-68 -14.5875156 10.8515971 0.586897 0.165819 124.75138 1.5000 6.65
-68 -7.6472594 -3.8905882 0.650377 0.190278 138.08826 1.5000 6.39
-68 6.9882215 -13.5565861 0.532227 0.226926 129.66115 1.5000 6.42
-69 2.6039489 3.6745080 2.930790 0.046937 31.86086 1.5000 11.45
-69 10.7289822 -11.1572676 1.939554 0.057465 58.34173 1.5000 11.68
-69 -10.7738999 12.8109455 1.862324 0.067971 46.64140 1.5000 11.54
-70 6.0935140 -16.6090186 1.581300 0.092701 33.50078 1.5000 7.39
-70 -14.6184750 6.8227033 1.472476 0.122279 63.75621 1.5000 7.17
-70 -10.5002861 -2.6646404 1.105215 0.117458 175.23476 1.5000 7.52
-71 16.5121525 -3.5196819 0.091809 0.001033 78.06127 1.5000 12.52
-71 0.8381572 15.7962753 0.075629 0.004756 7.95312 1.5000 11.07
-71 4.3548701 14.0324012 0.037846 0.008976 53.77684 1.5000 11.13
-72 6.1939611 7.6041748 3.148253 0.781497 148.17302 1.5000 5.26
-72 12.4841128 -7.6487509 3.288863 0.829458 100.17621 1.5000 5.15
-72 -8.4859528 15.2787147 2.403395 1.040799 125.67280 1.5000 5.24
-73 1.2664786 2.7720698 0.338377 0.058125 35.35947 1.5000 13.41
-73 10.2287300 -12.1803663 0.215495 0.074500 76.11091 1.5000 13.63
-73 -11.3726996 11.9228981 0.200816 0.095615 73.42239 1.5000 13.44
-74 14.2660661 -10.8918824 1.461158 0.020111 48.44266 1.5000 6.97
-74 -4.2386057 12.1078222 1.845690 0.043224 19.30371 1.5000 5.89
-75 1.6682920 3.3763413 0.732986 0.064828 53.95586 1.5000 8.12
-75 -10.8262629 12.0695095 0.581357 0.084547 19.10402 1.5000 8.08
-75 10.6503513 -12.0025208 0.473102 0.082742 29.42353 1.5000 8.33
-76 1.1545882 -1.0972397 0.741401 0.013959 30.83257 1.5000 11.72
-76 -0.6725456 -3.0000851 0.466250 0.018268 64.73455 1.5000 11.93
-76 -15.2962124 15.3637616 0.236516 0.012636 0.67847 1.5000 13.07
-76 6.1748702 -7.1198104 0.347677 0.024568 175.97714 1.5000 11.93
-77 -2.6973417 0.2940260 0.482623 0.019284 42.39270 1.5000 9.32
-77 -10.6855249 10.9131743 0.241891 0.018527 161.38476 1.5000 10.12
-77 10.6884502 -13.1698494 0.170452 0.017441 157.84230 1.5000 10.56
-78 2.3094589 17.6823759 1.660124 0.394168 164.03633 1.5000 11.59
-78 9.2854380 13.3080509 2.051221 0.492479 140.77914 1.5000 11.12
-78 17.2189593 3.0433873 1.046768 0.369445 86.13514 1.5000 12.16
-79 16.7834979 -1.7923174 2.029341 0.021093 80.57855 1.5000 11.19
-79 1.6352172 16.1829209 1.262533 0.108044 10.85484 1.5000 9.93
-79 5.2103675 14.2777820 0.777049 0.179898 35.13864 1.5000 9.90
-80 -14.5875156 10.8515971 0.586897 0.165819 124.75138 1.5000 6.65
-80 -7.6472594 -3.8905882 0.650377 0.190278 138.08826 1.5000 6.39
-80 6.9882215 -13.5565861 0.532227 0.226926 129.66115 1.5000 6.42
-81 2.6039489 3.6745080 2.930790 0.046937 31.86086 1.5000 11.45
-81 10.7289822 -11.1572676 1.939554 0.057465 58.34173 1.5000 11.68
-81 -10.7738999 12.8109455 1.862324 0.067971 46.64140 1.5000 11.54
+1 2.6786414 2.9037137 34.846746 0.044513 62.43569 1.5000 7.31
+1 -11.8484678 13.3774063 41.530243 0.058299 177.79654 1.5000 6.83
+1 4.0661877 3.9744763 29.457874 0.062625 88.13611 1.5000 7.12
+1 9.9037156 -10.4750928 34.634902 0.075692 176.98499 1.5000 6.74
+1 4.0126317 3.7839766 8.056275 0.092167 29.31648 1.5000 8.11
+2 -2.6672578 18.3345087 41.997267 0.110967 173.97058 1.5000 6.12
+2 10.9447278 9.1407162 21.764719 0.201950 53.10684 1.5000 6.18
+2 15.5032593 1.7617723 21.466000 0.359263 2.07308 1.5000 5.57
+3 2.5286812 2.1814455 8.492000 0.162661 47.99553 1.5000 7.44
+3 -12.8723561 13.6577943 6.819286 0.276050 22.52507 1.5000 7.10
+3 9.0081277 -10.0961549 6.817756 0.353202 31.25532 1.5000 6.84
+4 0.9931065 -0.6166611 30.665228 0.055918 49.21350 1.5000 7.21
+4 2.6731229 0.5936478 29.130127 0.064254 22.85212 1.5000 7.11
+4 -1.3357837 -2.8791104 23.057323 0.057944 27.85221 1.5000 7.48
+4 6.8406945 -8.1120884 20.574529 0.084858 91.24217 1.5000 7.19
+4 -14.8580101 14.8551082 11.821152 0.067028 86.10016 1.5000 8.04
+5 3.0332009 -14.6506316 44.886755 0.021107 16.70178 1.5000 7.85
+5 -16.7625927 10.1664732 28.244146 0.021254 39.18535 1.5000 8.35
+5 -8.4051550 -5.9872594 23.055451 0.027855 44.51158 1.5000 8.27
+6 -11.5124820 15.9970530 40.696426 0.013649 36.64913 1.5000 8.43
+6 10.1123828 -6.2419359 56.155555 0.018888 63.79066 1.5000 7.73
+6 6.8681008 4.5905394 38.939423 0.020371 11.59956 1.5000 8.05
+7 24.4507881 35.4519041 12.030865 0.300848 1.46355 1.5000 6.39
+8 0.4294985 0.3368270 17.875083 0.032617 41.79900 1.5000 8.38
+8 1.6975175 1.4121062 17.112202 0.036880 38.61027 1.5000 8.30
+8 -2.6232074 -2.2017016 12.024122 0.043511 37.79858 1.5000 8.50
+8 8.4787047 -10.5710816 7.510928 0.077451 57.12033 1.5000 8.38
+8 -13.4279371 13.3064944 6.248360 0.071115 48.69075 1.5000 8.68
+9 -12.2729411 23.8607339 20.470974 0.136267 2.27812 1.5000 6.68
+10 8.3810986 7.9105612 16.010085 0.133434 131.96675 1.5000 6.96
+10 13.6029426 -3.8834435 17.642796 0.159210 104.80731 1.5000 6.66
+10 -6.8153894 16.9665646 9.442848 0.187264 133.56470 1.5000 7.17
+11 -5.1498533 15.3209108 173.625956 0.031134 14.17097 1.5000 5.95
+11 14.4489784 -7.0793624 140.119597 0.032200 49.70945 1.5000 6.15
+11 5.5393836 10.3171987 119.187919 0.048467 84.58593 1.5000 5.88
+13 7.1884621 14.6438184 8.810482 0.197648 149.64280 1.5000 7.19
+13 3.4752171 16.8591921 5.944757 0.243989 146.03286 1.5000 7.39
+13 17.4302447 2.1956866 2.643565 0.166403 100.26995 1.5000 8.68
+15 15.5981812 -6.8318188 25.200594 0.009008 63.02757 1.5000 9.40
+15 -2.4643003 15.0794959 31.326621 0.012990 15.38269 1.5000 8.77
+15 3.9723885 12.3621230 21.562778 0.021541 103.31672 1.5000 8.63
+18 -30.7098819 -12.9315068 6.970062 0.009213 81.11684 1.5000 10.77
+19 3.6379252 -17.6227227 7.146992 0.110295 19.29253 1.5000 8.04
+19 -16.0320557 5.4916474 5.593125 0.221276 35.01589 1.5000 7.55
+19 -11.8452654 -3.3467879 4.615328 0.203070 43.22065 1.5000 7.85
+21 -20.5880555 8.9330217 0.336281 0.012562 64.84746 1.5000 9.16
+21 -6.9429627 -12.0159266 1.166746 0.070850 153.61893 1.5000 5.93
+21 -6.0959847 -12.6680059 0.871997 0.092417 127.23807 1.5000 5.96
+22 17.3720423 2.1882518 0.401602 0.028619 88.59134 1.5000 7.71
+22 7.5287043 14.3642186 0.621280 0.062720 157.01625 1.5000 6.38
+22 3.0851649 17.0240832 0.263579 0.118444 129.67642 1.5000 6.62
+23 2.8607712 -0.5076926 0.213277 0.013634 42.04891 1.5000 11.25
+23 1.6688003 -1.2566873 0.211058 0.014425 24.74250 1.5000 11.20
+23 -15.6866333 15.6846331 0.055011 0.007139 161.67661 1.5000 13.42
+23 5.5056116 -6.4105492 0.103301 0.014270 154.36472 1.5000 11.99
+23 -0.1804322 -3.2186461 0.079899 0.021747 92.27222 1.5000 11.81
+24 16.5900066 -0.7754914 0.712519 0.072540 87.02948 1.5000 8.20
+24 7.3401048 12.9906345 0.760423 0.122910 154.27397 1.5000 7.56
+24 0.0386506 17.0817865 0.360577 0.204579 144.39844 1.5000 7.82
+25 -0.5583018 18.1080905 0.501994 0.040558 170.87683 1.5000 12.51
+25 10.3473234 11.1996532 0.595960 0.072388 123.22053 1.5000 11.69
+25 16.4071453 2.0674751 0.195809 0.113742 29.24962 1.5000 12.41
+26 -3.6567136 -0.9732294 0.613570 0.136539 20.52350 1.5000 8.39
+26 9.7897822 -12.7521828 0.387438 0.165176 60.03851 1.5000 8.68
+26 -11.8833211 11.3909941 0.388550 0.197227 54.76414 1.5000 8.49
+28 -8.8303495 -10.2295939 6.278282 0.234335 139.37523 1.5000 6.45
+28 -4.0458078 -14.0620768 4.997050 0.265757 137.81228 1.5000 6.56
+28 -20.3075095 8.7187072 1.783383 0.229330 104.42889 1.5000 7.84
+29 -15.4000922 14.6578471 0.554882 0.080878 168.20716 1.5000 7.98
+29 6.1406271 -8.3810986 0.850477 0.127398 160.62243 1.5000 7.02
+29 -1.9455793 -3.6482572 0.607310 0.104516 75.02961 1.5000 7.60
+30 16.9039460 -0.7408146 0.359729 0.006926 79.75543 1.5000 10.59
+30 1.7452270 16.5074530 0.192209 0.032517 12.43361 1.5000 9.59
+30 5.9670452 14.1843178 0.142793 0.047678 18.69459 1.5000 9.50
+31 16.5139533 -3.8451131 0.953037 0.000768 78.65420 1.5000 13.92
+31 1.3122749 15.4732685 0.736143 0.005785 8.96780 1.5000 12.01
+31 3.6766341 14.3003407 0.396902 0.009968 50.91748 1.5000 12.09
+32 -19.3763875 6.4541278 0.128487 0.031001 54.30275 1.5000 10.09
+32 -1.0479213 -16.8624332 0.128406 0.034213 18.77784 1.5000 9.98
+32 -11.6234921 -7.1394340 0.093792 0.032315 27.56423 1.5000 10.38
+33 4.9678189 6.5523306 1.035937 0.007901 58.63553 1.5000 12.89
+33 4.8916122 6.4043332 0.134215 0.001541 67.64338 1.5000 16.89
+33 4.4835950 5.8871981 0.940972 0.014050 53.91057 1.5000 12.37
+33 11.5214718 -9.5055470 1.026647 0.019262 61.39067 1.5000 11.94
+33 -9.7698615 14.0765994 0.908140 0.022506 44.63464 1.5000 11.90
+34 -2.2449807 0.8443678 2.901961 0.044099 35.99226 1.5000 11.38
+34 -1.7488676 1.1591862 2.977897 0.049058 29.24244 1.5000 11.24
+34 -10.3350314 10.6553444 0.528620 0.034140 157.42648 1.5000 13.51
+34 10.9154015 -13.3958844 0.355410 0.030735 151.89374 1.5000 14.05
+34 -0.1693480 2.6923842 0.460071 0.085692 92.53445 1.5000 12.66
+35 4.9563019 -16.6557199 1.930376 0.327434 22.06554 1.5000 7.96
+35 -15.3681458 7.0270758 1.355189 0.533500 46.15385 1.5000 7.81
+35 -10.6478399 -3.3934448 1.012171 0.499596 26.32483 1.5000 8.20
+36 3.2872076 -12.0751123 14.331348 0.001907 13.82550 1.5000 11.15
+36 1.6738457 -11.2374340 15.677721 0.002940 137.74766 1.5000 10.58
+36 3.1567351 -11.1001942 0.296112 0.000065 147.14197 1.5000 19.02
+36 -5.1645512 -6.8981440 6.859944 0.001718 126.59146 1.5000 12.06
+36 2.7773357 -10.1800027 9.300518 0.002411 40.70501 1.5000 11.36
+36 4.6978370 -10.8714327 10.212730 0.003100 139.05250 1.5000 10.99
+36 -16.6487889 13.0453302 5.206013 0.001663 156.29122 1.5000 12.40
+37 -2.1269981 17.7617793 0.839156 0.146064 4.33835 1.5000 10.16
+37 15.9185386 -0.2597829 0.821566 0.184455 67.47268 1.5000 9.93
+37 9.4646705 10.7235961 0.648079 0.368730 111.52299 1.5000 9.44
+38 -6.2209643 -15.8705378 0.937813 0.024682 150.95781 1.5000 12.26
+38 -13.2037550 -10.2148779 0.700964 0.031005 127.62004 1.5000 12.33
+38 -22.1953978 5.4721115 0.242842 0.040746 173.16383 1.5000 13.18
+39 3.0311854 -15.8087389 1.252695 0.234004 3.24979 1.5000 8.95
+39 -9.7784399 -5.3083389 0.738281 0.286702 83.26631 1.5000 9.30
+39 -16.7501370 8.6159058 0.771763 0.333544 8.36668 1.5000 9.09
+40 14.4943918 -4.2787566 0.899056 0.075603 70.33570 1.5000 7.07
+40 -5.2517705 16.7018885 0.654067 0.082554 26.17358 1.5000 7.32
+40 7.7887960 9.4447620 0.461531 0.148342 29.43871 1.5000 7.06
+41 -21.3994996 7.4374779 0.805568 0.051520 45.86021 1.5000 10.84
+41 -6.4802050 -13.9904609 1.839320 0.135517 163.28270 1.5000 8.89
+41 -9.5515104 -11.4940659 1.284799 0.250315 107.79707 1.5000 8.61
+42 -13.1218077 -2.4426109 2.740340 0.076499 108.90562 1.5000 7.35
+42 -15.9520107 3.4705471 2.675128 0.095901 126.65994 1.5000 7.13
+42 3.2708949 -18.5152726 1.415411 0.073095 157.24523 1.5000 8.12
+43 2.4040622 -15.7223625 1.788496 0.078973 3.45353 1.5000 11.06
+43 -17.1344297 8.7164684 1.119392 0.106268 14.08819 1.5000 11.25
+43 -9.7153852 -5.7390496 1.035435 0.100594 75.58401 1.5000 11.39
+44 6.0645759 12.9687435 3.176308 0.076161 149.54933 1.5000 9.26
+44 16.2815141 -2.9996537 1.610420 0.078817 101.50110 1.5000 9.96
+44 -0.7771551 16.5017779 2.063999 0.102766 149.23566 1.5000 9.40
+45 9.4626240 -6.2544876 1.954645 0.014537 64.97744 1.5000 11.32
+45 -12.2638492 15.9438916 1.317881 0.010645 39.84224 1.5000 12.09
+45 6.5910318 3.9202741 1.399586 0.015605 4.33088 1.5000 11.61
+46 -3.0687338 -17.7236024 0.356506 0.042735 171.67764 1.5000 9.74
+46 -13.9649554 -8.2787881 0.235762 0.056968 21.99600 1.5000 9.88
+46 -20.6025935 4.0245670 0.274450 0.072585 38.35190 1.5000 9.45
+47 7.8688735 -16.4945455 2.140230 0.052530 36.56310 1.5000 5.67
+47 -13.2316933 6.5044726 2.135390 0.089221 58.89319 1.5000 5.10
+47 -9.9045947 -1.3440890 2.118815 0.101371 171.47927 1.5000 4.97
+48 -11.8766946 -5.4321795 1.277907 0.000346 122.06413 1.5000 6.19
+48 -17.8993441 5.9518384 1.256311 0.000450 122.67094 1.5000 5.92
+48 1.1980351 -17.4474217 1.104809 0.000408 143.89240 1.5000 6.17
+49 15.1040854 2.3571428 20.405080 0.181247 112.06946 1.5000 5.63
+49 -3.2499829 18.5619013 9.661079 0.185628 143.94899 1.5000 6.41
+49 11.0148414 8.7776491 7.289542 0.142045 128.95513 1.5000 7.01
+50 4.0302499 1.6646014 1.472434 0.043041 52.02926 1.5000 11.25
+50 -14.3130702 15.3193587 0.715350 0.038570 19.59706 1.5000 12.15
+50 7.4496735 -7.3164411 0.854058 0.072139 27.68499 1.5000 11.28
+51 -3.3306276 15.3707635 1.049646 0.198237 3.90516 1.5000 8.52
+51 4.9105775 11.6912780 0.876254 0.255630 120.78365 1.5000 8.44
+51 15.2979283 -6.5626243 0.655293 0.224586 52.14704 1.5000 8.90
+53 12.4238186 -9.5226407 2.782816 0.158592 59.54959 1.5000 5.89
+53 4.4818211 7.4433151 4.129262 0.264990 70.81067 1.5000 4.90
+53 -8.4071710 13.9983285 2.582502 0.205370 39.71256 1.5000 5.69
+55 13.7062119 1.7502342 5.328491 0.056648 109.47776 1.5000 6.36
+55 10.8979365 7.0403753 3.840084 0.048798 128.49434 1.5000 6.88
+55 -5.6997638 18.4884777 1.410663 0.058820 140.27986 1.5000 7.76
+56 6.0934996 -16.6090512 1.581300 0.092701 33.50068 1.5000 7.39
+56 -14.6185041 6.8228040 1.472465 0.122278 63.75571 1.5000 7.17
+56 -10.5003167 -2.6646358 1.105217 0.117459 175.23463 1.5000 7.52
+57 16.5121298 -3.5196512 0.091809 0.001033 78.06134 1.5000 12.52
+57 0.8380873 15.7963787 0.075625 0.004756 7.95489 1.5000 11.07
+57 4.3549332 14.0322052 0.037846 0.008974 53.77328 1.5000 11.13
+58 6.1940435 7.6041197 3.148249 0.781489 148.17202 1.5000 5.26
+58 12.4840475 -7.6487274 3.288870 0.829461 100.17624 1.5000 5.15
+58 -8.4859894 15.2787222 2.403389 1.040798 125.67266 1.5000 5.24
+59 -1.3843044 0.2825942 0.978418 0.055405 37.83465 1.5000 12.31
+59 -2.1529569 -0.2200877 0.944927 0.060551 28.59354 1.5000 12.25
+59 1.2664764 2.7720168 0.338380 0.058125 35.35961 1.5000 13.41
+59 10.2287475 -12.1804134 0.215495 0.074500 76.11086 1.5000 13.63
+59 -11.3728101 11.9229549 0.200815 0.095614 73.42242 1.5000 13.44
+60 -24.2385721 12.1078788 1.845691 0.043224 19.30370 1.5000 5.89
+62 1.6683259 3.3763712 0.732983 0.064828 53.95580 1.5000 8.12
+62 -10.8262474 12.0695496 0.581357 0.084546 19.10396 1.5000 8.08
+62 10.6503443 -12.0025523 0.473102 0.082742 29.42349 1.5000 8.33
+63 -18.3247480 11.1333391 0.506447 0.006385 64.26762 1.5000 12.98
+63 -6.8051055 -8.5551042 0.515073 0.012408 8.22585 1.5000 12.24
+63 -0.6725802 -13.0001017 0.466238 0.018267 64.73402 1.5000 11.93
+64 -2.6972584 0.2939139 0.482636 0.019284 42.39216 1.5000 9.32
+64 -1.1736250 1.2720707 0.550505 0.024027 27.21757 1.5000 8.94
+64 -10.6854789 10.9131908 0.241892 0.018527 161.38484 1.5000 10.12
+64 10.6884103 -13.1698056 0.170452 0.017441 157.84233 1.5000 10.56
+64 0.0847059 2.5684164 0.261767 0.035008 78.77792 1.5000 9.34
+65 2.3094964 17.6823656 1.660138 0.394169 164.03616 1.5000 11.59
+65 9.2855248 13.3080092 2.051222 0.492480 140.77908 1.5000 11.12
+65 17.2189383 3.0433702 1.046772 0.369444 86.13513 1.5000 12.16
+66 16.7835267 -1.7922896 2.029349 0.021093 80.57857 1.5000 11.19
+66 1.6353393 16.1830764 1.262461 0.108044 10.85672 1.5000 9.93
+66 5.2102839 14.2779135 0.777039 0.179918 35.14046 1.5000 9.90
+68 -14.5875060 10.8516456 0.586895 0.165819 124.75149 1.5000 6.65
+68 -7.6472299 -3.8905736 0.650382 0.190282 138.08841 1.5000 6.39
+68 6.9882137 -13.5564786 0.532229 0.226926 129.66042 1.5000 6.42
+69 2.6039461 3.6745087 2.930789 0.046937 31.86090 1.5000 11.45
+69 10.7290626 -11.1572633 1.939552 0.057464 58.34179 1.5000 11.68
+69 -10.7739067 12.8109934 1.862320 0.067971 46.64136 1.5000 11.54
+70 6.0934996 -16.6090512 1.581300 0.092701 33.50068 1.5000 7.39
+70 -14.6185041 6.8228040 1.472465 0.122278 63.75571 1.5000 7.17
+70 -10.5003167 -2.6646358 1.105217 0.117459 175.23463 1.5000 7.52
+71 16.5121298 -3.5196512 0.091809 0.001033 78.06134 1.5000 12.52
+71 0.8380873 15.7963787 0.075625 0.004756 7.95489 1.5000 11.07
+71 4.3549332 14.0322052 0.037846 0.008974 53.77328 1.5000 11.13
+72 6.1940435 7.6041197 3.148249 0.781489 148.17202 1.5000 5.26
+72 12.4840475 -7.6487274 3.288870 0.829461 100.17624 1.5000 5.15
+72 -8.4859894 15.2787222 2.403389 1.040798 125.67266 1.5000 5.24
+73 -1.3843044 0.2825942 0.978418 0.055405 37.83465 1.5000 12.31
+73 -2.1529569 -0.2200877 0.944927 0.060551 28.59354 1.5000 12.25
+73 1.2664764 2.7720168 0.338380 0.058125 35.35961 1.5000 13.41
+73 10.2287475 -12.1804134 0.215495 0.074500 76.11086 1.5000 13.63
+73 -11.3728101 11.9229549 0.200815 0.095614 73.42242 1.5000 13.44
+74 -40.2385139 12.1078204 1.845692 0.043224 19.30365 1.5000 5.89
+75 1.6683259 3.3763712 0.732983 0.064828 53.95580 1.5000 8.12
+75 -10.8262474 12.0695496 0.581357 0.084546 19.10396 1.5000 8.08
+75 10.6503443 -12.0025523 0.473102 0.082742 29.42349 1.5000 8.33
+76 5.4470051 -16.7065589 0.821138 0.012368 8.82319 1.5000 11.74
+76 -10.6725714 -3.0001379 0.466236 0.018267 64.73442 1.5000 11.93
+76 -15.0255834 6.8065365 0.504194 0.023781 2.47909 1.5000 11.56
+77 -2.6972584 0.2939139 0.482636 0.019284 42.39216 1.5000 9.32
+77 -1.1736250 1.2720707 0.550505 0.024027 27.21757 1.5000 8.94
+77 -10.6854789 10.9131908 0.241892 0.018527 161.38484 1.5000 10.12
+77 10.6884103 -13.1698056 0.170452 0.017441 157.84233 1.5000 10.56
+77 0.0847059 2.5684164 0.261767 0.035008 78.77792 1.5000 9.34
+78 2.3094964 17.6823656 1.660138 0.394169 164.03616 1.5000 11.59
+78 9.2855248 13.3080092 2.051222 0.492480 140.77908 1.5000 11.12
+78 17.2189383 3.0433702 1.046772 0.369444 86.13513 1.5000 12.16
+79 16.7835267 -1.7922896 2.029349 0.021093 80.57857 1.5000 11.19
+79 1.6353393 16.1830764 1.262461 0.108044 10.85672 1.5000 9.93
+79 5.2102839 14.2779135 0.777039 0.179918 35.14046 1.5000 9.90
+80 -14.5875060 10.8516456 0.586895 0.165819 124.75149 1.5000 6.65
+80 -7.6472299 -3.8905736 0.650382 0.190282 138.08841 1.5000 6.39
+80 6.9882137 -13.5564786 0.532229 0.226926 129.66042 1.5000 6.42
+81 2.6039461 3.6745087 2.930789 0.046937 31.86090 1.5000 11.45
+81 10.7290626 -11.1572633 1.939552 0.057464 58.34179 1.5000 11.68
+81 -10.7739067 12.8109934 1.862320 0.067971 46.64136 1.5000 11.54
diff --git a/Benchmarks/ConfigFiles/m1931_modif_210Img.txt b/Benchmarks/ConfigFiles/m1931_modif_210Img_sec.txt
similarity index 90%
copy from Benchmarks/ConfigFiles/m1931_modif_210Img.txt
copy to Benchmarks/ConfigFiles/m1931_modif_210Img_sec.txt
index 535aeac..e174318 100644
--- a/Benchmarks/ConfigFiles/m1931_modif_210Img.txt
+++ b/Benchmarks/ConfigFiles/m1931_modif_210Img_sec.txt
@@ -1,208 +1,201 @@
#REFERENCE 0 0.0000000 0.0000000
1 2.6786527 2.9037603 34.846982 0.044513 62.43534 1.5000 7.31
1 -11.8484881 13.3774383 41.530182 0.058299 177.79656 1.5000 6.83
1 9.9036907 -10.4750539 34.634955 0.075692 176.98494 1.5000 6.74
2 -2.6672309 18.3345021 41.997432 0.110967 173.97052 1.5000 6.12
2 10.9447338 9.1407289 21.765016 0.201954 53.10618 1.5000 6.18
2 15.5033588 1.7617761 21.466066 0.359257 2.07257 1.5000 5.57
3 2.5286856 2.1814895 8.491896 0.162661 47.99557 1.5000 7.44
3 -12.8724137 13.6578294 6.819269 0.276049 22.52509 1.5000 7.10
3 9.0081475 -10.0961433 6.817751 0.353202 31.25535 1.5000 6.84
4 -1.3358012 -2.8790666 23.057318 0.057944 27.85238 1.5000 7.48
4 6.8407717 -8.1121706 20.574407 0.084857 91.24205 1.5000 7.19
4 -14.8580323 14.8551029 11.821152 0.067028 86.10019 1.5000 8.04
5 3.0331801 -14.6505707 44.887170 0.021107 16.70145 1.5000 7.85
5 -16.7625653 10.1664425 28.244176 0.021254 39.18539 1.5000 8.35
5 -8.4052010 -5.9872161 23.055401 0.027855 44.51135 1.5000 8.27
6 -11.5124223 15.9970478 40.696439 0.013649 36.64909 1.5000 8.43
6 10.1123445 -6.2419290 56.155712 0.018888 63.79062 1.5000 7.73
6 6.8681360 4.5905742 38.939242 0.020371 11.59938 1.5000 8.05
7 24.4507222 35.4518592 12.030922 0.300848 1.46343 1.5000 6.39
8 0.4294637 0.3368837 17.875325 0.032617 41.79915 1.5000 8.38
8 1.6975484 1.4120697 17.112294 0.036880 38.61008 1.5000 8.30
8 -2.6231902 -2.2017254 12.024044 0.043511 37.79877 1.5000 8.50
8 8.4786219 -10.5709866 7.510960 0.077451 57.12046 1.5000 8.38
8 -13.4278497 13.3064510 6.248381 0.071115 48.69072 1.5000 8.68
9 -12.2729271 23.8607406 20.470970 0.136267 2.27813 1.5000 6.68
10 8.3811315 7.9105283 16.010105 0.133434 131.96659 1.5000 6.96
10 13.6029638 -3.8834001 17.642846 0.159210 104.80730 1.5000 6.66
10 -6.8153537 16.9665578 9.442873 0.187264 133.56479 1.5000 7.17
11 -5.1498385 15.3208303 173.626665 0.031134 14.17090 1.5000 5.95
11 14.4490163 -7.0793369 140.119512 0.032200 49.70954 1.5000 6.15
11 5.5394577 10.3171621 119.187491 0.048467 84.58532 1.5000 5.88
-13 7.1884598 14.6438521 8.810685 0.197648 149.64275 1.5000 7.19
-13 3.4752244 16.8592447 5.944581 0.243989 146.03268 1.5000 7.39
-13 17.4301938 2.1957167 2.643589 0.166403 100.27012 1.5000 8.68
+12 7.1884598 14.6438521 8.810685 0.197648 149.64275 1.5000 7.19
+12 3.4752244 16.8592447 5.944581 0.243989 146.03268 1.5000 7.39
+12 17.4301938 2.1957167 2.643589 0.166403 100.27012 1.5000 8.68
15 15.5982332 -6.8318229 25.200524 0.009008 63.02755 1.5000 9.40
15 -2.4642730 15.0794757 31.326845 0.012990 15.38250 1.5000 8.77
15 3.9724585 12.3620713 21.562202 0.021541 103.31489 1.5000 8.63
-18 1.0479069 -1.1044685 10.768333 0.006847 22.52448 1.5000 10.62
-18 -15.2217938 15.3201427 5.277769 0.004422 164.44813 1.5000 11.87
-18 6.2890313 -7.2140096 8.200618 0.007841 155.73619 1.5000 10.77
-18 -0.7098203 -2.9315479 6.970110 0.009213 81.11879 1.5000 10.77
19 3.6379121 -17.6226603 7.147027 0.110295 19.29235 1.5000 8.04
19 -16.0320537 5.4916741 5.593126 0.221276 35.01584 1.5000 7.55
19 -11.8452915 -3.3468708 4.615314 0.203069 43.22132 1.5000 7.85
21 -20.5880303 8.9329582 0.336281 0.012562 64.84756 1.5000 9.16
22 17.3720134 2.1882537 0.401604 0.028619 88.59146 1.5000 7.71
22 7.5287582 14.3644280 0.621351 0.062721 157.01493 1.5000 6.38
22 3.0852536 17.0240154 0.263595 0.118443 129.67801 1.5000 6.62
23 -15.6866276 15.6845567 0.055011 0.007139 161.67654 1.5000 13.42
23 5.5055808 -6.4106321 0.103302 0.014270 154.36511 1.5000 11.99
23 -0.1804981 -3.2185423 0.079898 0.021746 92.26818 1.5000 11.81
24 16.5899706 -0.7754639 0.712518 0.072540 87.02951 1.5000 8.20
24 7.3400718 12.9906877 0.760430 0.122910 154.27397 1.5000 7.56
24 0.0386740 17.0817512 0.360582 0.204579 144.39873 1.5000 7.82
25 -0.5582720 18.1081624 0.501990 0.040558 170.87681 1.5000 12.51
25 10.3472270 11.1997347 0.595960 0.072387 123.22072 1.5000 11.69
25 16.4070668 2.0675156 0.195808 0.113744 29.24929 1.5000 12.41
26 -3.6566538 -0.9732215 0.613580 0.136538 20.52408 1.5000 8.39
26 9.7897135 -12.7521977 0.387439 0.165176 60.03832 1.5000 8.68
26 -11.8832764 11.3910034 0.388550 0.197227 54.76396 1.5000 8.49
28 -8.8301584 -10.2296004 6.278205 0.234335 139.37575 1.5000 6.45
28 -4.0456315 -14.0621632 4.996734 0.265758 137.81187 1.5000 6.56
28 -20.3074952 8.7186627 1.783389 0.229330 104.42894 1.5000 7.84
29 -15.4001065 14.6578526 0.554882 0.080878 168.20715 1.5000 7.98
29 6.1406896 -8.3811360 0.850465 0.127398 160.62235 1.5000 7.02
29 -1.9455502 -3.6483147 0.607322 0.104519 75.02974 1.5000 7.60
30 16.9039033 -0.7407439 0.359729 0.006926 79.75539 1.5000 10.59
30 1.7452545 16.5076313 0.192199 0.032516 12.43573 1.5000 9.59
30 5.9668771 14.1843117 0.142788 0.047680 18.69737 1.5000 9.50
31 16.5139312 -3.8451362 0.953037 0.000768 78.65420 1.5000 13.92
31 1.3123810 15.4732942 0.736126 0.005786 8.96831 1.5000 12.01
31 3.6765900 14.3002919 0.396902 0.009968 50.91878 1.5000 12.09
32 -19.3763552 6.4541388 0.128487 0.031001 54.30277 1.5000 10.09
32 -1.0478628 -16.8624882 0.128406 0.034213 18.77814 1.5000 9.98
32 -11.6234609 -7.1394199 0.093792 0.032315 27.56403 1.5000 10.38
33 4.4835691 5.8872015 0.940979 0.014051 53.91017 1.5000 12.37
33 11.5215383 -9.5055949 1.026644 0.019262 61.39061 1.5000 11.94
33 -9.7699119 14.0765786 0.908140 0.022506 44.63471 1.5000 11.90
34 -10.3351451 10.6553738 0.528612 0.034140 157.42627 1.5000 13.51
34 10.9153487 -13.3958584 0.355411 0.030735 151.89381 1.5000 14.05
34 -0.1692590 2.6923973 0.460052 0.085688 92.53567 1.5000 12.66
35 4.9563344 -16.6556871 1.930383 0.327434 22.06552 1.5000 7.96
35 -15.3680550 7.0269874 1.355196 0.533507 46.15432 1.5000 7.81
35 -10.6479548 -3.3933905 1.012174 0.499609 26.32416 1.5000 8.20
+36 -5.1645974 -6.8981448 6.859959 0.001718 126.59104 1.5000 12.06
+36 4.6977819 -10.8714692 10.212929 0.003100 139.05114 1.5000 10.99
+36 -16.6488044 13.0453261 5.206010 0.001663 156.29120 1.5000 12.40
37 -2.1269883 17.7618777 0.839149 0.146064 4.33839 1.5000 10.16
37 15.9185406 -0.2597843 0.821566 0.184455 67.47268 1.5000 9.93
37 9.4646689 10.7235580 0.648071 0.368732 111.52255 1.5000 9.44
38 -6.2210215 -15.8704320 0.937844 0.024682 150.95757 1.5000 12.26
38 -13.2038524 -10.2148699 0.700936 0.031004 127.61932 1.5000 12.33
38 -22.1953724 5.4720839 0.242843 0.040746 173.16403 1.5000 13.18
39 3.0312580 -15.8087360 1.252698 0.234004 3.25013 1.5000 8.95
39 -9.7784925 -5.3082813 0.738278 0.286702 83.26595 1.5000 9.30
39 -16.7500612 8.6158344 0.771764 0.333547 8.36658 1.5000 9.09
40 14.4943275 -4.2786351 0.899064 0.075604 70.33597 1.5000 7.07
40 -5.2517624 16.7019063 0.654067 0.082554 26.17356 1.5000 7.32
40 7.7887901 9.4448365 0.461530 0.148344 29.43926 1.5000 7.06
41 -21.3995271 7.4375134 0.805568 0.051520 45.86022 1.5000 10.84
41 -6.4802887 -13.9904447 1.839369 0.135518 163.28206 1.5000 8.89
41 -9.5515612 -11.4941910 1.284929 0.250311 107.80091 1.5000 8.61
42 -13.1218821 -2.4426522 2.740368 0.076499 108.90595 1.5000 7.35
42 -15.9520465 3.4705114 2.675136 0.095900 126.65988 1.5000 7.13
42 3.2709225 -18.5152424 1.415411 0.073096 157.24543 1.5000 8.12
43 2.4040917 -15.7224455 1.788478 0.078973 3.45380 1.5000 11.06
43 -17.1344846 8.7164655 1.119390 0.106268 14.08818 1.5000 11.25
43 -9.7154364 -5.7390332 1.035433 0.100593 75.58388 1.5000 11.39
44 6.0647302 12.9685543 3.176132 0.076161 149.54889 1.5000 9.26
44 16.2815642 -2.9996360 1.610419 0.078817 101.50112 1.5000 9.96
44 -0.7772337 16.5018612 2.063908 0.102766 149.23518 1.5000 9.40
45 9.4626405 -6.2544580 1.954653 0.014537 64.97761 1.5000 11.32
45 -12.2639022 15.9438784 1.317880 0.010645 39.84226 1.5000 12.09
45 6.5910853 3.9203102 1.399579 0.015605 4.33053 1.5000 11.61
46 -3.0687036 -17.7236144 0.356505 0.042735 171.67775 1.5000 9.74
46 -13.9649508 -8.2788025 0.235761 0.056968 21.99631 1.5000 9.88
46 -20.6025554 4.0245322 0.274445 0.072585 38.35050 1.5000 9.45
47 7.8688411 -16.4944952 2.140231 0.052530 36.56283 1.5000 5.67
47 -13.2317046 6.5044642 2.135389 0.089221 58.89322 1.5000 5.10
47 -9.9045632 -1.3440293 2.118868 0.101374 171.47875 1.5000 4.97
48 -11.8767244 -5.4322162 1.277910 0.000346 122.06425 1.5000 6.19
48 -17.8993462 5.9518549 1.256308 0.000450 122.67098 1.5000 5.92
48 1.1979541 -17.4473601 1.104821 0.000408 143.89249 1.5000 6.17
49 15.1040611 2.3571116 20.405094 0.181247 112.06935 1.5000 5.63
49 -3.2500200 18.5620044 9.660994 0.185627 143.94883 1.5000 6.41
49 11.0148057 8.7776669 7.289802 0.142045 128.95353 1.5000 7.01
50 4.0302890 1.6646021 1.472455 0.043042 52.02903 1.5000 11.25
50 -14.3131270 15.3193885 0.715348 0.038570 19.59704 1.5000 12.15
50 7.4497300 -7.3163268 0.854058 0.072140 27.68554 1.5000 11.28
51 -3.3306534 15.3708033 1.049641 0.198236 3.90524 1.5000 8.52
51 4.9104924 11.6912630 0.876250 0.255629 120.78366 1.5000 8.44
51 15.2979744 -6.5626779 0.655289 0.224586 52.14670 1.5000 8.90
53 12.4238247 -9.5226491 2.782814 0.158592 59.54957 1.5000 5.89
53 4.4819744 7.4432173 4.128941 0.264958 70.80977 1.5000 4.90
53 -8.4071970 13.9983475 2.582500 0.205369 39.71261 1.5000 5.69
55 13.7061891 1.7502961 5.328571 0.056648 109.47795 1.5000 6.36
55 10.8979905 7.0403978 3.840162 0.048798 128.49451 1.5000 6.88
55 -5.6997319 18.4884138 1.410669 0.058820 140.27996 1.5000 7.76
56 6.0935140 -16.6090186 1.581300 0.092701 33.50078 1.5000 7.39
56 -14.6184750 6.8227033 1.472476 0.122279 63.75621 1.5000 7.17
56 -10.5002861 -2.6646404 1.105215 0.117458 175.23476 1.5000 7.52
57 16.5121525 -3.5196819 0.091809 0.001033 78.06127 1.5000 12.52
57 0.8381572 15.7962753 0.075629 0.004756 7.95312 1.5000 11.07
57 4.3548701 14.0324012 0.037846 0.008976 53.77684 1.5000 11.13
58 6.1939611 7.6041748 3.148253 0.781497 148.17302 1.5000 5.26
58 12.4841128 -7.6487509 3.288863 0.829458 100.17621 1.5000 5.15
58 -8.4859528 15.2787147 2.403395 1.040799 125.67280 1.5000 5.24
59 1.2664786 2.7720698 0.338377 0.058125 35.35947 1.5000 13.41
59 10.2287300 -12.1803663 0.215495 0.074500 76.11091 1.5000 13.63
59 -11.3726996 11.9228981 0.200816 0.095615 73.42239 1.5000 13.44
-60 14.2660661 -10.8918824 1.461158 0.020111 48.44266 1.5000 6.97
-60 -4.2386057 12.1078222 1.845690 0.043224 19.30371 1.5000 5.89
62 1.6682920 3.3763413 0.732986 0.064828 53.95586 1.5000 8.12
62 -10.8262629 12.0695095 0.581357 0.084547 19.10402 1.5000 8.08
62 10.6503513 -12.0025208 0.473102 0.082742 29.42353 1.5000 8.33
-63 1.1545882 -1.0972397 0.741401 0.013959 30.83257 1.5000 11.72
-63 -0.6725456 -3.0000851 0.466250 0.018268 64.73455 1.5000 11.93
-63 -15.2962124 15.3637616 0.236516 0.012636 0.67847 1.5000 13.07
-63 6.1748702 -7.1198104 0.347677 0.024568 175.97714 1.5000 11.93
+63 -18.3247786 11.1333764 0.506446 0.006385 64.26758 1.5000 12.98
+63 -6.8051024 -8.5551780 0.515075 0.012408 8.22559 1.5000 12.24
+63 -0.6726177 -13.0001072 0.466237 0.018267 64.73408 1.5000 11.93
64 -2.6973417 0.2940260 0.482623 0.019284 42.39270 1.5000 9.32
64 -10.6855249 10.9131743 0.241891 0.018527 161.38476 1.5000 10.12
64 10.6884502 -13.1698494 0.170452 0.017441 157.84230 1.5000 10.56
65 2.3094589 17.6823759 1.660124 0.394168 164.03633 1.5000 11.59
65 9.2854380 13.3080509 2.051221 0.492479 140.77914 1.5000 11.12
65 17.2189593 3.0433873 1.046768 0.369445 86.13514 1.5000 12.16
66 16.7834979 -1.7923174 2.029341 0.021093 80.57855 1.5000 11.19
66 1.6352172 16.1829209 1.262533 0.108044 10.85484 1.5000 9.93
66 5.2103675 14.2777820 0.777049 0.179898 35.13864 1.5000 9.90
68 -14.5875156 10.8515971 0.586897 0.165819 124.75138 1.5000 6.65
68 -7.6472594 -3.8905882 0.650377 0.190278 138.08826 1.5000 6.39
68 6.9882215 -13.5565861 0.532227 0.226926 129.66115 1.5000 6.42
69 2.6039489 3.6745080 2.930790 0.046937 31.86086 1.5000 11.45
69 10.7289822 -11.1572676 1.939554 0.057465 58.34173 1.5000 11.68
69 -10.7738999 12.8109455 1.862324 0.067971 46.64140 1.5000 11.54
70 6.0935140 -16.6090186 1.581300 0.092701 33.50078 1.5000 7.39
70 -14.6184750 6.8227033 1.472476 0.122279 63.75621 1.5000 7.17
70 -10.5002861 -2.6646404 1.105215 0.117458 175.23476 1.5000 7.52
71 16.5121525 -3.5196819 0.091809 0.001033 78.06127 1.5000 12.52
71 0.8381572 15.7962753 0.075629 0.004756 7.95312 1.5000 11.07
71 4.3548701 14.0324012 0.037846 0.008976 53.77684 1.5000 11.13
72 6.1939611 7.6041748 3.148253 0.781497 148.17302 1.5000 5.26
72 12.4841128 -7.6487509 3.288863 0.829458 100.17621 1.5000 5.15
72 -8.4859528 15.2787147 2.403395 1.040799 125.67280 1.5000 5.24
73 1.2664786 2.7720698 0.338377 0.058125 35.35947 1.5000 13.41
73 10.2287300 -12.1803663 0.215495 0.074500 76.11091 1.5000 13.63
73 -11.3726996 11.9228981 0.200816 0.095615 73.42239 1.5000 13.44
-74 14.2660661 -10.8918824 1.461158 0.020111 48.44266 1.5000 6.97
-74 -4.2386057 12.1078222 1.845690 0.043224 19.30371 1.5000 5.89
75 1.6682920 3.3763413 0.732986 0.064828 53.95586 1.5000 8.12
75 -10.8262629 12.0695095 0.581357 0.084547 19.10402 1.5000 8.08
75 10.6503513 -12.0025208 0.473102 0.082742 29.42353 1.5000 8.33
-76 1.1545882 -1.0972397 0.741401 0.013959 30.83257 1.5000 11.72
-76 -0.6725456 -3.0000851 0.466250 0.018268 64.73455 1.5000 11.93
-76 -15.2962124 15.3637616 0.236516 0.012636 0.67847 1.5000 13.07
-76 6.1748702 -7.1198104 0.347677 0.024568 175.97714 1.5000 11.93
+76 5.4470677 -16.7065330 0.821139 0.012368 8.82318 1.5000 11.74
+76 -10.6726371 -3.0000599 0.466240 0.018267 64.73439 1.5000 11.93
+76 -15.0255655 6.8065642 0.504194 0.023781 2.47917 1.5000 11.56
77 -2.6973417 0.2940260 0.482623 0.019284 42.39270 1.5000 9.32
77 -10.6855249 10.9131743 0.241891 0.018527 161.38476 1.5000 10.12
77 10.6884502 -13.1698494 0.170452 0.017441 157.84230 1.5000 10.56
78 2.3094589 17.6823759 1.660124 0.394168 164.03633 1.5000 11.59
78 9.2854380 13.3080509 2.051221 0.492479 140.77914 1.5000 11.12
78 17.2189593 3.0433873 1.046768 0.369445 86.13514 1.5000 12.16
79 16.7834979 -1.7923174 2.029341 0.021093 80.57855 1.5000 11.19
79 1.6352172 16.1829209 1.262533 0.108044 10.85484 1.5000 9.93
79 5.2103675 14.2777820 0.777049 0.179898 35.13864 1.5000 9.90
80 -14.5875156 10.8515971 0.586897 0.165819 124.75138 1.5000 6.65
80 -7.6472594 -3.8905882 0.650377 0.190278 138.08826 1.5000 6.39
80 6.9882215 -13.5565861 0.532227 0.226926 129.66115 1.5000 6.42
81 2.6039489 3.6745080 2.930790 0.046937 31.86086 1.5000 11.45
81 10.7289822 -11.1572676 1.939554 0.057465 58.34173 1.5000 11.68
-81 -10.7738999 12.8109455 1.862324 0.067971 46.64140 1.5000 11.54
+81 -10.7738999 12.8109455 1.862324 0.067971 46.64140 1.5000 11.54
\ No newline at end of file
diff --git a/Benchmarks/ConfigFiles/m1931_modif_217Pot.par b/Benchmarks/ConfigFiles/m1931_modif_217Pot.par
index bedcbb3..0121a70 100644
--- a/Benchmarks/ConfigFiles/m1931_modif_217Pot.par
+++ b/Benchmarks/ConfigFiles/m1931_modif_217Pot.par
@@ -1,3061 +1,3062 @@
#The user defines which mode of lenstool he wants to run here
runmode
reference 0 0.0 0.0
image 1 ../ConfigFiles/m1931_modif_210Img.txt
- nbgridcells 2000
- chi2 1
+ nbgridcells 3000
+ chi2 1 #Bug lenstool possible si image multiple n'est pas multiple!
+ verbose 0
end
image
multfile 1 ../ConfigFiles/m1931_modif_210Img.txt
forme -1
end
grille
nombre 128
polaire 0
nlentille 217
nlens_opt 0
end
cosmologie
model 1
H0 70.000
omegaM 0.300
omegaX 0.700
omegaK 0.
wX -1.000
wa 0.000
end
champ
dmax 150.000
end
potentiel O1
profil 81
x_centre -3.947
y_centre -0.275
ellipticite 0.594
angle_pos 27.637
core_radius 22.158
core_radius_kpc 141.291
cut_radius 156.824
cut_radius_kpc 1000.000
v_disp 1011.884
z_lens 0.5440
end
potentiel O2
profil 81
x_centre -28.172
y_centre -36.312
ellipticite 0.443
angle_pos 34.760
core_radius 4.069
core_radius_kpc 25.949
cut_radius 156.824
cut_radius_kpc 1000.000
v_disp 383.278
z_lens 0.5440
end
potentiel O3
profil 81
x_centre 21.380
y_centre 43.405
ellipticite 0.168
angle_pos 51.971
core_radius 11.652
core_radius_kpc 74.299
cut_radius 156.824
cut_radius_kpc 1000.000
v_disp 622.549
z_lens 0.5440
end
potentiel O4
profil 81
x_centre -17.407
y_centre 102.179
ellipticite 0.506
angle_pos 129.205
core_radius 0.821
core_radius_kpc 5.236
cut_radius 156.824
cut_radius_kpc 1000.000
v_disp 381.988
z_lens 0.5440
end
potentiel O5
profil 81
x_centre -41.433
y_centre -49.160
ellipticite 0.462
angle_pos 175.686
core_radius 15.145
core_radius_kpc 96.573
cut_radius 156.824
cut_radius_kpc 1000.000
v_disp 287.477
z_lens 0.5440
end
potentiel O6
profil 81
x_centre 0.000
y_centre 0.000
ellipticite 0.200
angle_pos 34.000
core_radius 2.022
core_radius_kpc 12.896
cut_radius 23.940
cut_radius_kpc 152.653
v_disp 396.956
z_lens 0.5440
end
potentiel O7
profil 81
x_centre 3.160
y_centre -11.102
ellipticite 0.429
angle_pos 134.000
core_radius 0.022
core_radius_kpc 0.142
cut_radius 6.741
cut_radius_kpc 42.985
v_disp 151.585
z_lens 0.5440
end
potentiel 1824
profil 81
x_centre -17.171
y_centre 101.181
ellipticite 0.319
angle_pos -58.200
core_radius 0.051
core_radius_kpc 0.326
cut_radius 19.414
cut_radius_kpc 123.793
v_disp 211.648
mag 18.960
z_lens 0.5440
end
potentiel 521
profil 81
x_centre -25.711
y_centre -32.225
ellipticite 0.095
angle_pos 43.500
core_radius 0.044
core_radius_kpc 0.279
cut_radius 16.600
cut_radius_kpc 105.851
v_disp 195.710
mag 19.300
z_lens 0.5440
end
potentiel 363
profil 81
x_centre -44.856
y_centre -54.794
ellipticite 0.061
angle_pos 52.300
core_radius 0.036
core_radius_kpc 0.231
cut_radius 13.744
cut_radius_kpc 87.639
v_disp 178.079
mag 19.710
z_lens 0.5440
end
potentiel 263
profil 81
x_centre 67.459
y_centre -60.867
ellipticite 0.811
angle_pos 81.400
core_radius 0.034
core_radius_kpc 0.216
cut_radius 12.886
cut_radius_kpc 82.167
v_disp 172.430
mag 19.850
z_lens 0.5440
end
potentiel 706
profil 81
x_centre 30.887
y_centre -16.616
ellipticite 0.301
angle_pos -29.600
core_radius 0.032
core_radius_kpc 0.203
cut_radius 12.081
cut_radius_kpc 77.036
v_disp 166.960
mag 19.990
z_lens 0.5440
end
potentiel 1046
profil 81
x_centre -29.238
y_centre 16.291
ellipticite 0.699
angle_pos 76.600
core_radius 0.032
core_radius_kpc 0.201
cut_radius 11.970
cut_radius_kpc 76.330
v_disp 166.193
mag 20.010
z_lens 0.5440
end
potentiel 1285
profil 81
x_centre -50.255
y_centre 35.707
ellipticite 0.069
angle_pos -66.500
core_radius 0.032
core_radius_kpc 0.201
cut_radius 11.970
cut_radius_kpc 76.330
v_disp 166.193
mag 20.010
z_lens 0.5440
end
potentiel 1205
profil 81
x_centre 112.259
y_centre 28.363
ellipticite 0.296
angle_pos 31.500
core_radius 0.031
core_radius_kpc 0.200
cut_radius 11.915
cut_radius_kpc 75.979
v_disp 165.811
mag 20.020
z_lens 0.5440
end
potentiel 1401
profil 81
x_centre 19.267
y_centre 48.038
ellipticite 0.232
angle_pos -77.300
core_radius 0.031
core_radius_kpc 0.197
cut_radius 11.752
cut_radius_kpc 74.937
v_disp 164.670
mag 20.050
z_lens 0.5440
end
potentiel 1413
profil 81
x_centre 16.332
y_centre 47.380
ellipticite 0.385
angle_pos -44.600
core_radius 0.031
core_radius_kpc 0.196
cut_radius 11.698
cut_radius_kpc 74.593
v_disp 164.291
mag 20.060
z_lens 0.5440
end
potentiel 953
profil 81
x_centre -62.420
y_centre 4.573
ellipticite 0.074
angle_pos -6.600
core_radius 0.030
core_radius_kpc 0.193
cut_radius 11.484
cut_radius_kpc 73.231
v_disp 162.785
mag 20.100
z_lens 0.5440
end
potentiel 798
profil 81
x_centre -16.125
y_centre -7.725
ellipticite 0.600
angle_pos 77.200
core_radius 0.030
core_radius_kpc 0.191
cut_radius 11.379
cut_radius_kpc 72.560
v_disp 162.037
mag 20.120
z_lens 0.5440
end
potentiel 864
profil 81
x_centre 56.711
y_centre -5.554
ellipticite 0.699
angle_pos 12.100
core_radius 0.029
core_radius_kpc 0.186
cut_radius 11.069
cut_radius_kpc 70.582
v_disp 159.813
mag 20.180
z_lens 0.5440
end
potentiel 904
profil 81
x_centre 4.019
y_centre 3.592
ellipticite 0.501
angle_pos 59.900
core_radius 0.029
core_radius_kpc 0.182
cut_radius 10.867
cut_radius_kpc 69.294
v_disp 158.348
mag 20.220
z_lens 0.5440
end
potentiel 315
profil 81
x_centre -1.846
y_centre -61.590
ellipticite 0.192
angle_pos 49.600
core_radius 0.028
core_radius_kpc 0.182
cut_radius 10.817
cut_radius_kpc 68.976
v_disp 157.984
mag 20.230
z_lens 0.5440
end
potentiel 1390
profil 81
x_centre 126.886
y_centre 55.159
ellipticite 0.819
angle_pos -67.600
core_radius 0.028
core_radius_kpc 0.177
cut_radius 10.522
cut_radius_kpc 67.096
v_disp 155.816
mag 20.290
z_lens 0.5440
end
potentiel 1776
profil 81
x_centre -1.086
y_centre 86.173
ellipticite 0.899
angle_pos 28.700
core_radius 0.027
core_radius_kpc 0.175
cut_radius 10.426
cut_radius_kpc 66.481
v_disp 155.101
mag 20.310
z_lens 0.5440
end
potentiel 1377
profil 81
x_centre 90.400
y_centre 42.432
ellipticite 0.257
angle_pos -85.900
core_radius 0.027
core_radius_kpc 0.173
cut_radius 10.330
cut_radius_kpc 65.871
v_disp 154.388
mag 20.330
z_lens 0.5440
end
potentiel 806
profil 81
x_centre -3.619
y_centre -4.564
ellipticite 0.498
angle_pos 62.600
core_radius 0.027
core_radius_kpc 0.171
cut_radius 10.188
cut_radius_kpc 64.967
v_disp 153.325
mag 20.360
z_lens 0.5440
end
potentiel 1010
profil 81
x_centre -84.443
y_centre 10.517
ellipticite 0.205
angle_pos 53.900
core_radius 0.027
core_radius_kpc 0.169
cut_radius 10.095
cut_radius_kpc 64.372
v_disp 152.621
mag 20.380
z_lens 0.5440
end
potentiel 1309
profil 81
x_centre 35.328
y_centre 35.538
ellipticite 0.415
angle_pos -29.600
core_radius 0.027
core_radius_kpc 0.169
cut_radius 10.095
cut_radius_kpc 64.372
v_disp 152.621
mag 20.380
z_lens 0.5440
end
potentiel 569
profil 81
x_centre -73.383
y_centre -30.375
ellipticite 0.237
angle_pos -47.200
core_radius 0.027
core_radius_kpc 0.169
cut_radius 10.095
cut_radius_kpc 64.372
v_disp 152.621
mag 20.380
z_lens 0.5440
end
potentiel 569
profil 81
x_centre -73.383
y_centre -30.375
ellipticite 0.237
angle_pos -47.200
core_radius 0.027
core_radius_kpc 0.169
cut_radius 10.095
cut_radius_kpc 64.372
v_disp 152.621
mag 20.380
z_lens 0.5440
end
potentiel 555
profil 81
x_centre -66.266
y_centre -32.185
ellipticite 0.724
angle_pos 51.500
core_radius 0.026
core_radius_kpc 0.168
cut_radius 10.002
cut_radius_kpc 63.782
v_disp 151.919
mag 20.400
z_lens 0.5440
end
potentiel 703
profil 81
x_centre 48.615
y_centre -16.848
ellipticite 0.095
angle_pos 59.000
core_radius 0.026
core_radius_kpc 0.166
cut_radius 9.911
cut_radius_kpc 63.197
v_disp 151.221
mag 20.420
z_lens 0.5440
end
potentiel 1009
profil 81
x_centre -61.953
y_centre 9.350
ellipticite 0.819
angle_pos 34.500
core_radius 0.026
core_radius_kpc 0.166
cut_radius 9.865
cut_radius_kpc 62.906
v_disp 150.874
mag 20.430
z_lens 0.5440
end
potentiel 1631
profil 81
x_centre 6.295
y_centre 68.144
ellipticite 0.153
angle_pos 67.000
core_radius 0.026
core_radius_kpc 0.163
cut_radius 9.730
cut_radius_kpc 62.043
v_disp 149.835
mag 20.460
z_lens 0.5440
end
potentiel 659
profil 81
x_centre 20.192
y_centre -20.872
ellipticite 0.385
angle_pos -50.700
core_radius 0.026
core_radius_kpc 0.163
cut_radius 9.730
cut_radius_kpc 62.043
v_disp 149.835
mag 20.460
z_lens 0.5440
end
potentiel 1552
profil 81
x_centre -3.558
y_centre 58.940
ellipticite 0.816
angle_pos 34.900
core_radius 0.025
core_radius_kpc 0.162
cut_radius 9.641
cut_radius_kpc 61.475
v_disp 149.147
mag 20.480
z_lens 0.5440
end
potentiel 661
profil 81
x_centre 30.071
y_centre -20.960
ellipticite 0.166
angle_pos 44.200
core_radius 0.025
core_radius_kpc 0.162
cut_radius 9.641
cut_radius_kpc 61.475
v_disp 149.147
mag 20.480
z_lens 0.5440
end
potentiel 339
profil 81
x_centre -43.741
y_centre -57.996
ellipticite 0.800
angle_pos 66.000
core_radius 0.025
core_radius_kpc 0.160
cut_radius 9.508
cut_radius_kpc 60.631
v_disp 148.120
mag 20.510
z_lens 0.5440
end
potentiel 1012
profil 81
x_centre 2.586
y_centre 9.067
ellipticite 0.105
angle_pos -52.100
core_radius 0.025
core_radius_kpc 0.157
cut_radius 9.378
cut_radius_kpc 59.799
v_disp 147.100
mag 20.540
z_lens 0.5440
end
potentiel 1585
profil 81
x_centre -2.326
y_centre 60.886
ellipticite 0.257
angle_pos 29.500
core_radius 0.024
core_radius_kpc 0.154
cut_radius 9.164
cut_radius_kpc 58.438
v_disp 145.416
mag 20.590
z_lens 0.5440
end
potentiel 465
profil 81
x_centre -39.868
y_centre -45.479
ellipticite 0.087
angle_pos 21.200
core_radius 0.024
core_radius_kpc 0.152
cut_radius 9.080
cut_radius_kpc 57.902
v_disp 144.748
mag 20.610
z_lens 0.5440
end
potentiel 500
profil 81
x_centre 46.797
y_centre -41.394
ellipticite 0.415
angle_pos 14.400
core_radius 0.024
core_radius_kpc 0.152
cut_radius 9.080
cut_radius_kpc 57.902
v_disp 144.748
mag 20.610
z_lens 0.5440
end
potentiel 888
profil 81
x_centre -7.647
y_centre -2.333
ellipticite 0.105
angle_pos 22.300
core_radius 0.024
core_radius_kpc 0.151
cut_radius 8.997
cut_radius_kpc 57.371
v_disp 144.083
mag 20.630
z_lens 0.5440
end
potentiel 880
profil 81
x_centre 19.815
y_centre -1.565
ellipticite 0.550
angle_pos -71.800
core_radius 0.023
core_radius_kpc 0.149
cut_radius 8.874
cut_radius_kpc 56.584
v_disp 143.091
mag 20.660
z_lens 0.5440
end
potentiel 1181
profil 81
x_centre 25.415
y_centre 22.961
ellipticite 0.198
angle_pos 20.700
core_radius 0.023
core_radius_kpc 0.146
cut_radius 8.712
cut_radius_kpc 55.551
v_disp 141.779
mag 20.700
z_lens 0.5440
end
potentiel 829
profil 81
x_centre 147.582
y_centre -8.551
ellipticite 0.246
angle_pos -20.900
core_radius 0.023
core_radius_kpc 0.146
cut_radius 8.672
cut_radius_kpc 55.296
v_disp 141.453
mag 20.710
z_lens 0.5440
end
potentiel 1073
profil 81
x_centre -95.435
y_centre 14.332
ellipticite 0.105
angle_pos 5.600
core_radius 0.022
core_radius_kpc 0.143
cut_radius 8.513
cut_radius_kpc 54.287
v_disp 140.157
mag 20.750
z_lens 0.5440
end
potentiel 1531
profil 81
x_centre 7.113
y_centre 53.268
ellipticite 0.471
angle_pos -16.900
core_radius 0.022
core_radius_kpc 0.143
cut_radius 8.513
cut_radius_kpc 54.287
v_disp 140.157
mag 20.750
z_lens 0.5440
end
potentiel 987
profil 81
x_centre -8.215
y_centre 6.494
ellipticite 0.385
angle_pos 73.500
core_radius 0.022
core_radius_kpc 0.140
cut_radius 8.320
cut_radius_kpc 53.051
v_disp 138.552
mag 20.800
z_lens 0.5440
end
potentiel 595
profil 81
x_centre 23.467
y_centre -28.479
ellipticite 0.280
angle_pos -68.100
core_radius 0.022
core_radius_kpc 0.138
cut_radius 8.206
cut_radius_kpc 52.323
v_disp 137.598
mag 20.830
z_lens 0.5440
end
potentiel 884
profil 81
x_centre -105.029
y_centre -4.239
ellipticite 0.492
angle_pos -5.400
core_radius 0.021
core_radius_kpc 0.136
cut_radius 8.130
cut_radius_kpc 51.844
v_disp 136.966
mag 20.850
z_lens 0.5440
end
potentiel 1756
profil 81
x_centre 78.520
y_centre 79.622
ellipticite 0.424
angle_pos 12.100
core_radius 0.021
core_radius_kpc 0.135
cut_radius 8.056
cut_radius_kpc 51.368
v_disp 136.337
mag 20.870
z_lens 0.5440
end
potentiel 672
profil 81
x_centre 5.966
y_centre -19.973
ellipticite 0.451
angle_pos 83.700
core_radius 0.021
core_radius_kpc 0.135
cut_radius 8.056
cut_radius_kpc 51.368
v_disp 136.337
mag 20.870
z_lens 0.5440
end
potentiel 1170
profil 81
x_centre -16.595
y_centre 21.765
ellipticite 0.220
angle_pos -57.600
core_radius 0.021
core_radius_kpc 0.134
cut_radius 7.982
cut_radius_kpc 50.897
v_disp 135.710
mag 20.890
z_lens 0.5440
end
potentiel 895
profil 81
x_centre -22.333
y_centre -2.388
ellipticite 0.117
angle_pos -39.900
core_radius 0.021
core_radius_kpc 0.134
cut_radius 7.982
cut_radius_kpc 50.897
v_disp 135.710
mag 20.890
z_lens 0.5440
end
potentiel 414
profil 81
x_centre -38.966
y_centre -52.176
ellipticite 0.198
angle_pos 15.500
core_radius 0.021
core_radius_kpc 0.133
cut_radius 7.945
cut_radius_kpc 50.663
v_disp 135.398
mag 20.900
z_lens 0.5440
end
potentiel 1288
profil 81
x_centre 1.920
y_centre 31.968
ellipticite 0.342
angle_pos 9.200
core_radius 0.021
core_radius_kpc 0.131
cut_radius 7.800
cut_radius_kpc 49.739
v_disp 134.157
mag 20.940
z_lens 0.5440
end
potentiel 763
profil 81
x_centre 37.630
y_centre -13.753
ellipticite 0.600
angle_pos -43.800
core_radius 0.021
core_radius_kpc 0.131
cut_radius 7.800
cut_radius_kpc 49.739
v_disp 134.157
mag 20.940
z_lens 0.5440
end
potentiel 1412
profil 81
x_centre 49.349
y_centre 42.314
ellipticite 0.000
angle_pos 53.500
core_radius 0.020
core_radius_kpc 0.130
cut_radius 7.764
cut_radius_kpc 49.510
v_disp 133.848
mag 20.950
z_lens 0.5440
end
potentiel 1062
profil 81
x_centre 114.777
y_centre 10.203
ellipticite 0.704
angle_pos 8.100
core_radius 0.020
core_radius_kpc 0.130
cut_radius 7.729
cut_radius_kpc 49.283
v_disp 133.541
mag 20.960
z_lens 0.5440
end
potentiel 1574
profil 81
x_centre -63.476
y_centre 59.932
ellipticite 0.600
angle_pos -65.100
core_radius 0.020
core_radius_kpc 0.127
cut_radius 7.588
cut_radius_kpc 48.383
v_disp 132.316
mag 21.000
z_lens 0.5440
end
potentiel 1994
profil 81
x_centre 23.977
y_centre 118.125
ellipticite 0.528
angle_pos 52.100
core_radius 0.020
core_radius_kpc 0.127
cut_radius 7.588
cut_radius_kpc 48.383
v_disp 132.316
mag 21.000
z_lens 0.5440
end
potentiel 92
profil 81
x_centre 24.242
y_centre -98.757
ellipticite 0.342
angle_pos -37.200
core_radius 0.020
core_radius_kpc 0.127
cut_radius 7.588
cut_radius_kpc 48.383
v_disp 132.316
mag 21.000
z_lens 0.5440
end
potentiel 580
profil 81
x_centre -37.932
y_centre -32.860
ellipticite 0.658
angle_pos -24.900
core_radius 0.019
core_radius_kpc 0.124
cut_radius 7.381
cut_radius_kpc 47.065
v_disp 130.501
mag 21.060
z_lens 0.5440
end
potentiel 1059
profil 81
x_centre 25.824
y_centre 11.355
ellipticite 0.133
angle_pos 18.500
core_radius 0.018
core_radius_kpc 0.118
cut_radius 7.016
cut_radius_kpc 44.740
v_disp 127.237
mag 21.170
z_lens 0.5440
end
potentiel 1001
profil 81
x_centre 46.608
y_centre 7.062
ellipticite 0.000
angle_pos -84.700
core_radius 0.018
core_radius_kpc 0.117
cut_radius 6.984
cut_radius_kpc 44.534
v_disp 126.944
mag 21.180
z_lens 0.5440
end
potentiel 1440
profil 81
x_centre 72.603
y_centre 47.167
ellipticite 0.658
angle_pos -50.000
core_radius 0.018
core_radius_kpc 0.116
cut_radius 6.888
cut_radius_kpc 43.923
v_disp 126.070
mag 21.210
z_lens 0.5440
end
potentiel 495
profil 81
x_centre -17.713
y_centre -43.740
ellipticite 0.600
angle_pos 22.000
core_radius 0.018
core_radius_kpc 0.116
cut_radius 6.888
cut_radius_kpc 43.923
v_disp 126.070
mag 21.210
z_lens 0.5440
end
potentiel 558
profil 81
x_centre -16.515
y_centre -33.903
ellipticite 0.220
angle_pos -50.300
core_radius 0.018
core_radius_kpc 0.116
cut_radius 6.888
cut_radius_kpc 43.923
v_disp 126.070
mag 21.210
z_lens 0.5440
end
potentiel 1546
profil 81
x_centre 35.359
y_centre 54.831
ellipticite 0.280
angle_pos -47.000
core_radius 0.018
core_radius_kpc 0.114
cut_radius 6.794
cut_radius_kpc 43.321
v_disp 125.203
mag 21.240
z_lens 0.5440
end
potentiel 1249
profil 81
x_centre 30.816
y_centre 28.247
ellipticite 0.246
angle_pos -48.800
core_radius 0.018
core_radius_kpc 0.113
cut_radius 6.731
cut_radius_kpc 42.924
v_disp 124.627
mag 21.260
z_lens 0.5440
end
potentiel 461
profil 81
x_centre -63.940
y_centre -46.818
ellipticite 0.324
angle_pos 6.600
core_radius 0.018
core_radius_kpc 0.113
cut_radius 6.731
cut_radius_kpc 42.924
v_disp 124.627
mag 21.260
z_lens 0.5440
end
potentiel 817
profil 81
x_centre -12.932
y_centre -8.450
ellipticite 0.198
angle_pos -31.200
core_radius 0.018
core_radius_kpc 0.113
cut_radius 6.731
cut_radius_kpc 42.924
v_disp 124.627
mag 21.260
z_lens 0.5440
end
potentiel 1529
profil 81
x_centre -87.434
y_centre 52.897
ellipticite 0.867
angle_pos -29.100
core_radius 0.017
core_radius_kpc 0.111
cut_radius 6.609
cut_radius_kpc 42.140
v_disp 123.485
mag 21.300
z_lens 0.5440
end
potentiel 1458
profil 81
x_centre 18.628
y_centre 46.026
ellipticite 0.117
angle_pos -74.700
core_radius 0.017
core_radius_kpc 0.109
cut_radius 6.518
cut_radius_kpc 41.562
v_disp 122.635
mag 21.330
z_lens 0.5440
end
potentiel 749
profil 81
x_centre -52.959
y_centre -15.151
ellipticite 0.133
angle_pos -67.900
core_radius 0.017
core_radius_kpc 0.109
cut_radius 6.518
cut_radius_kpc 41.562
v_disp 122.635
mag 21.330
z_lens 0.5440
end
potentiel 402
profil 81
x_centre -51.055
y_centre -53.688
ellipticite 0.117
angle_pos -79.300
core_radius 0.017
core_radius_kpc 0.108
cut_radius 6.458
cut_radius_kpc 41.181
v_disp 122.071
mag 21.350
z_lens 0.5440
end
potentiel 849
profil 81
x_centre 19.478
y_centre -5.130
ellipticite 0.471
angle_pos 31.000
core_radius 0.017
core_radius_kpc 0.108
cut_radius 6.458
cut_radius_kpc 41.181
v_disp 122.071
mag 21.350
z_lens 0.5440
end
potentiel 715
profil 81
x_centre 94.758
y_centre -19.879
ellipticite 0.385
angle_pos 6.700
core_radius 0.017
core_radius_kpc 0.108
cut_radius 6.428
cut_radius_kpc 40.992
v_disp 121.790
mag 21.360
z_lens 0.5440
end
potentiel 2145
profil 81
x_centre 26.418
y_centre 152.034
ellipticite 0.424
angle_pos -61.300
core_radius 0.017
core_radius_kpc 0.106
cut_radius 6.311
cut_radius_kpc 40.243
v_disp 120.674
mag 21.400
z_lens 0.5440
end
potentiel 2145
profil 81
x_centre 26.418
y_centre 152.034
ellipticite 0.424
angle_pos -61.300
core_radius 0.017
core_radius_kpc 0.106
cut_radius 6.311
cut_radius_kpc 40.243
v_disp 120.674
mag 21.400
z_lens 0.5440
end
potentiel 70
profil 81
x_centre 25.255
y_centre -102.010
ellipticite 0.342
angle_pos -47.800
core_radius 0.017
core_radius_kpc 0.106
cut_radius 6.311
cut_radius_kpc 40.243
v_disp 120.674
mag 21.400
z_lens 0.5440
end
potentiel 211
profil 81
x_centre 41.229
y_centre -78.527
ellipticite 0.658
angle_pos -3.200
core_radius 0.017
core_radius_kpc 0.105
cut_radius 6.282
cut_radius_kpc 40.059
v_disp 120.396
mag 21.410
z_lens 0.5440
end
potentiel 1688
profil 81
x_centre -0.894
y_centre 71.496
ellipticite 0.133
angle_pos 77.900
core_radius 0.016
core_radius_kpc 0.105
cut_radius 6.253
cut_radius_kpc 39.874
v_disp 120.119
mag 21.420
z_lens 0.5440
end
potentiel 937
profil 81
x_centre 37.904
y_centre 1.546
ellipticite 0.424
angle_pos -73.300
core_radius 0.016
core_radius_kpc 0.104
cut_radius 6.225
cut_radius_kpc 39.691
v_disp 119.843
mag 21.430
z_lens 0.5440
end
potentiel 2067
profil 81
x_centre 36.025
y_centre 138.420
ellipticite 0.087
angle_pos 35.700
core_radius 0.016
core_radius_kpc 0.104
cut_radius 6.196
cut_radius_kpc 39.509
v_disp 119.568
mag 21.440
z_lens 0.5440
end
potentiel 836
profil 81
x_centre 139.232
y_centre -6.932
ellipticite 0.670
angle_pos -35.800
core_radius 0.016
core_radius_kpc 0.104
cut_radius 6.196
cut_radius_kpc 39.509
v_disp 119.568
mag 21.440
z_lens 0.5440
end
potentiel 1524
profil 81
x_centre -37.456
y_centre 53.916
ellipticite 0.471
angle_pos 79.200
core_radius 0.016
core_radius_kpc 0.103
cut_radius 6.139
cut_radius_kpc 39.147
v_disp 119.018
mag 21.460
z_lens 0.5440
end
potentiel 1264
profil 81
x_centre -31.987
y_centre 33.080
ellipticite 0.679
angle_pos -82.100
core_radius 0.016
core_radius_kpc 0.103
cut_radius 6.111
cut_radius_kpc 38.967
v_disp 118.744
mag 21.470
z_lens 0.5440
end
potentiel 669
profil 81
x_centre -16.413
y_centre -23.715
ellipticite 0.324
angle_pos 20.000
core_radius 0.016
core_radius_kpc 0.103
cut_radius 6.111
cut_radius_kpc 38.967
v_disp 118.744
mag 21.470
z_lens 0.5440
end
potentiel 1970
profil 81
x_centre -20.691
y_centre 111.618
ellipticite 0.385
angle_pos 1.600
core_radius 0.016
core_radius_kpc 0.102
cut_radius 6.083
cut_radius_kpc 38.788
v_disp 118.471
mag 21.480
z_lens 0.5440
end
potentiel 1207
profil 81
x_centre 73.817
y_centre 24.093
ellipticite 0.342
angle_pos -81.100
core_radius 0.016
core_radius_kpc 0.101
cut_radius 6.027
cut_radius_kpc 38.432
v_disp 117.927
mag 21.500
z_lens 0.5440
end
potentiel 1969
profil 81
x_centre 6.830
y_centre 112.550
ellipticite 0.385
angle_pos -68.300
core_radius 0.016
core_radius_kpc 0.101
cut_radius 5.999
cut_radius_kpc 38.256
v_disp 117.656
mag 21.510
z_lens 0.5440
end
potentiel 1363
profil 81
x_centre 15.580
y_centre 35.661
ellipticite 0.528
angle_pos 13.600
core_radius 0.015
core_radius_kpc 0.098
cut_radius 5.863
cut_radius_kpc 37.385
v_disp 116.309
mag 21.560
z_lens 0.5440
end
potentiel 1573
profil 81
x_centre 13.329
y_centre 57.151
ellipticite 0.342
angle_pos 3.300
core_radius 0.015
core_radius_kpc 0.098
cut_radius 5.863
cut_radius_kpc 37.385
v_disp 116.309
mag 21.560
z_lens 0.5440
end
potentiel 448
profil 81
x_centre 49.083
y_centre -49.641
ellipticite 0.324
angle_pos 57.100
core_radius 0.015
core_radius_kpc 0.098
cut_radius 5.863
cut_radius_kpc 37.385
v_disp 116.309
mag 21.560
z_lens 0.5440
end
potentiel 1590
profil 81
x_centre 50.581
y_centre 59.075
ellipticite 0.153
angle_pos -47.300
core_radius 0.015
core_radius_kpc 0.098
cut_radius 5.836
cut_radius_kpc 37.213
v_disp 116.041
mag 21.570
z_lens 0.5440
end
potentiel 683
profil 81
x_centre -58.310
y_centre -22.309
ellipticite 0.324
angle_pos 68.200
core_radius 0.015
core_radius_kpc 0.098
cut_radius 5.809
cut_radius_kpc 37.042
v_disp 115.775
mag 21.580
z_lens 0.5440
end
potentiel 1583
profil 81
x_centre 11.630
y_centre 59.053
ellipticite 0.133
angle_pos 45.700
core_radius 0.015
core_radius_kpc 0.097
cut_radius 5.782
cut_radius_kpc 36.872
v_disp 115.508
mag 21.590
z_lens 0.5440
end
potentiel 1790
profil 81
x_centre 99.193
y_centre 85.405
ellipticite 0.550
angle_pos -82.600
core_radius 0.015
core_radius_kpc 0.096
cut_radius 5.729
cut_radius_kpc 36.534
v_disp 114.978
mag 21.610
z_lens 0.5440
end
potentiel 508
profil 81
x_centre -26.682
y_centre -39.870
ellipticite 0.541
angle_pos -71.300
core_radius 0.015
core_radius_kpc 0.095
cut_radius 5.677
cut_radius_kpc 36.199
v_disp 114.449
mag 21.630
z_lens 0.5440
end
potentiel 1191
profil 81
x_centre -14.413
y_centre 21.081
ellipticite 0.438
angle_pos -7.600
core_radius 0.015
core_radius_kpc 0.094
cut_radius 5.573
cut_radius_kpc 35.538
v_disp 113.400
mag 21.670
z_lens 0.5440
end
potentiel 335
profil 81
x_centre 27.945
y_centre -63.954
ellipticite 0.153
angle_pos -0.700
core_radius 0.014
core_radius_kpc 0.092
cut_radius 5.472
cut_radius_kpc 34.890
v_disp 112.360
mag 21.710
z_lens 0.5440
end
potentiel 1423
profil 81
x_centre 66.962
y_centre 41.650
ellipticite 0.280
angle_pos -57.900
core_radius 0.014
core_radius_kpc 0.091
cut_radius 5.396
cut_radius_kpc 34.411
v_disp 111.587
mag 21.740
z_lens 0.5440
end
potentiel 374
profil 81
x_centre -19.978
y_centre -58.703
ellipticite 0.600
angle_pos -86.500
core_radius 0.014
core_radius_kpc 0.090
cut_radius 5.372
cut_radius_kpc 34.253
v_disp 111.330
mag 21.750
z_lens 0.5440
end
potentiel 902
profil 81
x_centre 141.528
y_centre -2.009
ellipticite 0.724
angle_pos 47.400
core_radius 0.014
core_radius_kpc 0.089
cut_radius 5.298
cut_radius_kpc 33.783
v_disp 110.564
mag 21.780
z_lens 0.5440
end
potentiel 732
profil 81
x_centre 9.619
y_centre -18.079
ellipticite 0.153
angle_pos -35.700
core_radius 0.014
core_radius_kpc 0.087
cut_radius 5.177
cut_radius_kpc 33.014
v_disp 109.298
mag 21.830
z_lens 0.5440
end
potentiel 827
profil 81
x_centre -27.188
y_centre -9.758
ellipticite 0.000
angle_pos 26.000
core_radius 0.013
core_radius_kpc 0.086
cut_radius 5.106
cut_radius_kpc 32.561
v_disp 108.546
mag 21.860
z_lens 0.5440
end
potentiel 1065
profil 81
x_centre -0.437
y_centre 11.831
ellipticite 0.600
angle_pos -84.600
core_radius 0.013
core_radius_kpc 0.085
cut_radius 5.083
cut_radius_kpc 32.411
v_disp 108.296
mag 21.870
z_lens 0.5440
end
potentiel 1700
profil 81
x_centre -8.444
y_centre 73.209
ellipticite 0.246
angle_pos -77.800
core_radius 0.013
core_radius_kpc 0.085
cut_radius 5.083
cut_radius_kpc 32.411
v_disp 108.296
mag 21.870
z_lens 0.5440
end
potentiel 174
profil 81
x_centre -22.535
y_centre -83.316
ellipticite 0.280
angle_pos 72.600
core_radius 0.013
core_radius_kpc 0.085
cut_radius 5.083
cut_radius_kpc 32.411
v_disp 108.296
mag 21.870
z_lens 0.5440
end
potentiel 1481
profil 81
x_centre 35.590
y_centre 47.814
ellipticite 0.133
angle_pos 81.900
core_radius 0.013
core_radius_kpc 0.085
cut_radius 5.060
cut_radius_kpc 32.262
v_disp 108.047
mag 21.880
z_lens 0.5440
end
potentiel 1892
profil 81
x_centre 43.731
y_centre 98.992
ellipticite 0.438
angle_pos 26.000
core_radius 0.013
core_radius_kpc 0.084
cut_radius 5.013
cut_radius_kpc 31.967
v_disp 107.551
mag 21.900
z_lens 0.5440
end
potentiel 789
profil 81
x_centre 20.110
y_centre -15.000
ellipticite 0.324
angle_pos -5.400
core_radius 0.013
core_radius_kpc 0.084
cut_radius 5.013
cut_radius_kpc 31.967
v_disp 107.551
mag 21.900
z_lens 0.5440
end
potentiel 1691
profil 81
x_centre 79.279
y_centre 70.447
ellipticite 0.000
angle_pos 47.200
core_radius 0.013
core_radius_kpc 0.084
cut_radius 4.990
cut_radius_kpc 31.820
v_disp 107.303
mag 21.910
z_lens 0.5440
end
potentiel 1502
profil 81
x_centre 15.052
y_centre 50.166
ellipticite 0.438
angle_pos 84.600
core_radius 0.013
core_radius_kpc 0.083
cut_radius 4.944
cut_radius_kpc 31.528
v_disp 106.810
mag 21.930
z_lens 0.5440
end
potentiel 1608
profil 81
x_centre -34.474
y_centre 61.106
ellipticite 0.133
angle_pos -26.800
core_radius 0.013
core_radius_kpc 0.083
cut_radius 4.944
cut_radius_kpc 31.528
v_disp 106.810
mag 21.930
z_lens 0.5440
end
potentiel 739
profil 81
x_centre 36.713
y_centre -17.194
ellipticite 0.280
angle_pos -68.000
core_radius 0.013
core_radius_kpc 0.083
cut_radius 4.944
cut_radius_kpc 31.528
v_disp 106.810
mag 21.930
z_lens 0.5440
end
potentiel 2021
profil 81
x_centre 56.653
y_centre 123.183
ellipticite 0.246
angle_pos 42.300
core_radius 0.013
core_radius_kpc 0.082
cut_radius 4.876
cut_radius_kpc 31.095
v_disp 106.075
mag 21.960
z_lens 0.5440
end
potentiel 2032
profil 81
x_centre -12.053
y_centre 126.509
ellipticite 0.180
angle_pos 59.700
core_radius 0.013
core_radius_kpc 0.082
cut_radius 4.876
cut_radius_kpc 31.095
v_disp 106.075
mag 21.960
z_lens 0.5440
end
potentiel 2032
profil 81
x_centre -12.053
y_centre 126.509
ellipticite 0.180
angle_pos 59.700
core_radius 0.013
core_radius_kpc 0.082
cut_radius 4.876
cut_radius_kpc 31.095
v_disp 106.075
mag 21.960
z_lens 0.5440
end
potentiel 1160
profil 81
x_centre 28.569
y_centre 19.335
ellipticite 0.508
angle_pos 84.200
core_radius 0.012
core_radius_kpc 0.078
cut_radius 4.657
cut_radius_kpc 29.696
v_disp 103.661
mag 22.060
z_lens 0.5440
end
potentiel 1480
profil 81
x_centre 5.642
y_centre 46.905
ellipticite 0.508
angle_pos 55.000
core_radius 0.012
core_radius_kpc 0.078
cut_radius 4.657
cut_radius_kpc 29.696
v_disp 103.661
mag 22.060
z_lens 0.5440
end
potentiel 1865
profil 81
x_centre 47.592
y_centre 96.187
ellipticite 0.153
angle_pos 81.000
core_radius 0.012
core_radius_kpc 0.078
cut_radius 4.657
cut_radius_kpc 29.696
v_disp 103.661
mag 22.060
z_lens 0.5440
end
potentiel 917
profil 81
x_centre -50.140
y_centre -2.214
ellipticite 0.438
angle_pos -36.900
core_radius 0.012
core_radius_kpc 0.077
cut_radius 4.614
cut_radius_kpc 29.424
v_disp 103.184
mag 22.080
z_lens 0.5440
end
potentiel 1569
profil 81
x_centre 32.821
y_centre 56.235
ellipticite 0.385
angle_pos 48.000
core_radius 0.012
core_radius_kpc 0.077
cut_radius 4.593
cut_radius_kpc 29.288
v_disp 102.947
mag 22.090
z_lens 0.5440
end
potentiel 816
profil 81
x_centre 36.345
y_centre -11.277
ellipticite 0.220
angle_pos 79.300
core_radius 0.012
core_radius_kpc 0.076
cut_radius 4.530
cut_radius_kpc 28.887
v_disp 102.238
mag 22.120
z_lens 0.5440
end
potentiel 714
profil 81
x_centre 85.269
y_centre -19.902
ellipticite 0.280
angle_pos -45.000
core_radius 0.012
core_radius_kpc 0.075
cut_radius 4.489
cut_radius_kpc 28.622
v_disp 101.768
mag 22.140
z_lens 0.5440
end
potentiel 1954
profil 81
x_centre -30.008
y_centre 108.414
ellipticite 0.000
angle_pos 47.500
core_radius 0.012
core_radius_kpc 0.075
cut_radius 4.468
cut_radius_kpc 28.490
v_disp 101.534
mag 22.150
z_lens 0.5440
end
potentiel 589
profil 81
x_centre 0.604
y_centre -31.459
ellipticite 0.753
angle_pos 69.000
core_radius 0.012
core_radius_kpc 0.074
cut_radius 4.386
cut_radius_kpc 27.970
v_disp 100.604
mag 22.190
z_lens 0.5440
end
potentiel 968
profil 81
x_centre 108.819
y_centre 1.631
ellipticite 0.690
angle_pos -41.000
core_radius 0.011
core_radius_kpc 0.073
cut_radius 4.326
cut_radius_kpc 27.586
v_disp 99.911
mag 22.220
z_lens 0.5440
end
potentiel 797
profil 81
x_centre 72.595
y_centre -14.067
ellipticite 0.180
angle_pos -4.100
core_radius 0.011
core_radius_kpc 0.072
cut_radius 4.287
cut_radius_kpc 27.333
v_disp 99.452
mag 22.240
z_lens 0.5440
end
potentiel 948
profil 81
x_centre 54.183
y_centre 0.479
ellipticite 0.438
angle_pos 66.400
core_radius 0.011
core_radius_kpc 0.072
cut_radius 4.267
cut_radius_kpc 27.208
v_disp 99.223
mag 22.250
z_lens 0.5440
end
potentiel 1441
profil 81
x_centre -12.633
y_centre 42.951
ellipticite 0.153
angle_pos 64.900
core_radius 0.011
core_radius_kpc 0.072
cut_radius 4.267
cut_radius_kpc 27.208
v_disp 99.223
mag 22.250
z_lens 0.5440
end
potentiel 235
profil 81
x_centre -13.923
y_centre -74.222
ellipticite 0.724
angle_pos -67.700
core_radius 0.011
core_radius_kpc 0.071
cut_radius 4.247
cut_radius_kpc 27.083
v_disp 98.995
mag 22.260
z_lens 0.5440
end
potentiel 1882
profil 81
x_centre -1.020
y_centre 96.964
ellipticite 0.438
angle_pos -5.500
core_radius 0.011
core_radius_kpc 0.071
cut_radius 4.228
cut_radius_kpc 26.958
v_disp 98.767
mag 22.270
z_lens 0.5440
end
potentiel 903
profil 81
x_centre -44.195
y_centre -3.324
ellipticite 0.385
angle_pos 83.400
core_radius 0.011
core_radius_kpc 0.071
cut_radius 4.208
cut_radius_kpc 26.835
v_disp 98.540
mag 22.280
z_lens 0.5440
end
potentiel 1549
profil 81
x_centre 40.294
y_centre 54.422
ellipticite 0.180
angle_pos 40.700
core_radius 0.011
core_radius_kpc 0.070
cut_radius 4.189
cut_radius_kpc 26.711
v_disp 98.314
mag 22.290
z_lens 0.5440
end
potentiel 1630
profil 81
x_centre 1.771
y_centre 63.893
ellipticite 0.508
angle_pos 78.200
core_radius 0.011
core_radius_kpc 0.070
cut_radius 4.189
cut_radius_kpc 26.711
v_disp 98.314
mag 22.290
z_lens 0.5440
end
potentiel 1611
profil 81
x_centre -66.449
y_centre 61.932
ellipticite 0.153
angle_pos -63.700
core_radius 0.011
core_radius_kpc 0.070
cut_radius 4.189
cut_radius_kpc 26.711
v_disp 98.314
mag 22.290
z_lens 0.5440
end
potentiel 2002
profil 81
x_centre 5.683
y_centre 116.983
ellipticite 0.438
angle_pos -5.700
core_radius 0.011
core_radius_kpc 0.069
cut_radius 4.132
cut_radius_kpc 26.345
v_disp 97.637
mag 22.320
z_lens 0.5440
end
potentiel 926
profil 81
x_centre 60.578
y_centre -2.288
ellipticite 0.000
angle_pos -62.800
core_radius 0.011
core_radius_kpc 0.069
cut_radius 4.094
cut_radius_kpc 26.103
v_disp 97.188
mag 22.340
z_lens 0.5440
end
potentiel 536
profil 81
x_centre 10.785
y_centre -39.088
ellipticite 0.600
angle_pos 31.700
core_radius 0.011
core_radius_kpc 0.068
cut_radius 4.056
cut_radius_kpc 25.864
v_disp 96.742
mag 22.360
z_lens 0.5440
end
potentiel 1049
profil 81
x_centre 11.516
y_centre 9.079
ellipticite 0.000
angle_pos 77.800
core_radius 0.011
core_radius_kpc 0.068
cut_radius 4.037
cut_radius_kpc 25.745
v_disp 96.519
mag 22.370
z_lens 0.5440
end
potentiel 1485
profil 81
x_centre 34.668
y_centre 47.406
ellipticite 0.000
angle_pos 73.800
core_radius 0.011
core_radius_kpc 0.068
cut_radius 4.037
cut_radius_kpc 25.745
v_disp 96.519
mag 22.370
z_lens 0.5440
end
potentiel 812
profil 81
x_centre -39.633
y_centre -12.253
ellipticite 0.220
angle_pos -50.200
core_radius 0.010
core_radius_kpc 0.065
cut_radius 3.856
cut_radius_kpc 24.586
v_disp 94.322
mag 22.470
z_lens 0.5440
end
potentiel 321
profil 81
x_centre -39.084
y_centre -65.841
ellipticite 0.508
angle_pos -54.900
core_radius 0.010
core_radius_kpc 0.063
cut_radius 3.768
cut_radius_kpc 24.027
v_disp 93.242
mag 22.520
z_lens 0.5440
end
potentiel 1478
profil 81
x_centre 90.819
y_centre 46.175
ellipticite 0.385
angle_pos 37.300
core_radius 0.010
core_radius_kpc 0.062
cut_radius 3.716
cut_radius_kpc 23.697
v_disp 92.600
mag 22.550
z_lens 0.5440
end
potentiel 1733
profil 81
x_centre 38.025
y_centre 75.604
ellipticite 0.600
angle_pos 72.800
core_radius 0.010
core_radius_kpc 0.061
cut_radius 3.648
cut_radius_kpc 23.265
v_disp 91.752
mag 22.590
z_lens 0.5440
end
potentiel 1024
profil 81
x_centre 4.888
y_centre 6.388
ellipticite 0.324
angle_pos -55.600
core_radius 0.010
core_radius_kpc 0.061
cut_radius 3.648
cut_radius_kpc 23.265
v_disp 91.752
mag 22.590
z_lens 0.5440
end
potentiel 1212
profil 81
x_centre 83.163
y_centre 22.050
ellipticite 0.000
angle_pos -6.900
core_radius 0.009
core_radius_kpc 0.060
cut_radius 3.549
cut_radius_kpc 22.631
v_disp 90.493
mag 22.650
z_lens 0.5440
end
potentiel 1735
profil 81
x_centre 35.189
y_centre 74.888
ellipticite 0.508
angle_pos -31.800
core_radius 0.009
core_radius_kpc 0.059
cut_radius 3.533
cut_radius_kpc 22.527
v_disp 90.285
mag 22.660
z_lens 0.5440
end
potentiel 844
profil 81
x_centre -43.565
y_centre -8.940
ellipticite 0.180
angle_pos 67.000
core_radius 0.009
core_radius_kpc 0.059
cut_radius 3.516
cut_radius_kpc 22.423
v_disp 90.077
mag 22.670
z_lens 0.5440
end
potentiel 1007
profil 81
x_centre -120.848
y_centre 6.127
ellipticite 0.324
angle_pos -86.700
core_radius 0.009
core_radius_kpc 0.059
cut_radius 3.500
cut_radius_kpc 22.320
v_disp 89.870
mag 22.680
z_lens 0.5440
end
potentiel 949
profil 81
x_centre -5.561
y_centre -1.157
ellipticite 0.000
angle_pos 49.000
core_radius 0.009
core_radius_kpc 0.058
cut_radius 3.468
cut_radius_kpc 22.115
v_disp 89.457
mag 22.700
z_lens 0.5440
end
potentiel 962
profil 81
x_centre 46.152
y_centre 1.135
ellipticite 0.385
angle_pos 69.800
core_radius 0.009
core_radius_kpc 0.057
cut_radius 3.405
cut_radius_kpc 21.712
v_disp 88.637
mag 22.740
z_lens 0.5440
end
potentiel 477
profil 81
x_centre -17.963
y_centre -47.532
ellipticite 0.000
angle_pos 49.600
core_radius 0.009
core_radius_kpc 0.056
cut_radius 3.358
cut_radius_kpc 21.414
v_disp 88.026
mag 22.770
z_lens 0.5440
end
potentiel 1553
profil 81
x_centre 41.109
y_centre 54.635
ellipticite 0.385
angle_pos -81.600
core_radius 0.009
core_radius_kpc 0.056
cut_radius 3.327
cut_radius_kpc 21.218
v_disp 87.622
mag 22.790
z_lens 0.5440
end
potentiel 1952
profil 81
x_centre 78.745
y_centre 107.303
ellipticite 0.220
angle_pos -80.500
core_radius 0.009
core_radius_kpc 0.055
cut_radius 3.297
cut_radius_kpc 21.023
v_disp 87.219
mag 22.810
z_lens 0.5440
end
potentiel 1279
profil 81
x_centre 49.006
y_centre 29.323
ellipticite 0.180
angle_pos -19.500
core_radius 0.009
core_radius_kpc 0.055
cut_radius 3.267
cut_radius_kpc 20.830
v_disp 86.819
mag 22.830
z_lens 0.5440
end
potentiel 1622
profil 81
x_centre 66.244
y_centre 62.535
ellipticite 0.220
angle_pos -63.300
core_radius 0.008
core_radius_kpc 0.053
cut_radius 3.163
cut_radius_kpc 20.170
v_disp 85.431
mag 22.900
z_lens 0.5440
end
potentiel 967
profil 81
x_centre -9.951
y_centre 0.986
ellipticite 0.280
angle_pos 49.100
core_radius 0.008
core_radius_kpc 0.052
cut_radius 3.091
cut_radius_kpc 19.710
v_disp 84.453
mag 22.950
z_lens 0.5440
end
potentiel 1291
profil 81
x_centre 52.459
y_centre 30.262
ellipticite 0.385
angle_pos 40.000
core_radius 0.008
core_radius_kpc 0.052
cut_radius 3.091
cut_radius_kpc 19.710
v_disp 84.453
mag 22.950
z_lens 0.5440
end
potentiel 95
profil 81
x_centre -6.777
y_centre -100.726
ellipticite 0.508
angle_pos 7.900
core_radius 0.008
core_radius_kpc 0.051
cut_radius 3.063
cut_radius_kpc 19.530
v_disp 84.065
mag 22.970
z_lens 0.5440
end
potentiel 95
profil 81
x_centre -6.777
y_centre -100.726
ellipticite 0.508
angle_pos 7.900
core_radius 0.008
core_radius_kpc 0.051
cut_radius 3.063
cut_radius_kpc 19.530
v_disp 84.065
mag 22.970
z_lens 0.5440
end
potentiel 1973
profil 81
x_centre 1.035
y_centre 110.933
ellipticite 0.000
angle_pos 10.200
core_radius 0.008
core_radius_kpc 0.051
cut_radius 3.049
cut_radius_kpc 19.440
v_disp 83.871
mag 22.980
z_lens 0.5440
end
potentiel 804
profil 81
x_centre 63.285
y_centre -12.358
ellipticite 0.835
angle_pos -67.300
core_radius 0.008
core_radius_kpc 0.051
cut_radius 3.049
cut_radius_kpc 19.440
v_disp 83.871
mag 22.980
z_lens 0.5440
end
potentiel 1098
profil 81
x_centre 16.388
y_centre 13.604
ellipticite 0.220
angle_pos 62.400
core_radius 0.008
core_radius_kpc 0.051
cut_radius 3.035
cut_radius_kpc 19.351
v_disp 83.678
mag 22.990
z_lens 0.5440
end
potentiel 1071
profil 81
x_centre 21.843
y_centre 9.908
ellipticite 0.280
angle_pos 79.900
core_radius 0.008
core_radius_kpc 0.051
cut_radius 3.021
cut_radius_kpc 19.262
v_disp 83.486
mag 23.000
z_lens 0.5440
end
potentiel 1005
profil 81
x_centre -21.408
y_centre 4.476
ellipticite 0.000
angle_pos 43.600
core_radius 0.008
core_radius_kpc 0.050
cut_radius 2.993
cut_radius_kpc 19.085
v_disp 83.102
mag 23.020
z_lens 0.5440
end
potentiel 1041
profil 81
x_centre 74.313
y_centre 6.903
ellipticite 0.220
angle_pos -7.600
core_radius 0.007
core_radius_kpc 0.046
cut_radius 2.768
cut_radius_kpc 17.648
v_disp 79.912
mag 23.190
z_lens 0.5440
end
potentiel 1671
profil 81
x_centre 49.952
y_centre 66.839
ellipticite 0.280
angle_pos 49.700
core_radius 0.007
core_radius_kpc 0.045
cut_radius 2.680
cut_radius_kpc 17.088
v_disp 78.635
mag 23.260
z_lens 0.5440
end
potentiel 1410
profil 81
x_centre -9.710
y_centre 39.387
ellipticite 0.220
angle_pos -68.200
core_radius 0.007
core_radius_kpc 0.045
cut_radius 2.655
cut_radius_kpc 16.931
v_disp 78.273
mag 23.280
z_lens 0.5440
end
potentiel 641
profil 81
x_centre -43.458
y_centre -28.904
ellipticite 0.000
angle_pos 71.600
core_radius 0.007
core_radius_kpc 0.045
cut_radius 2.655
cut_radius_kpc 16.931
v_disp 78.273
mag 23.280
z_lens 0.5440
end
potentiel 891
profil 81
x_centre 54.843
y_centre -6.060
ellipticite 0.220
angle_pos 19.100
core_radius 0.007
core_radius_kpc 0.043
cut_radius 2.583
cut_radius_kpc 16.470
v_disp 77.199
mag 23.340
z_lens 0.5440
end
potentiel 1634
profil 81
x_centre 1.410
y_centre 62.955
ellipticite 0.220
angle_pos 49.600
core_radius 0.007
core_radius_kpc 0.043
cut_radius 2.571
cut_radius_kpc 16.394
v_disp 77.022
mag 23.350
z_lens 0.5440
end
potentiel 756
profil 81
x_centre 88.150
y_centre -17.539
ellipticite 0.000
angle_pos -27.300
core_radius 0.007
core_radius_kpc 0.042
cut_radius 2.501
cut_radius_kpc 15.948
v_disp 75.965
mag 23.410
z_lens 0.5440
end
potentiel 728
profil 81
x_centre -88.528
y_centre -19.880
ellipticite 0.000
angle_pos 77.300
core_radius 0.006
core_radius_kpc 0.041
cut_radius 2.444
cut_radius_kpc 15.585
v_disp 75.095
mag 23.460
z_lens 0.5440
end
potentiel 1976
profil 81
x_centre -31.112
y_centre 111.352
ellipticite 0.220
angle_pos 65.200
core_radius 0.006
core_radius_kpc 0.040
cut_radius 2.411
cut_radius_kpc 15.371
v_disp 74.578
mag 23.490
z_lens 0.5440
end
potentiel 674
profil 81
x_centre -22.082
y_centre -24.760
ellipticite 0.471
angle_pos -8.000
core_radius 0.006
core_radius_kpc 0.040
cut_radius 2.399
cut_radius_kpc 15.300
v_disp 74.407
mag 23.500
z_lens 0.5440
end
potentiel 197
profil 81
x_centre -12.279
y_centre -80.981
ellipticite 0.220
angle_pos 22.800
core_radius 0.006
core_radius_kpc 0.038
cut_radius 2.291
cut_radius_kpc 14.612
v_disp 72.713
mag 23.600
z_lens 0.5440
end
potentiel 2082
profil 81
x_centre -4.956
y_centre 137.993
ellipticite 0.280
angle_pos -2.000
core_radius 0.006
core_radius_kpc 0.037
cut_radius 2.188
cut_radius_kpc 13.954
v_disp 71.058
mag 23.700
z_lens 0.5440
end
potentiel 2056
profil 81
x_centre 46.931
y_centre 131.403
ellipticite 0.471
angle_pos 39.000
core_radius 0.006
core_radius_kpc 0.036
cut_radius 2.168
cut_radius_kpc 13.826
v_disp 70.732
mag 23.720
z_lens 0.5440
end
potentiel 727
profil 81
x_centre 66.622
y_centre -20.593
ellipticite 0.280
angle_pos -17.200
core_radius 0.006
core_radius_kpc 0.036
cut_radius 2.138
cut_radius_kpc 13.636
v_disp 70.245
mag 23.750
z_lens 0.5440
end
potentiel 170
profil 81
x_centre 49.237
y_centre -86.439
ellipticite 0.471
angle_pos 26.400
core_radius 0.006
core_radius_kpc 0.035
cut_radius 2.109
cut_radius_kpc 13.449
v_disp 69.761
mag 23.780
z_lens 0.5440
end
potentiel 261
profil 81
x_centre 68.762
y_centre -74.096
ellipticite 0.724
angle_pos 57.400
core_radius 0.006
core_radius_kpc 0.035
cut_radius 2.090
cut_radius_kpc 13.326
v_disp 69.441
mag 23.800
z_lens 0.5440
end
potentiel 1696
profil 81
x_centre 108.320
y_centre 70.296
ellipticite 0.000
angle_pos 89.300
core_radius 0.005
core_radius_kpc 0.035
cut_radius 2.080
cut_radius_kpc 13.265
v_disp 69.281
mag 23.810
z_lens 0.5440
end
potentiel 431
profil 81
x_centre -57.808
y_centre -53.469
ellipticite 0.220
angle_pos -55.500
core_radius 0.005
core_radius_kpc 0.034
cut_radius 2.024
cut_radius_kpc 12.903
v_disp 68.330
mag 23.870
z_lens 0.5440
end
potentiel 1218
profil 81
x_centre 79.068
y_centre 21.544
ellipticite 0.280
angle_pos 29.600
core_radius 0.005
core_radius_kpc 0.031
cut_radius 1.871
cut_radius_kpc 11.932
v_disp 65.707
mag 24.040
z_lens 0.5440
end
potentiel 1920
profil 81
x_centre -40.004
y_centre 102.934
ellipticite 0.600
angle_pos 62.000
core_radius 0.005
core_radius_kpc 0.031
cut_radius 1.871
cut_radius_kpc 11.932
v_disp 65.707
mag 24.040
z_lens 0.5440
end
potentiel 510
profil 81
x_centre -66.402
y_centre -43.523
ellipticite 0.280
angle_pos -50.200
core_radius 0.005
core_radius_kpc 0.031
cut_radius 1.845
cut_radius_kpc 11.768
v_disp 65.255
mag 24.070
z_lens 0.5440
end
potentiel 1575
profil 81
x_centre 28.066
y_centre 55.607
ellipticite 0.471
angle_pos -55.600
core_radius 0.005
core_radius_kpc 0.030
cut_radius 1.795
cut_radius_kpc 11.447
v_disp 64.360
mag 24.130
z_lens 0.5440
end
potentiel 1159
profil 81
x_centre -78.664
y_centre 17.127
ellipticite 0.280
angle_pos 60.900
core_radius 0.005
core_radius_kpc 0.030
cut_radius 1.762
cut_radius_kpc 11.238
v_disp 63.770
mag 24.170
z_lens 0.5440
end
potentiel 1103
profil 81
x_centre 73.550
y_centre 12.772
ellipticite 0.000
angle_pos -3.400
core_radius 0.004
core_radius_kpc 0.029
cut_radius 1.699
cut_radius_kpc 10.832
v_disp 62.606
mag 24.250
z_lens 0.5440
end
potentiel 646
profil 81
x_centre 127.339
y_centre -28.463
ellipticite 0.849
angle_pos -9.800
core_radius 0.004
core_radius_kpc 0.027
cut_radius 1.637
cut_radius_kpc 10.440
v_disp 61.463
mag 24.330
z_lens 0.5440
end
potentiel 2083
profil 81
x_centre -4.479
y_centre 138.544
ellipticite 0.600
angle_pos 78.500
core_radius 0.004
core_radius_kpc 0.027
cut_radius 1.630
cut_radius_kpc 10.392
v_disp 61.322
mag 24.340
z_lens 0.5440
end
potentiel 2116
profil 81
x_centre 11.719
y_centre 155.161
ellipticite 0.600
angle_pos -83.100
core_radius 0.004
core_radius_kpc 0.025
cut_radius 1.486
cut_radius_kpc 9.478
v_disp 58.562
mag 24.540
z_lens 0.5440
end
potentiel 1193
profil 81
x_centre -134.574
y_centre 20.007
ellipticite 0.600
angle_pos 82.600
core_radius 0.004
core_radius_kpc 0.025
cut_radius 1.473
cut_radius_kpc 9.391
v_disp 58.293
mag 24.560
z_lens 0.5440
end
potentiel 901
profil 81
x_centre 11.621
y_centre -5.548
ellipticite 0.000
angle_pos -27.100
core_radius 0.004
core_radius_kpc 0.024
cut_radius 1.459
cut_radius_kpc 9.305
v_disp 58.025
mag 24.580
z_lens 0.5440
end
potentiel 611
profil 81
x_centre -19.114
y_centre -32.280
ellipticite 0.000
angle_pos -2.100
core_radius 0.004
core_radius_kpc 0.024
cut_radius 1.459
cut_radius_kpc 9.305
v_disp 58.025
mag 24.580
z_lens 0.5440
end
potentiel 1086
profil 81
x_centre -123.761
y_centre 11.173
ellipticite 0.385
angle_pos 17.200
core_radius 0.004
core_radius_kpc 0.023
cut_radius 1.387
cut_radius_kpc 8.845
v_disp 56.574
mag 24.690
z_lens 0.5440
end
potentiel 152
profil 81
x_centre -18.145
y_centre -90.572
ellipticite 0.385
angle_pos 28.100
core_radius 0.004
core_radius_kpc 0.023
cut_radius 1.381
cut_radius_kpc 8.804
v_disp 56.443
mag 24.700
z_lens 0.5440
end
potentiel 1669
profil 81
x_centre 15.752
y_centre 65.861
ellipticite 0.385
angle_pos 61.400
core_radius 0.004
core_radius_kpc 0.023
cut_radius 1.343
cut_radius_kpc 8.564
v_disp 55.669
mag 24.760
z_lens 0.5440
end
potentiel 1277
profil 81
x_centre -127.867
y_centre 28.270
ellipticite 0.724
angle_pos -47.300
core_radius 0.003
core_radius_kpc 0.022
cut_radius 1.300
cut_radius_kpc 8.293
v_disp 54.779
mag 24.830
z_lens 0.5440
end
potentiel 1951
profil 81
x_centre -36.339
y_centre 106.443
ellipticite 0.600
angle_pos 58.900
core_radius 0.003
core_radius_kpc 0.021
cut_radius 1.242
cut_radius_kpc 7.919
v_disp 53.532
mag 24.930
z_lens 0.5440
end
potentiel 200
profil 81
x_centre 54.676
y_centre -80.840
ellipticite 0.000
angle_pos -25.900
core_radius 0.003
core_radius_kpc 0.020
cut_radius 1.192
cut_radius_kpc 7.598
v_disp 52.434
mag 25.020
z_lens 0.5440
end
potentiel 390
profil 81
x_centre 86.720
y_centre -60.453
ellipticite 0.385
angle_pos -10.000
core_radius 0.003
core_radius_kpc 0.020
cut_radius 1.186
cut_radius_kpc 7.563
v_disp 52.313
mag 25.030
z_lens 0.5440
end
potentiel 1750
profil 81
x_centre -75.915
y_centre 75.846
ellipticite 0.600
angle_pos 24.500
core_radius 0.003
core_radius_kpc 0.019
cut_radius 1.159
cut_radius_kpc 7.391
v_disp 51.715
mag 25.080
z_lens 0.5440
end
potentiel 1601
profil 81
x_centre -16.092
y_centre 57.591
ellipticite 0.385
angle_pos 73.900
core_radius 0.003
core_radius_kpc 0.019
cut_radius 1.138
cut_radius_kpc 7.256
v_disp 51.240
mag 25.120
z_lens 0.5440
end
potentiel 973
profil 81
x_centre 15.636
y_centre 0.609
ellipticite 0.385
angle_pos -17.000
core_radius 0.003
core_radius_kpc 0.019
cut_radius 1.138
cut_radius_kpc 7.256
v_disp 51.240
mag 25.120
z_lens 0.5440
end
potentiel 258
profil 81
x_centre 19.172
y_centre -74.899
ellipticite 0.385
angle_pos 44.200
core_radius 0.003
core_radius_kpc 0.019
cut_radius 1.117
cut_radius_kpc 7.124
v_disp 50.771
mag 25.160
z_lens 0.5440
end
potentiel 231
profil 81
x_centre -31.227
y_centre -77.756
ellipticite 0.000
angle_pos 11.100
core_radius 0.003
core_radius_kpc 0.018
cut_radius 1.047
cut_radius_kpc 6.679
v_disp 49.160
mag 25.300
z_lens 0.5440
end
potentiel 833
profil 81
x_centre 41.778
y_centre -11.943
ellipticite 0.000
angle_pos -49.100
core_radius 0.003
core_radius_kpc 0.017
cut_radius 1.019
cut_radius_kpc 6.497
v_disp 48.486
mag 25.360
z_lens 0.5440
end
potentiel 501
profil 81
x_centre 108.469
y_centre -44.859
ellipticite 0.385
angle_pos 57.400
core_radius 0.003
core_radius_kpc 0.017
cut_radius 1.000
cut_radius_kpc 6.378
v_disp 48.041
mag 25.400
z_lens 0.5440
end
potentiel 853
profil 81
x_centre -98.920
y_centre -10.610
ellipticite 0.385
angle_pos -31.500
core_radius 0.003
core_radius_kpc 0.016
cut_radius 0.977
cut_radius_kpc 6.233
v_disp 47.491
mag 25.450
z_lens 0.5440
end
potentiel 675
profil 81
x_centre 12.833
y_centre -25.358
ellipticite 0.000
angle_pos -72.400
core_radius 0.003
core_radius_kpc 0.016
cut_radius 0.969
cut_radius_kpc 6.176
v_disp 47.273
mag 25.470
z_lens 0.5440
end
potentiel 1594
profil 81
x_centre -29.240
y_centre 56.619
ellipticite 0.000
angle_pos 33.900
core_radius 0.002
core_radius_kpc 0.015
cut_radius 0.887
cut_radius_kpc 5.658
v_disp 45.249
mag 25.660
z_lens 0.5440
end
fini
\ No newline at end of file
diff --git a/src/chi_computation.cpp b/src/chi_computation.cpp
index 9ec9701..e56b960 100644
--- a/src/chi_computation.cpp
+++ b/src/chi_computation.cpp
@@ -1,209 +1,296 @@
//
#include <string.h>
#include "structure_hpc.hpp"
#include "delense_CPU_utils.hpp"
#ifdef __WITH_MPI
#include <mpi.h>
#include "mpi_check.h"
#endif
#define MAX(a,b) (((a)>(b))?(a):(b))
void
-chi_computation(double* chi, int* images_found, int* numimagesfound, struct point* imagesposition, const int MAXIMPERSOURCE, const int *nimages_strongLensing, const galaxy* images, const int nsets)
+chi_computation (double* chi, int* images_found, int* numimagesfound, struct point* imagesposition, const int MAXIMPERSOURCE, const int *nimages_strongLensing,
+ const galaxy* images, const int nsets)
{
- /*
- int MAXIMPERSOURCE;
- int numImages = 0;
- int maxImagesPerSource = 0;
- for( int source_id = 0; source_id < nsets; source_id ++)
- {
- numImages += nimages_strongLensing[source_id];
- maxImagesPerSource = MAX(nimages_strongLensing[source_id], maxImagesPerSource);
- }
- MAXIMPERSOURCE = maxImagesPerSource;
- */
- int nimagesfound [nsets][MAXIMPERSOURCE]; // number of images found from the theoretical sources
- struct point tim [MAXIMPERSOURCE]; // theoretical images (computed from sources)
- int index = 0;
- //
- for( int source_id = 0; source_id < nsets; source_id ++)
+ /*
+ int MAXIMPERSOURCE;
+ int numImages = 0;
+ int maxImagesPerSource = 0;
+ for( int source_id = 0; source_id < nsets; source_id ++)
+ {
+ numImages += nimages_strongLensing[source_id];
+ maxImagesPerSource = MAX(nimages_strongLensing[source_id], maxImagesPerSource);
+ }
+ MAXIMPERSOURCE = maxImagesPerSource;
+ */
+ int nimagesfound[nsets][MAXIMPERSOURCE]; // number of images found from the theoretical sources
+ int prediction_assigned[MAXIMPERSOURCE]; // if prediction has already been assigned
+ struct point tim[MAXIMPERSOURCE]; // theoretical images (computed from sources)
+ int index = 0;
+ //
+ for (int source_id = 0; source_id < nsets; source_id++)
+ {
+ unsigned short int nimages = nimages_strongLensing[source_id];
+ //
+ //______________________Initialisation______________________
+ //
+ //for (unsigned short int i = 0; i < nimages; ++i)
+ // nimagesfound[source_id][i] = 0;
+ memset (&nimagesfound[source_id][0], 0, nimages * sizeof(int));
+ memset (&prediction_assigned[0], 0, MAXIMPERSOURCE * sizeof(int));
+ //
+ //@todo if(numimagesfound[source_id] > 0) ---> check against no img needed
+
+ // Creating image_distance matrix
+ double image_dist[MAXIMPERSOURCE][MAXIMPERSOURCE]; //[constraints][predictions]
+ struct point prediction_position;
+ int image_index;
+ //looping through predictions
+ for (int ii = 0; ii < numimagesfound[source_id]; ++ii)
{
- unsigned short int nimages = nimages_strongLensing[source_id];
- //
- //______________________Initialisation______________________
- //
- //for (unsigned short int i = 0; i < nimages; ++i)
- // nimagesfound[source_id][i] = 0;
- memset(&nimagesfound[source_id][0], 0, nimages*sizeof(int));
- //
- //for( unsigned short int image_id = 0; image_id < nimages; image_id++)
- //{
- //#endif //collapsing
- //
- double image_dist[MAXIMPERSOURCE];
- //
- //printf(" Source %d, number of images found %d\n", source_id, numimagesfound[source_id]);
- //
- struct point image_position;
- int image_index;
- //
- for (int ii = 0; ii < /*loc*/numimagesfound[source_id]; ++ii)
- {
- //
- //
- int image_index = 0;
- //
- image_position = imagesposition[source_id*MAXIMPERSOURCE + ii];
- image_dist[0] = mychi_dist(image_position, images[index + 0].center); // get the distance to the real image
- //@@printf("\nsource = %d n images found = %d, n images = %d,\n --> pos x y = %f %f\n", source_id, numimagesfound[source_id], nimages_strongLensing[source_id], image_position.x, image_position.y); fflush(stdout);
- //@@ printf(" ix iy = %f %f, dist = %f\n", images[index].center.x, images[index].center.y, mychi_dist(image_position, images[index].center));
- for(int i = 1; i < nimages_strongLensing[source_id]; i++)
- {
- // get the distance to each real image and keep the index of the closest real image
- image_dist[i] = mychi_dist(image_position, images[index + i].center);
- //@@printf(" ix iy = %f %f, dist = %f\n", images[index+i].center.x, images[index+i].center.y, image_dist[i]);
- if (image_dist[i] < image_dist[image_index])
- {
- image_index = i;
- }
- }
- //@@printf(" --> image_index = %d\n", image_index); fflush(stdout);
- //
- // we should exit loops here
- //
- // p1_time += myseconds();
- //
- int skip_image = 0;
- // Sometimes due to the numerical errors at the centerpoint,
- // for SIE potentials an additional image will appear at the center of the Potential.
- // This is due to the fact that it is not possible to simulate an infinity value
- // at the center correctly, Check that sis correspond to Nlens[0]
- /*
- for (int iterator = 0; iterator < runmode->Nlens[0]; ++iterator)
- {
- if ( fabs(image_position.x - lens[0].position_x[iterator]) <= dx/2. and fabs(image_position.y - lens[0].position_y[iterator]) <= dx/2.)
- {
- skip_image = 1;
- printf("WARNING: You are using SIE potentials. An image to close to one of the potential centers has been classified as numerical error and removed \n");
- }
- }*/
- //printf("%d %d %d %d\n", x_id, y_id,thread_found_image, skip_image);
- if (!skip_image)
- {
- //#pragma omp atomic
- (*images_found)++;
- struct point temp;
- //printf(" source %d, image %d, index %d, Images found: %d\n", source_id, image_id, image_index, nimagesfound[source_id][image_id][image_index]);
- //checking whether a closest image has already been found
- if (nimagesfound[source_id][image_index] == 0)
- { // if no image found up to now
+ (*images_found)++;
+ prediction_assigned[ii] = 0;
+ //
+ prediction_position = imagesposition[source_id * MAXIMPERSOURCE + ii];
+ image_dist[0][ii] = mychi_dist (prediction_position, images[index + 0].center);
+ printf ("\nsource = %d n images found = %d, n images = %d,\n --> pos x y = %f %f\n", source_id, numimagesfound[source_id],
+ nimages_strongLensing[source_id], prediction_position.x, prediction_position.y);
+ fflush (stdout);
+ //looping through constraints
+ printf (" ix iy = %f %f, dist = %f\n", images[index + 0].center.x, images[index + 0].center.y, image_dist[0][ii]);
+ for (int jj = 1; jj < nimages; jj++)
+ {
+ // get the distance to each real image and keep the index of the closest real image
+ image_dist[jj][ii] = mychi_dist (prediction_position, images[index + jj].center);
+ printf (" ix iy = %f %f, dist = %f\n", images[index + jj].center.x, images[index + jj].center.y, image_dist[jj][ii]);
+ }
+ }
- //image position is allocated to theoretical image
- //#pragma omp critical
- tim[image_index] = image_position;
- //#pragma omp atomic
- nimagesfound[source_id][image_index]++;
- }
- else if (nimagesfound[source_id][image_index] > 0)
- { // if we have already found an image
- // If the new image we found is closer than the previous image
- //printf(" tim2: %f %f\n", image_dist[image_index], mychi_dist(images[index + image_index].center, tim[image_index]));
- if (image_dist[image_index] < mychi_dist(images[index + image_index].center, tim[image_index]))
- {
- temp = tim[image_index]; // we store the position of the old image in temp
- //#pragma omp critical
- tim[image_index] = image_position; // we link the observed image with the image we just found
- //printf("tim2 %d %d = %f %f\n", image_id, image_index, image_position.x, image_position.y);
- }
- else
- {
- temp = image_position; // we store the position of the image we just found in temp
- }
- // initialising second_closest_id to the highest value
- // Loop over all images in the set except the closest one
- // and initialize to the furthest away image
- int second_closest_id = 0;
- for (int i = 1; i < nimages_strongLensing[source_id] && (i != image_index); i++)
- {
- if(image_dist[i] > image_dist[second_closest_id]) second_closest_id=i;
- }
- ///////////////////////////////////////////////////////////////
- // Loop over all images in the set that are not yet allocated to a theoretical image
- // and allocate the closest one
- // we search for an observed image not already linked (nimagesfound=0)
- for(int i = 0; i < nimages_strongLensing[source_id] && nimagesfound[source_id][i] == 0; i++)
- {
- if(image_dist[i] < image_dist[second_closest_id])
- {
- second_closest_id = i;
- // im_index value changes only if we found a not linked yet image
- image_index = i;
- //#pragma omp critical
- tim[image_index] = temp; // if we found an observed and not already linked image, we allocate the theoretical image temp
- }
- }
- // increasing the total number of images found (If we find more than 1 theoretical image linked to 1 real image,
- // these theoretical
- //#pragma omp atomic
- nimagesfound[source_id][image_index]++;
- // images are included in this number)
- }
- }
- //#pragma omp atomic
- //loc_images_found++;
- //thread_found_image = 0; // for next iteration
+ //assigning preditions to contraints
+ //looping through constraints
+ for (int jj = 0; jj < nimages; jj++)
+ {
+ int prediction_index = -1;
+ //find first non assigned prediction
+ for (int ii = 0; ii < numimagesfound[source_id]; ++ii)
+ {
+ printf(" cont %d predic %d : prediction_index %d %d\n", jj, ii, prediction_index,prediction_assigned[ii]);
+ if (prediction_assigned[ii] == 0)
+ {
+ prediction_index = ii;
+ prediction_assigned[ii] = 1;
+ break;
}
- //
- //}
- //#pragma omp barrier
- //____________________________ end of image loop
- //
- //____________________________ computing the local chi square
- //
- double chiimage;
- //
- int _nimages = nimages_strongLensing[source_id];
- //
- for( int iter = 0; iter < _nimages; iter++)
+ }
+ //If all predictions have not already been assigned
+ if (prediction_index != -1)
+ {
+ //looping through predictions
+ for (int ii = 1; ii < numimagesfound[source_id]; ++ii)
+ {
+ printf(" cont %d predic %d : prediction_index %d %d %f %f\n", jj, ii, prediction_index,prediction_assigned[ii] , image_dist[jj][ii], image_dist[jj][prediction_index]);
+ if (image_dist[jj][ii] < image_dist[jj][prediction_index] and prediction_assigned[ii] == 0)
+ {
+
+ /* Uncertain about if that should be done ...
+ //checking if prediction is indeed the closest // wrong should be if not already assigned
+ bool none_closer = true;
+ for (int kk = 0; kk < nimages; kk++)
+ {
+ printf(" cont %d predic %d %d : %f %f\n", jj, ii, kk, image_dist[jj][ii], image_dist[kk][ii]);
+ if (image_dist[jj][ii] > image_dist[kk][ii])
+ {
+ none_closer = false;
+ }
+ }
+ if(none_closer == true){
+ prediction_index = ii ;
+ prediction_assigned[ii] +=1;
+ }*/
+ prediction_assigned[prediction_index] = 0;
+ prediction_index = ii;
+ prediction_assigned[ii] = 1;
+ }
+ }
+ //write position
+ tim[jj] = imagesposition[source_id * MAXIMPERSOURCE + prediction_index];
+ nimagesfound[source_id][jj]++;
+ printf (" --> tim%d = %f %f\n", jj, tim[jj].x, tim[jj].y);
+ }
+ }
+
+#if 0
+
+ struct point image_position;
+ int image_index;
+ //
+ for (int ii = 0; ii < /*loc*/numimagesfound[source_id]; ++ii)
+ {
+ //
+ //
+ int image_index = 0;
+ //
+ image_position = imagesposition[source_id*MAXIMPERSOURCE + ii];
+ image_dist[0] = mychi_dist(image_position, images[index + 0].center);// get the distance to the real image
+ printf("\nsource = %d n images found = %d, n images = %d,\n --> pos x y = %f %f\n", source_id, numimagesfound[source_id], nimages_strongLensing[source_id], image_position.x, image_position.y); fflush(stdout);
+ printf(" ix iy = %f %f, dist = %f\n", images[index].center.x, images[index].center.y, mychi_dist(image_position, images[index].center));
+ for(int i = 1; i < nimages_strongLensing[source_id]; i++)
+ {
+ // get the distance to each real image and keep the index of the closest real image
+ image_dist[i] = mychi_dist(image_position, images[index + i].center);
+ printf(" ix iy = %f %f, dist = %f\n", images[index+i].center.x, images[index+i].center.y, image_dist[i]);
+ if (image_dist[i] < image_dist[image_index])
{
- int i = iter;
- //int i=iter/nimages_strongLensing[source_id];
- //int j=iter%nimages_strongLensing[source_id];
+ image_index = i;
+ }
+ }
+ //@@printf(" --> image_index = %d\n", image_index); fflush(stdout);
+ //
+ // we should exit loops here
+ //
+ // p1_time += myseconds();
+ //
+ int skip_image = 0;
+ // Sometimes due to the numerical errors at the centerpoint,
+ // for SIE potentials an additional image will appear at the center of the Potential.
+ // This is due to the fact that it is not possible to simulate an infinity value
+ // at the center correctly, Check that sis correspond to Nlens[0]
+ /*
+ for (int iterator = 0; iterator < runmode->Nlens[0]; ++iterator)
+ {
+ if ( fabs(image_position.x - lens[0].position_x[iterator]) <= dx/2. and fabs(image_position.y - lens[0].position_y[iterator]) <= dx/2.)
+ {
+ skip_image = 1;
+ printf("WARNING: You are using SIE potentials. An image to close to one of the potential centers has been classified as numerical error and removed \n");
+ }
+ }*/
+ //printf("%d %d %d %d\n", x_id, y_id,thread_found_image, skip_image);
+ if (!skip_image)
+ {
+ //#pragma omp atomic
+ (*images_found)++;
+ struct point temp;
+ printf(" nimagesfound %d\n", source_id, image_index, nimagesfound[source_id][image_index]);
+ //checking whether a closest image has already been found
+ if (nimagesfound[source_id][image_index] == 0)
+ { // if no image found up to now
+
+ //image position is allocated to theoretical image
+ //#pragma omp critical
+ tim[image_index] = image_position;
+ //#pragma omp atomic
+ nimagesfound[source_id][image_index]++;
+ }
+ else if (nimagesfound[source_id][image_index] > 0)
+ { // if we have already found an image
+
+ // If the new image we found is closer than the previous image
+ printf(" tim2: %f %f\n", image_dist[image_index], mychi_dist(images[index + image_index].center, tim[image_index]));
+ if (image_dist[image_index] < mychi_dist(images[index + image_index].center, tim[image_index]))
+ {
+ temp = tim[image_index]; // we store the position of the old image in temp
+ //#pragma omp critical
+ tim[image_index] = image_position;// we link the observed image with the image we just found
+ printf("on remplace la vieille img %d = %f %f\n", image_index, temp.x, temp.y);
+ //il faut recalculer les distances pr la vieille image
+ image_dist[0] = mychi_dist(temp, images[index + 0].center);// get the distance to the real image
+ printf(" ix iy = %f %f, dist = %f\n", images[index].center.x, images[index].center.y, image_dist[0]);
+ for(int i = 1; i < nimages_strongLensing[source_id]; i++)
+ {
+ // get the distance to each real image and keep the index of the closest real image
+ image_dist[i] = mychi_dist(temp, images[index + i].center);
+ printf(" ix iy = %f %f, dist = %f\n", images[index+i].center.x, images[index+i].center.y, image_dist[i]);
+ }
+ }
+ else
+ {
+ printf("on garde la vieille img %d = %f %f\n", image_index, tim[image_index].x, tim[image_index].y);
+ temp = image_position; // we store the position of the image we just found in temp
+ }
- //if(i != j)
- //{
- if(nimagesfound[source_id][i] > 0)
+ // initialising second_closest_id to the highest value
+ // Loop over all images in the set except the closest one
+ // and initialize to the furthest away image
+ int second_closest_id = 0;
+ ///////////////////////////////////////////////////////////////
+ // Loop over all images in the set that are not yet allocated to a theoretical image
+ // and allocate the closest one
+ for(int i = 0; i < nimages_strongLensing[source_id]; i++)
+ {
+ printf("second_closest_id %d de %f %f %d %d %d %f %f\n", second_closest_id, temp.x, temp.y, i, nimages_strongLensing[source_id],image_index,image_dist[i] , image_dist[second_closest_id]);
+ if(image_dist[i] <= image_dist[second_closest_id] and i != image_index)
{
- double pow1 = images[index + i].center.x - tim[i].x;
- double pow2 = images[index + i].center.y - tim[i].y;
- //
- //chiimage = pow(images[index + j].center.x - tim[i][j].x, 2) + pow(images[index + j].center.y - tim[i][j].y, 2); // compute the chi2
- chiimage = pow1*pow1 + pow2*pow2; // compute the chi2
- //printf("%d %d = %.15f\n", i, j, chiimage);
- (*chi) += chiimage;
+ second_closest_id = i;
+ // im_index value changes only if we found a not linked yet image
+ image_index = i;
+ //#pragma omp critical
+ tim[image_index] = temp;// if we found an observed and not already linked image, we allocate the theoretical image temp
}
- else
- //if(nimagesfound[source_id][i] == 0)
- {
- // If we do not find a correpsonding image, we add a big value to the chi2 to disfavor the model
- (*chi) += 1000.*nimages_strongLensing[source_id];
- printf(" source_id = %d, image number = %d has not found an image, chi = %f\n", source_id, i, *chi);
- }
- //}
+ }
+ // increasing the total number of images found (If we find more than 1 theoretical image linked to 1 real image,
+ // these theoretical
+ //#pragma omp atomic
+ nimagesfound[source_id][image_index]++;
+ // images are included in this number)
}
- //
- //____________________________ end of computing the local chi square
- //
- //printf("%d: chi = %.15f\n", source_id, *chi);
- /*
- for (int i=0; i < nimages_strongLensing[source_id]; ++i){
- for (int j=0; j < nimages_strongLensing[source_id]; ++j){
- printf(" %d",nimagesfound[i][j]);
- }
- printf("\n");
- }*/
+ }
+ //#pragma omp atomic
+ //loc_images_found++;
+ //thread_found_image = 0; // for next iteration
+ }
+#endif
+ //
+ //}
+ //#pragma omp barrier
+ //____________________________ end of image loop
+ //
+ //____________________________ computing the local chi square
+ //
+ double chiimage;
+ //
+ int _nimages = nimages_strongLensing[source_id];
+ //
+ for (int iter = 0; iter < _nimages; iter++)
+ {
+ int i = iter;
+ //int i=iter/nimages_strongLensing[source_id];
+ //int j=iter%nimages_strongLensing[source_id];
- //Incrementing Index: Images already treated by previous source_id
- index += nimages_strongLensing[source_id];
+ //if(i != j)
+ //{
+ if (nimagesfound[source_id][i] > 0)
+ {
+ double pow1 = images[index + i].center.x - tim[i].x;
+ double pow2 = images[index + i].center.y - tim[i].y;
+ //
+ //chiimage = pow(images[index + j].center.x - tim[i][j].x, 2) + pow(images[index + j].center.y - tim[i][j].y, 2); // compute the chi2
+ chiimage = pow1 * pow1 + pow2 * pow2; // compute the chi2
+ //printf("%d %d = %.15f\n", i, j, chiimage);
+ (*chi) += chiimage;
+ printf (" -> chi = %f\n", *chi);
+ }
+ else
+ //if(nimagesfound[source_id][i] == 0)
+ {
+ // If we do not find a correpsonding image, we add a big value to the chi2 to disfavor the model
+ (*chi) += 1000. * nimages_strongLensing[source_id];
+ printf (" source_id = %d, image number = %d has not found an image, chi = %f\n", source_id, i, *chi);
+ }
+ //}
}
+ //
+ //____________________________ end of computing the local chi square
+ //
+ //printf("%d: chi = %.15f\n", source_id, *chi);
+ /*
+ for (int i=0; i < nimages_strongLensing[source_id]; ++i){
+ for (int j=0; j < nimages_strongLensing[source_id]; ++j){
+ printf(" %d",nimagesfound[i][j]);
+ }
+ printf("\n");
+ }*/
+
+ //Incrementing Index: Images already treated by previous source_id
+ index += nimages_strongLensing[source_id];
+ }
}