29295 ответов в теме

[Karmael]

зависит от железа ессно, а так ограничения прописываются в конфиг симуляции, пр дефолту 8 часов, потом останов.

тонкость заключается в том, что есть понятие снапшота, сиречь снимка текущего состояния. чем чаще снапшоты тем больше дискретность, состояния межу снапштоами - не сохраняются.

В самый мелкий возможный шаг на сколько я понял 0.0025 последнее творчество генерил с шагом 0.01 для пушей плавности и нагенерил гигов 500 наверное состояний. В ранних видосах, я этого не делал, и в итоге у меня на z=1 было пять снапшотов, остальное это кручение камеры.

Так вот, z тоже задается в конфиге, по дефолту = 1, но это скучно, т.к. формирование галактик и прочего интересного начинается сильно позже. 

[Xiemargl]

Ну а я тут того, этого, оставшись со сраными бинарниками один на один, челендж окцептнул, и у меня новый СТРОНГЛИ ОРИДЖИНАЛ КОНТЕНТ:

 

Спойлер

ну после того, как размер картинки подрос до гигабайт, ессно ни о каких анимейтед гифах речи быть немогет, и посему ведоз!

 

При просмотре данного ведоза, заскучавший читатель вдруг встрепенеца и спросит, а хренли у тебя там дыр черных больше чем звезд?

И найдется что ответить автор, - нихрена. т.к. картинки строятся по плотности вещества, очевидная гиперплотность черных дыр, очевидна, и соотвественно они на изображениях - ярче. Что вот правда удивительно, что черные дыры начали образовываться раньше звезд, а это значит, что теория о "тихом" возникновении черных дыр при коллапсе газовых скоплений достаточной массы, - стала мейнстримом, и теперь в симуляции. Я тут надулся и даже график формирования объектов построил. А что бы было понятно, z- это красное смещение, по сути, скорость с которой разлетается вселенная после БП, и как суть единственное правильное время, потому что другова для нас нет. 

 

 

видоз можно и нужно смотреть на замедлении

поскольку нихрена не понятно, лучше скажи - что у тебя за исходные данные и как потом обрабатываются ?

[Karmael]

Ах да, симуляцию можно останавливать и запускать снова, с нужного снапшота.

более того, это иногда необходимо, в силу каких то сложных внутренних багов, симуляция некрасиво падает шелестя про то, время нонкондишен, баста. Перезапуск с последнего снапа будет падать так же, и надо пулять -1 снап, пока не перестанет падать.

К слову, законченную симуляцию можно продолжить  поправив конфиг

поскольку нихрена не понятно, лучше скажи - что у тебя за исходные данные и как потом обрабатываются ?

 

ну какие в симуляции исходные данные? сгенеренные ессно, специальным генератором, который наполняет твою коробку расстояний ровным слоем неровного вещества, в в зависимости от рандомсида. Тем самым моделируя состояние вселеной сразу после БП, когда вещество разлеталось ровным плотным слоем,

с нулевой по 10 кадр это достаточно хорошо показывает.

 

Обрабатывается это все по последнему дыху космологии, с учетом звездного ветра, сотывания вселенной, звезды формируются, накапливают тяж. элементы,  взрываются роняя газ на соседей и дыры, что умеют в испарение и аккреацию. Сами процессы колапсов на столь маштабной модели напоминаю 20 на 20 мегапарсеков, - не считаются, но просиходят. 

Сообщение отредактировал Karmael: 22 December 2018 - 20:45

[Atrosha]

Моего прова тут дудосят. Быстро не получится. Докачивает всякие зависимости и не качает одну из них уже хрен знает сколько раз пытаюсь. На половине падает коннект.

Вот это вот: pfft-1.0.8-alpha2-fftw3.tar.gz

гигов 500 наверное

У меня с этим затык. Всего-то 100 гигов свободно...

[Karmael]

лол, ну с почином, читай что пишет и хочет, - удовлетворяй!

