Crp.
4 ГОСТ 25645.104—83
( 4 r \) = h (?*)
M ’V
h
)
U r „ ),J H n
I »
Л » !
......
’V - l,
где (rt — номер пересекаемой лучом поверхности.
К„
(70, Б)
(70)
3.3.5. Из полученных индексов неопределенности для точки гп отбирают относящиеся к К~А
зоне и формируют совокупность векторов неопределенности(г„) (/(=!,....
3.3.6. Определяют последовательность номеров зон(л~1......Лг0>, в которых располо
жены отрезки луча l„ (r„, Q), путем отыскания одинаковых (равных)векторов среди а* (гк) и
«к (r„) (К= 1
......
Отсутствие таких векторов для некоторой точки щ, (N0<N)
свидетельст вует о ее расположении вне объекта и процесс идентификации отрезков t „Q) для
n>N0 пре
кращают.г
3.3.7. Вычисляют количество вещества на пути луча g (го, Q) по формуле
—* -►Wl—1 *►*♦—у —г
<(f«i Q)=2 /я(г0,2)-р* (г0, $2).
(б)
Л-1
4 7
3.4. Функции самоэкраниронамности ш, (gu г0) и экранированности защитой ш2 (£г, г„) долж
ны быть рассчитаны аналогично w (g, г0), причем для вычисления ig, (/£, Q) и g2 (ъ, Q) осущест
вляют раздельное суммирование расстояний, пройденных лучом в зонах фантома и защиты, умножен
ных на плотность вещества всоответствующих зонах.
3.5. Возможный способ реализации алгоритма приведен в рекомендуемом приложении.
ПРИЛОЖЕНИЕ
Рекомендуемое
СПОСОБ РЕАЛИЗАЦИИ АЛГОРИТМА
(ПРОГРАММА OPTIC)
I. Описание программы OPTIC
COMMO областей и путемдания формальн х параметров
•lift Программа OPTIC предназначена длн^расчета функций экранированности и- (g, ~г„),саыоэкрапнрованноегн
*t (5ьЛ>) и экранированное™ защитой г0) точек г0 о объектах сложной геометрической конфигурации с непостоян
нымн физическими свойствами но объему. Вычисление функции зкраиироаанносги, определяемой выражением (2) насто
ящего стандарта, осуществляется методом Монте-Карло. Программа написана на языке Фортран и ориентирована на ЭВМ
тнпа ЕС или БЭСМ-6. Время счета одного варианта зависит от сложности объекта и требуемой точности вычисляемого функ
ционала. Блок-схема вычисления функционалов и>.», впредставлена на черт. 1 (в левом углу блоков указаны номера,
соответствующие пунктам описания программы).
ния
Передача
N
информации между
за
подпрограммами
ы
и связь их с
.
управляющей программой осуществляется и виде описа
1.2. О п и сан и е COM MON - о б л а ст ей
1.2.1. COMMON /AG/ UR. VV. \VR. A (50. 10). RO (30),
где UR. VV. WR — рабочие ячейхи;
A (50, 10)— массив коэффициентов,оивсываюшлх поверхности (задастся в соответствии с таблицейнастоящего
стандарта;
RO (30) — плотность вещества а аоие
1.2.2. COMMON AG1/N, NZON, ГГ (50), NCF (50). MI (30). IPZ (30,6),ID (30,6). KPN (30).
где N<50 — количество поверхностей, применяемое для задания объекта;
NZON<30 — количество зон, применяемое для задания объекта (включая пустоты);
IT (1) <14,1-*), . . . . N —тип поверхности;
NCF (1)<Ю, I - I
...........
N - максимальное количество коэффициентов, необходимое для задания поверхности I;
М1(К)<6. К=1, . . . NZON — количество поверхностей, ограничивающих зону К;
IPZ (К, J)<50— порядковыйномер поверхности, ограничивающей К-ю зону ( К = > NZON; J~ 1........ М1(К)).
KPN (К) — индекс материала в зоне К (предполагается, что индекс KPN —1-, имеет вещество фантома).
1.2.3. COMMON /AG2/ IDI (100). IP (100). SP (100), КР (100).
где 1D1 (100) — массив рабочих ячеек;