Re: Сортировка комплексных чисел?

From

Evgenij M. Baldin ()

To

Andrey Belyaev

Date

2002-08-08T18:26:39Z

Area

RU.ALGORITHMS

From: "Evgenij M. Baldin" <baldin@Slon.inp.nsk.su>

Добрый день

Andrey Belyaev <ols@metacom.ru> wrote:
> 1. Если я правильно понял для каждого набора в 20000 ты делаешь
> ищешь ближайшую к 20000 точек - то есть находишь ближайшую
> к каждой из точек (типа все со всеми) ?

есть два разных набора: один - набор событий, второй - матрица коррекции 
- подробно описано ниже

> 2. Не позволяет ли характер расположения точек произвести
> предварительный отброс процентов 95% точек формулой
> max(|x1-x2|+|y1-y2|) - то есть просто сложение разниц координат ?

как я понимаю основное время тратится на перебор - выключение вычисления 
корня дало прирост скорости только на 10% 

>> P.S. Следует ли подробно описать задачу - зачем? или это будет просто
>> избыточный шум?

> Физику можно не надо, а вот более подробно абстракцию - можно бы.

У меня есть наборы пар чисел (4 набора по 20-50 тыс после 
предварительного отбора), получаемых совершенно случайным образом.

Эти пары чисел, представляют из себя координаты (x,y) на плоскости 
(отсюда и аналогия с комплексными числами :) ) - каждое число пусть 
называется событием - это то что получается в жизни.

Процедура получения координат события требует коррекции. Для этой 
коррекции с помощью моделирования для различных точек плоскости (их было 
взято около 40 тыс) с помощью того же метода, что и в жизни был получен 
набор сопоставлений (x_sim,y_sim)->(x_r,y_r) - то есть, если 
удастся найти ближайшую точку из массива (x_r,y_r) к координате 
события (x,y), то получается как бы истинную координату события.

В идеале было бы вообще не находить ближайшую, а получить аппроксимацию 
от четырёх ближайших точек, например.

Матрица x_sim,y_sim регулярна, а вот восстанавливыемые координаты 
(x_r,y_r)  не очень :( - все попытки описать это какой-либо функцией 
окончились неудачей, расцепить координаты можно, но с потерей точности - 
кроме того симметрия в некоторых местах нарушается из-за нерегулярности 
структуры. 

Отображение абсолютно однозначно. 

С уважением
	Евгений


P.S. Координаты это координаты входа e+ и e- в детектор после рассеяния, а 
истинная координата нужна из-за сильной зависимости вероятности рассеяния 
от угла для поиска событий J/Psi->e+e- на фоне e+e- -> e+e- : 
эффект/фон=5%, поэтому всё надо вылизать


-- 
                  Budker Institute of Nuclear Physics
                     e-mail: E.M.Baldin@inp.nsk.su
	           WWW: http://www.inp.nsk.su/~baldin
--- ifmail v.2.15dev5
 * Origin: BINP, Novosibirsk, Russia (2:5020/400)