еще про начальные данные: Вот хайрез состояния на 0.0

В нем значится уложены газ, дарматтер и прочие штуки в количестве предполагаемом на этом этапе формировния вселенной.

т.к. насколько ровно распределено вещество было не ясно, то пакуем кучками разной плотности во имя рандома

 

У меня с этим затык. Всего-то 100 гигов свободно...

 

ну, если снапить не по барски, то хватит поиграться. 

К слову, лишний раз спешу напомнить, что интересное в этом всем начитается в интерпретации данных, извлечении и работы с ними. 

из советов по запуску, - лучшие результаты показала следующая схема:

1 mpi процесс на сокет

1 тред на кору

HT - off

Память жрет к слову тоже нормально при кубе 20x20х20 каждый mpi процесс сожрал по 9,8 Гб

[Atrosha]

Кармаэль он в папке /depends залей сюда или на гуглдиск или на яндекс, плиз. Он всего 4.3 MB весит.

Сообщение отредактировал Atrosha: 22 December 2018 - 21:54

[Karmael]

да собсвенно весь депендесес пожатый: https://dropmefiles.com/R3zRY

 

там два файла, каждый можно скачать отдельно нажав в "подробней"

Сообщение отредактировал Karmael: 22 December 2018 - 21:17

[Atrosha]

Спасипка.

Компильнулось.

выпадает с такой ошибкой. MPI и GSL установлены.
 

[atrosha-desktop:05120] [[INVALID],INVALID] ORTE_ERROR_LOG: A system-required executable either could not be found or was not executable by this user in file ess_singleton_module.c at line 458
[atrosha-desktop:05120] [[INVALID],INVALID] ORTE_ERROR_LOG: A system-required executable either could not be found or was not executable by this user in file ess_singleton_module.c at line 166
--------------------------------------------------------------------------
It looks like orte_init failed for some reason; your parallel process is
likely to abort.  There are many reasons that a parallel process can
fail during orte_init; some of which are due to configuration or
environment problems.  This failure appears to be an internal failure;
here's some additional information (which may only be relevant to an
Open MPI developer):

  orte_ess_init failed
  --> Returned value A system-required executable either could not be found or was not executable by this user (-126) instead of ORTE_SUCCESS
--------------------------------------------------------------------------
--------------------------------------------------------------------------
It looks like MPI_INIT failed for some reason; your parallel process is
likely to abort.  There are many reasons that a parallel process can
fail during MPI_INIT; some of which are due to configuration or environment
problems.  This failure appears to be an internal failure; here's some
additional information (which may only be relevant to an Open MPI
developer):

  ompi_mpi_init: ompi_rte_init failed
  --> Returned "A system-required executable either could not be found or was not executable by this user" (-126) instead of "Success" (0)
--------------------------------------------------------------------------
*** An error occurred in MPI_Init
*** on a NULL communicator
*** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort,
***    and potentially your MPI job)
[atrosha-desktop:5120] Local abort before MPI_INIT completed completed successfully, but am not able to aggregate error messages, and not able to guarantee that all other processes were killed!

Сообщение отредактировал Atrosha: 22 December 2018 - 21:31

[Karmael]

А сам по себе mpi работает? попробуй:

 

mpirun -np 2 echo "test"

 

должно два раза написать тест

[vmarkelov]

И тут врываюсь такой я, с удвоителем скорости интернета бесплатно, без смс

[Atrosha]

Разобрался. У меня не было установлено MPI. Жесть.

Залей сюда файл с параметрами. Я пока отключу HT и перезагружусь.

[Karmael]

Спойлер

 

paramfile.gadget

 

#  Relevant files
 
