Программа использует пакет colsys для решения обыкновенных дифференциальных уравнений (краевая задача). Colsys решатель черного ящика; пользователь должен предоставить оду, граничные условия и якобианы

The purpose of the program is to find numerical wormhole solutions to the EsGB equations when a selfinteraction potential for the scalar field is present.

Цель программы - найти численные решения кротовых нор для уравнений EsGB при наличии потенциала самовзаимодействия для скалярного поля.

The program uses the COLSYS package for solutions of ordinary differential equations (boundary value problem). COLSYS is a black box solver; the user has to provide the ODEs, the boundary conditions and the Jacobians. See colnix.info for more information.

Программа использует пакет COLSYS для решения обыкновенных дифференциальных уравнений (краевая задача). COLSYS - решатель черного ящика; пользователь должен предоставить ОДУ, граничные условия и якобианы. См. Colnix.info для получения дополнительной информации.

Files:

colnix1.f : COLSYS package

colnix.info : usage of COLSYS

Файлы:

colnix1.f: пакет COLSYS

colnix.info: использование COLSYS

esgbe2.f : main program esgbe20subs.f : subroutines for the ODEs, boundary conditions and Jacobians esgbdetasubs.f: subroutine to compute the determinant deta esgbe2.in : input parameters esgbe2.data: output

esgbe2.f: основная программа esgbe20subs.f: подпрограммы для ODE, граничные условия и якобианы esgbdetasubs.f: подпрограмма для вычисления определителя deta esgbe2.in: входные параметры esgbe2.data: output

esgbe20_alp=1.5_pt=-1.0_l=0.06 : solution file

esgbe20_alp = 1.5_pt = -1.0_l = 0.06: файл решения

Compile the subroutines with the intel FORTRAN compiler

ifort -c colnix1.f

ifort -c esgbdetasubs.f

ifort -c esgbe20subs.f

(any other compiler will do)

Скомпилируйте подпрограммы с помощью компилятора Intel FORTRAN

ifort -c colnix1.f

ifort -c esgbdetasubs.f

ifort -c esgbe20subs.f

(подойдет любой другой компилятор)

compilation of the main program and linking the subroutines ifort -o esgbe2 esgbe2.f esgbe20subs.o esgbdetasubs.o colnix1.o run the program with (linux)

./esgbe2

компиляция основной программы и связывание подпрограмм ifort -o esgbe2 esgbe2.f esgbe20subs.o esgbdetasubs.o colnix1.o запустить программу с (linux)

./esgbe2

The input file esgbe2.in contains the values of the parameters:

alp bet icnd icfc

1.5 1.5 3 1

Входной файл esgbe2.in содержит значения параметров:

alp bet icnd icfc

1,5 1,5 3 1

values of alp=$\alpha$, bet=$\beta$, icnd = indicator of boundary conditions at the throat/equator 1 f_0(0) = ft0, 2 f_1(0) = ft1, 3 phi(0) = pt icfc not needed

значения alp = $ \ alpha $, bet = $ \ beta $, icnd = индикатор граничных условий в горловине / экваторе 1 f_0 (0) = ft0, 2 f_1 (0) = ft1, 3 phi (0) = pt icfc не нужен

ft0 f1t pt

-22.0 1.025 -1.0

values of ft0,f1t,pt

xlm,c2,c4,c6

0.06 1.1 -2.0 1.0

ft0 f1t pt

-22,0 1,025 -1,0

значения ft0, f1t, pt

XLM, C2, C4, C6

0,06 1,1 -2,0 1,0

parameters of the potential tolerance

1d-2

ipar7,ipar9

-1 3

параметры потенциального допуска

1д-2

ipar7, ipar9

-1 3

COLSYS parameters (use ipar9 = 2 to increase the number of subintervals)

Параметры COLSYS (используйте ipar9 = 2 для увеличения количества подынтервалов)

Note that for fixed $\alpha$ the families of solutions are not uniquely determined by $f_0(0)$ or $\phi(0)$. Sometimes it is necessary to switch the indicator icnd.

Обратите внимание, что при фиксированном $ \ alpha $ семейства решений не определяются однозначно с помощью $ f_0 (0) $ или $ \ phi (0) $. Иногда бывает необходимо переключить индикатор icnd.

The file esgbe20_alp=1.5_pt=-1.0_l=0.06 contains solution for parameters $\alpha=1.5$, $\phi(0)=-1$ and selfinteraction potential To use it as initial guess:

Файл esgbe20_alp = 1.5_pt = -1.0_l = 0.06 содержит решение для параметров $ \ alpha = 1.5 $, $ \ phi (0) = - 1 $ и потенциал самовзаимодействия. Чтобы использовать его в качестве первоначального предположения:

cp esgbe20_alp=1.5_pt=-1.0_l=0.06 guess.data and run the program (with linux)

./esgbe2

cp esgbe20_alp = 1.5_pt = -1.0_l = 0.06 guess.data и запустите программу (с linux)

./esgbe2

In order to produce a family of solutions the output from the previous run can be used as initial guess. (but see Note above) cp result.data guess.data./esgbe2

Чтобы создать семейство решений, результат предыдущего прогона можно использовать в качестве начального предположения. (но см. примечание выше) cp result.data guess.data./esgbe2

Note that the output result.data is in a format ready to be used by COLSYS. For plotting use instead all_data.

Обратите внимание, что выходные данные result.data находятся в формате, готовом для использования COLSYS. Для построения вместо этого используйте all_data.

The output file esgbe2.data contains the properties of the solutions:

$\alpha$, $\lambda$, c2=$c_2$, c4=$c_4$, c6=$c_6$, $f_1(0)$,

$\alpha$, $f_0(1)$, $f_1(1)$, $f'_1(1)$ $\phi(1)$, $\phi'_1(1)$,

$\alpha$, $f_0(0)$, $f_1(0)$, $f'_1(0)$ $\phi(0)$, $\phi'_1(0)$,

deta_m, x_dm

Выходной файл esgbe2.data содержит свойства решений:

$ \ alpha $, $ \ lambda $, c2 = $ c_2 $, c4 = $ c_4 $, c6 = $ c_6 $, $ f_1 (0) $,

$ \ alpha $, $ f_0 (1) $, $ f_1 (1) $, $ f'_1 (1) $ $ \ phi (1) $, $ \ phi'_1 (1) $,

$ \ alpha $, $ f_0 (0) $, $ f_1 (0) $, $ f'_1 (0) $ $ \ phi (0) $, $ \ phi'_1 (0) $,

deta_m, x_dm

where deta_m is the minimum of deta and x_dm the location of the minimum. The vanishing of deta indicates the presence of a cusp singularity.

Anytime the program is run, the results are added to the file.

где deta_m - это минимум deta, а x_dm - расположение минимума. Исчезновение deta указывает на наличие особенности излома.

Каждый раз, когда программа запускается, результаты добавляются в файл.