InitCondFile = output/IC
OutputDir = output
TreeCoolFile = ../TREECOOL_fg_june11
# OutputList = 0.1,0.11,0.12,0,1.13,0.14,0.15,0.2,0.299,0.3,0.33333,0.4,0.45,0.5,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1,1.05,1.1,1.15,1.2,1.25,1.3,1.35,1.4,1.45,1.5,1.55,1.6,1.65,1.7,1.75,1.8,1.85,1.9,1.95,2.0
OutputList = 0.0,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0,0.08,0,0.09,0.1,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.19,0.2,0.21,0.22,0.23,0.24,0.25,0.26,0.27,0.28,0.29,0.3,0.31,0.32,0.33,0.34,0.35,0.36,0.37,0.38,0.39,0.4,0.41,0.42,0.43,0.44,0.45,0.46,0.47,0.48,0.49,0.5,0.51,0.52,0.53,0.54,0.55,0.56,0.57,0.58,0.59,0.6,0.61,0.62,0.63,0.64,0.65,0.66,0.67,0.68,0.69,0.7,0.71,0.72,0.73,0.74,0.75,0.76,0.77,0.78,0.79,0.8,0.81,0.82,0.83,0.84,0.85,0.86,0.87,0.88,0.89,0.9,0.91,0.92,0.93,0.94,0.95,0.96,0.97,0.98,0.99,1,1.01,1.02,1.03,1.04,1.05,1.06,1.07,1.08,1.09,1.1,1.12,1.13,1.14,1.15,1.16,1.17,1.18,1.19,1.2,1.21,1.22,1.23,1.24,1.25,1.26,1.27,1.28,1.29,1.3,1.31,1.32,1.33,1.34,1.35,1.36,1.37,1.38,1.39,1.4,1.41,1.42,1.43,1.44,1.45,1.46,1.47,1.48,1.49,1.5,1.51,1.52,1.53,1.54,1.55,1.56,1.57,1.58,1.59,1.6,1.61,1.62,1.63,1.64,1.65,1.66,1.67,1.68,1.69,1.7,1.71,1.72,1.73,1.74,1.75,1.76,1.77,1.78,1.79,1.8,1.81,1.82,1.83,1.84,1.85,1.86,1.87,1.88,1.89,1.9,1.91,1.92,1.93,1.94,1.95,1.96,1.97,1.98,1.99,2,2.01,2.02,2.03,2.04,2.05,2.06,2.07,2.08,2.09,2.1,2.12,2.13,2.14,2.15,2.16,2.17,2.18,2.19,2.2,2.21,2.22,2.23,2.24,2.25,2.26,2.27,2.28,2.29,2.3,2.31,2.32,2.33,2.34,2.35,2.36,2.37,2.38,2.39,2.4,2.41,2.42,2.43,2.44,2.45,2.46,2.47,2.48,2.49,2.5,2.51,2.52,2.53,2.54,2.55,2.56,2.57,2.58,2.59,2.6,2.61,2.62,2.63,2.64,2.65,2.66,2.67,2.68,2.69,2.7,2.71,2.72,2.73,2.74,2.75,2.76,2.77,2.78,2.79,2.8,2.81,2.82,2.83,2.84,2.85,2.86,2.87,2.88,2.89,2.9,2.91,2.92,2.93,2.94,2.95,2.96,2.97,2.98,2.99,3.0,3.01,3.02,3.03,3.04,3.05,3.06,3.07,3.08,3.09,3.1,3.12,3.13,3.14,3.15,3.16,3.17,3.18,3.19,3.2,3.21,3.22,3.23,3.24,3.25,3.26,3.27,3.28,3.29,3.3,3.31,3.32,3.33,3.34,3.35,3.36,3.37,3.38,3.39,3.4,3.41,3.42,3.43,3.44,3.45,3.46,3.47,3.48,3.49,3.5,3.51,3.52,3.53,3.54,3.55,3.56,3.57,3.58,3.59,3.6,3.61,3.62,3.63,3.64,3.65,3.66,3.67,3.68,3.69,3.7,3.71,3.72,3.73,3.74,3.75,3.76,3.77,3.78,3.79,3.8,3.81,3.82,3.83,3.84,3.85,3.86,3.87,3.88,3.89,3.9,3.91,3.92,3.93,3.94,3.95,3.96,3.97,2.98,3.99,4.0,4.01,4.02,4.03,4.04,4.05,4.06,4.07,0,4.08,0,4.09,4.1,4.11,4.12,4.13,4.14,4.15,4.16,4.17,4.18,4.19,4.2,4.21,4.22,4.23,4.24,4.25,4.26,4.27,4.28,4.29,4.3,4.31,4.32,4.33,4.34,4.35,4.36,4.37,4.38,4.39,4.4,4.41,4.42,4.43,4.44,4.45,4.46,4.47,4.48,4.49,4.5,4.51,4.52,4.53,4.54,4.55,4.56,4.57,4.58,4.59,4.6,4.61,4.62,4.63,4.64,4.65,4.66,4.67,4.68,4.69,4.7,4.71,4.72,4.73,4.74,4.75,4.76,4.77,4.78,4.79,4.8,4.81,4.82,4.83,4.84,4.85,4.86,4.87,4.88,4.89,4.9,4.91,4.92,4.93,4.94,4.95,4.96,4.97,4.98,4.99,5.0,5.01,5.02,5.03,5.04,5.05,5.06,5.07,0,5.08,0,5.09,5.1,5.11,5.12,5.13,5.14,5.15,5.16,5.17,5.18,5.19,5.2,5.21,5.22,5.23,5.24,5.25,5.26,5.27,5.28,5.29,5.3,5.31,5.32,5.33,5.34,5.35,5.36,5.37,5.38,5.39,5.4,5.41,5.42,5.43,5.44,5.45,5.46,5.47,5.48,5.49,5.5,5.51,5.52,5.53,5.54,5.55,5.56,5.57,5.58,5.59,5.6,5.61,5.62,5.63,5.64,5.65,5.66,5.67,5.68,5.69,5.7,5.71,5.72,5.73,5.74,5.75,5.76,5.77,5.78,5.79,5.8,5.81,5.82,5.83,5.84,5.85,5.86,5.87,5.88,5.89,5.9,5.91,5.92,5.93,5.94,5.95,5.96,5.97,5.98,5.99
 
Nmesh = 256
 
# CPU time -limit
 
TimeLimitCPU = 1043000 #= 8 hours
 
 
# Code options
 
#  Characteristics of run
 
TimeMax = 6.00000
 
Omega0 = 0.2814      # Total matter density  (at z=0)
OmegaLambda = 0.7186      # Cosmological constant (at z=0)
OmegaBaryon = 0.0464     # Baryon density        (at z=0)
HubbleParam = 0.697      # Hubble paramater (may be used for power spec parameterization)
 
CoolingOn = 1
StarformationOn = 1
BlackHoleOn = 1
HydroOn = 1
StarformationCriterion = density
RadiationOn = 1
MassiveNuLinRespOn = 0
WindOn = 1
MetalCoolFile = ../cooling_metal_UVB
 
# Accuracy of time integration
MaxSizeTimestep = 0.1
MinSizeTimestep = 0.00
 
SnapshotWithFOF = 1
FOFHaloLinkingLength = 0.2
FOFHaloMinLength = 32
 
#  Further parameters of SPH
DensityKernelType = quintic
 
DensityContrastLimit = 100   # max contrast for hydro force calculation
DensityResolutionEta = 1.0  # for Cubic spline 1.0 = 33
MaxNumNgbDeviation = 2
ArtBulkViscConst = 0.75
InitGasTemp = 580.0        # always ignored if set to 0
MinGasTemp = 5.0
 
 
# Memory allocation
 
PartAllocFactor = 2.0
 
# Softening lengths
 
MinGasHsmlFractional = 0.01
 
#----------------------BH Stuff-------------------------
BlackHoleFeedbackFactor = 0.05
BlackHoleFeedbackRadius = 0.
BlackHoleFeedbackRadiusMaxPhys = 0.
BlackHoleFeedbackMethod = spline | mass
SeedBlackHoleMass = 5.0e-5
BlackHoleAccretionFactor = 100.0
BlackHoleNgbFactor = 2.0
BlackHoleEddingtonFactor = 3.0
 
MinFoFMassForNewSeed = 1
TimeBetweenSeedingSearch 1.03
 
#----------------------SFR Stuff-------------------------
 
CritPhysDensity = 0       #  critical physical density for star formation in
#  hydrogen number density in cm^(-3)
 
CritOverDensity = 57.7   #  overdensity threshold value
 
QuickLymanAlphaProbability  0 # Set to 1.0 to turn dense gas directly into stars.
 
MaxSfrTimescale = 1.5     # in internal time units
 
 
TempSupernova = 1.0e8   #  in Kelvin
 
TempClouds = 1000.0   #  in Kelvin
FactorSN = 0.1
FactorEVP = 1000.0
 
WindModel = ofjt10,decouple
WindEfficiency = 2.0
WindEnergyFraction = 1.0
WindSigma0 = 353.0 #km/s
WindSpeedFactor = 3.7
 
WindFreeTravelLength = 20
WindFreeTravelDensFac = 0.1
 
 

 

в paramfile.genic я только boxSize поменял на 20000

 

 

в общем по уму делать надо так, копировать папку exampels/hydro куданить в exampels/hydro.orig

потом править файлы параметров. параметры запуска устанавливаются в run.sh

 

что делать если обвалилось? 

 

из папки где работаешь, делаешь

 

mpirun -np 2 ../../gadget/MP-gadget 1 - перезапускает симуляцию с крайнего снапшота

mpirun -np 2 ../../gadget/MP-gadget 2 <номер снапшота>  - перезапускает с конкретного снапшота

mpirun -np 2 ../../gadget/MP-gadget -1  - перезапускает с начала

[Karmael]

а я чот потёк на трехмерных построениях, и пойду в дотан штоль прям подеградирую

[Atrosha]

Чота считает. Уже 270 метров насчитала какого-то говна.

mpirun -np 2 ../../gadget/MP-gadget 1 - перезапускает симуляцию с крайнего снапшота
mpirun -np 2 ../../gadget/MP-gadget 2 <номер снапшота>  - перезапускает с конкретного снапшота
mpirun -np 2 ../../gadget/MP-gadget -1  - перезапускает с начала

У меня однопроцессорный же.

[Karmael]

ну запустить то можно и два, лел, но на самом деле верно, смотри по железу. 

на счет говна, в нем то вся соль. начинай изучать matplot.pyplot ну и bigfile https://github.com/rainwoodman/bigfile

ну или чем ты там рисовать собрался. без этого посмотреть что унутре говна не выйдет

есть еще https://nbodykit.rea...s.io/en/latest/ , он тоже умеет в BigFileCatalg (это формат этого самого говна) 

короче, самое по сути интересное

[Atrosha]

А 3D визуализатора нет? Поищи, мож замутим тут чонить.

Ага. Вот этот pyplot умеет. Сделаем. Только посчитаю чутку, чтобы было с чем работать.

[Crulod]

какой смысл тратить время на такую ерунду.

не ну я понимаю старкалк с каталогом САО... а вот это вот - чо?

вы понимаете шо играть в такие игры на коклютерах - развлечение не але... особенно глядя на то сколько сил затратить надо шобы это запустить.

вам шо? делать нехрен?

[Atrosha]

Тудыть твою маковку!

Не считает в бэкграунде нихрена. Всю ночь провисела без дела походу! Всего 13 ГБ насчитала.

[Crulod]

не ну а чо ты хотел от куска окаменевшего говна мамонта?

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

ты чего его схватил? типа станешь грамотней если твой коклютер изнасилдуешь и он нащитает тебе всякой никчемной херни.

а тошо оно - никчемная херня, не вооруженным глазом даже сквозь интернет видно.

...угу... 22 мегапарсека! после одын зет!

"ну придурки!"(тм)