Программирование. Рубрика в журнале - Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование

Статья научная
Данная статья посвящена проблеме подбора параметров при моделировании надежности программного обеспечения. В статье рассматривается процесс подбора параметров S-образной модели роста надежности ПО на основе данных об отказах в предыдущих релизах, проводится сравнение точности модели при разных способах подбора параметров. В качестве базовой модели, описанной в данной статье, используется модель Goel-Okumoto. Данная модель строится на основе данных об отказах в программе за определенный промежуток времени. Для того, чтобы получить адекватные оценки при помощи данной модели требуется определенное количество данных об отказах, которые недоступны до тех пор, пока система не была протестирована в течение достаточно долгого периода времени. В качестве экспериментальных данных используются статистические данные, собранные в течение трех последовательных релизов программного продукта промышленного масштаба. Оценка параметров модели была выполнена при помощи функции максимального правдоподобия.
Бесплатно

Статья научная
Рассматривается система линейных алгебраических уравнений с интервальной матрицей коэффициентов и интервальной правой частью. Для данных систем вводится понятие псевдорешения. Доказано существование псевдорешения для любых интервальных систем линейных уравнений, предложен способ поиска псевдорешения с помощью решения соответствующей задачи линейного программирования. Вследствие вырожденности полученной задачи для ее решения необходимо использовать вычисления, обеспечивающие точность, намного превышающую возможности стандартных типов данных языков программирования. Симплекс-метод в сочетании с безошибочными дробно-рациональными вычислениями дает решение задачи. Для реализации используется крупнозернистый параллелизм в распределенных системах на основе MPI. Для реализации безошибочных дробно-рациональных вычислений на GPU используется CUDA C.
Бесплатно

Статья научная
Проведено исследование влияния неопределенности в кинетических параметрах на результаты решения прямой и обратной задач химической кинетики. Кинетические данные представлены интервалами и рассматриваются как объекты интервального анализа. Разработан алгоритм решения прямой задачи методом интервального анализа чувствительности и алгоритм решения обратной задачи по вычислению областей неопределенности параметров. Проведен вычислительный эксперимент по поиску интервального вектора параметров на примере промышленно значимой реакции. Показано, что интервальное решение прямой задачи, полученное для математической модели реакции, удовлетворяет заданному предельно допустимому значению погрешности в экспериментальных замерах концентраций участвующих веществ.
Бесплатно

Статья научная
Задача планирования выполнения пакетов заданий в многостадийных системах при наличии ограничения на длительность интервалов времени ее функционирования является комплексной. Ее решение предполагает (с учетом требования формирования комплектов из результатов) определение составов пакетов, составов групп пакетов, выполняемых в течение временных интервалов заданной длительности, расписаний выполнения пакетов на приборах многостадийной системы. Для определения комплексных решений применен аппарат теории иерархических игр. Реализуется построение модели иерархической игры для принятия решений по составам пакетов, групп пакетов и расписаниям выполнения пакетов. В модели учтено требование формирования комплектов из результатов выполнения пакетов заданий. Задача определения составов групп пакетов является NP-трудной, поэтому для ее решения требуется применение приближенных методов оптимизации. Формулируются метод построения начального решения по составам групп пакетов и метод распределения результатов выполнения пакетов заданий по комплектам. Сформулирован способ построения новых решений по составам групп пакетов заданий. Введены условия, позволяющие определять исключаемые из групп пакеты на основе количества результатов выполнения заданий каждого типа, не включаемых в комплекты. Предложен метод локальной оптимизации решений по составам групп пакетов.
Бесплатно

Статья научная
В работе предложена неявная схема разрывного метода Галеркина для решения уравнений газовой динамики на неструктурированных сетках. Неявная схема основана на представлении системы сеточных уравнений в "дельта-форме". Для решения полученной в ходе аппроксимации исходных уравнений СЛАУ применяются решатели из библиотеки NVIDIA AmgX. Для верификации численного алгоритма был произведен расчет течения невязкого сжимаемого газа в плоском канале с клином и решена задача об обтекании симметричного профиля NACA0012. Проведено сравнение полученных результатов с результатами эксперимента и известными численными решениями представленных задач. Сделан вывод о хорошем совпадении численных и экспериментальных данных.
Бесплатно

Программирование графических процессов при помощи расширяемых языков
Статья научная
В статье рассматривается система программирования ГПУ NUDA, созданная на основе расширяемого языка Nemerle. Она дает программисту контроль над процессом переноса программы на ГПУ, переложив механическую работу на компилятор. Макросы и аннотации облегчают перенос программ на ГПУ и повышают их производительность без ущерба для размера и читаемости исходного кода. На ряде архитектур ГПУ и задач удалось добиться повышения производительности в несколько раз по сравнению с исходным вариантом.
Бесплатно

Статья научная
Задача определения механических свойств различных упругих материалов возникает во многих приложениях. В данной статье приводится описание программного комплекса SIMPLE, созданного для решения задачи определения эффективных параметров волновода по ультразвуковым сигналам, регистрируемым в отдельных точках поверхности. Данная версия программного комплекса предназначена для случая изотропного упругого слоя. Дается краткое описание математической модели, на которой базируется подход, и основных алгоритмов обработки экспериментальных данных и решения обратной задачи. Описывается также общая схема принципа функционирования SIMPLE на пользовательском и программном уровне. Приводятся примеры использования разработанного программного комплекса для определения модуля Юнга и коэффициента Пуассона образцов из стали и пластика.
Бесплатно

Программный комплекс TeDy для решения задач термодинамического моделирования
Статья научная
Для решения задач термодинамического моделирования РФЯЦ-ВНИИТФ им. академ. Е.И. Забабахина разрабатывает программный комплекс, обеспечивающий расчет химически равновесного состава в многокомпонентных многофазных системах, изменения термодинамических функций отдельных химических реакций, диаграмм состояния многофазных систем и доступ к базе термодинамических данных. В основе функции расчета равновесия лежит стехиометрический метод расчета суммарной энергии Гиббса системы в зависимости от координат реакций. Координаты минимума рассчитываются одним из численных методов оптимизации. На основе значений координат реакций, соответствующих минимуму энергии Гиббса, определяется равновесный вещественный состав. Расчет термодинамических функций отдельных реакций производится согласно закону Гесса. База данных программного комплекса обеспечивает хранение, просмотр и извлечение данных о термодинамических свойствах веществ и других данных, необходимых для расчетов. ПК TeDy нацелен на применение при проектировании и оптимизации широкого спектра технологических процессов, для проведения исследований и оценки поведения различных физико-химических систем. Программный комплекс успешно применялся при выполнении работ по моделированию ключевых технологических процессов замкнутого ядерного топливного цикла.
Бесплатно

Программный комплекс WAVES-l для моделирования и визуализации волновых процессов в упругом слое
Статья научная
Приводится описание программного комплекса WAVES-L, предназначенного для быстрого и наглядного моделирования волновых полей, возбуждаемых заданными поверхностными нагрузками в упругом слое с тремя вариантами закрепления его нижней грани. В основу алгоритмов, реализованных в программном комплексе, положены явные интегральные представления решения соответствующих краевых задач и выведенные из них асимптотические представления для бегущих волн в дальней от источника зоне. В ближней зоне, в том числе и под источником, амплитудно-частотные характеристики суммарного волнового поля определяются с помощью численного интегрирования несобственных контурных интегралов. Реализована также возможность расчета дисперсионных характеристик и собственных форм каждой из возбуждаемых бегущих волн. Удобный графический интерфейс позволяет легко менять параметры модели и визуализировать частотные и пространственные зависимости волновых характеристик как для отдельных волн, так и для суммарного волнового поля. Программный комплекс может применяться при решении задач ультразвукового неразрушающего контроля и мониторинга состояния тонкостенных инженерных конструкций, а также использоваться в учебном процессе в рамках лабораторных занятий по профильным предметам. Приводятся снимки экрана программы и графики сравнения дисперсионных кривых с экспериментальными измерениями.
Бесплатно

Проектирование новых конструкции тканевых бронепанелей с использованием суперкомпьютерных вычислений
Статья научная
Рассматривается проектирование новых конструкций тканевых бронепанелей с использованием суперкомпьютерных вычислений. Проведены численные эксперименты по исследованию масштабируемости задач динамического взаимодействия индентора с тканевыми защитными преградами, расположенными на регистрирующей среде с применением пакета программ LS-DYNA. Были даны рекомендации по увеличению защитных свойств бронепанелей.
Бесплатно

Разностная схема для решения уравнений роста опухоли с учетом ограничения потока
Статья научная
В статье исследована одномерная математическая модель роста раковой опухоли в квазилинейных уравениях параболического типа. В модели вводится ограничение на полный поток подвижных опухолевых клеток, что приводит к возможности вырождения системы уравнений в гиперболический тип и появлению разрывных (слабых) решений. Для нахождения слабых решений развитие опухоли трактовалось как появление новой фазы. В итоге решение задачи свелось к решению обобщенной (нелинейной) задачи Стефана. Предложена и реализована разностная схема для данной задачи с явным выделением подвижной границы фазового перехода. Показано, что данный подход позволяет описывать различные режимы опухолевого роста.
Бесплатно

Разработка и анализ быстрого псевдоспектрального метода решения нелинейных задач Дирихле
Статья научная
Разработан метод численного решения 1D, 2D и 3D краевых задач Дирихле для нелинейных уравнений эллиптического типа. Метод основан на применении чебышевских приближений искомой функции, не имеющих насыщения, и нового подхода к формированию и решению задач линейной алгебры при дискретизации исходных дифференциальных уравнений. При этом дифференциальные операторы аппроксимируются с помощью матриц, а само уравнение (в 2D и 3D случаях) - с помощью уравнения Сильвестра, либо его тензорного обобщения. В тестовых задачах с решениями различной степени гладкости показана зависимость порядка сходимости предложенного метода от гладкости искомого решения, строго соответствующая оценкам погрешности наилучших полиномиальных приближений. Указанные свойства свидетельствуют об отсутствии насыщения алгоритма и обеспечивают низкий расход памяти и машинного времени при численном анализе задач, решения которых имеют высокий порядок гладкости.
Бесплатно

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

Статья научная
Предложен, реализован и верифицирован новый высокоточный hp-вариант метода коллокации и наименьших квадратов (hp-МКНК) численного решения эллиптических задач в нерегулярных областях. При построении приближенного решения использовались граничные нерегулярные ячейки (н-ячейки), отсеченные границей области от ячеек прямоугольной сетки, и их законтурные части для записи уравнений коллокации и условий согласования. В малых и (или) вытянутых несамостоятельных н-ячейках отдельное решение не строилось, а продолжалось из соседних самостоятельных ячеек, в которых использовалась внешняя (и внутренняя в многосвязной области) часть границы области, заключенная в этих несамостоятельных н-ячейках, для записи краевых условий. Такой подход существенно упростил компьютерную реализацию разработанного hp-МКНК по сравнению с предыдущим хорошо зарекомендовавшим его вариантом, не потеряв при этом своей эффективности. Показана возможность уменьшения степени переопределения системы линейных алгебраических уравнений по сравнению с ее значениями в традиционных вариантах МКНК при решении бигармонического уравнения. Проведено сравнение с результатами других работ с демонстрацией преимуществ нового подхода. Приведены результаты расчетов кольцевых пластин различной толщины на изгиб в рамках теорий Кирхгофа - Лява и Рейсснера - Миндлина с помощью hp-МКНК, демонстрирующего отсутствие сдвигового запирания.
Бесплатно

Разработка параллельной СУБД на основе последовательной СУБД PostgreSQL с открытым исходным кодом
Статья научная
Статья посвящена разработке параллельной системы управления базами данных (СУБД) путем внедрения механизмов параллельной обработки запросов на основе концепции фрагментного параллелизма в свободно распространяемую на уровне исходных кодов СУБД PostgreSQL. Описана архитектура и принципы реализации параллельной СУБД PargreSQL для кластерных вычислительных систем, разрабатываемой на основе свободно распространяемой СУБД PostgreSQL. СУБД PostgreSQL является подсистемой в рамках системы PargreSQL. Описаны изменения, которые требуется внести в исходные тексты подсистем СУБД PostgreSQL. В исходные тексты PostgreSQL вносятся минимальные изменения. Изменения в структурах данных и алгоритмах инкапсулируются в новых файлах исходных текстов, подключаемых к исходным текстам PostgreSQL. Использование PargreSQL является прозрачным для пользовательских приложений. Подключение PargreSQL к прикладным программам, которые до этого использовали PostgreSQL, производится с минимальными изменениями в исходных кодах приложения. Параллельная СУБД PargreSQL, запущенная на одном вычислительном узле, работает так же, как последовательная СУБД PostgreSQL.
Бесплатно

Статья научная
В последнее время широкое распространение получили криптографические примитивы, которые базируются на эллиптических кривых над конечными полями. Основная причина этого заключается в том, что эллиптические кривые позволяют строить примеры конечных абелевых групп с хорошими, для криптографических целей параметрами. Кроме того, меняя характеристику поля можно легко повышать стойкость шифра. Существенную роль играет возможность удобной программной реализации. Разработан и программно реализован криптографический протокол цифровой подписи на основе эллиптических кривых. Протокол производит шифрование сообщения, формирование цифровой подписи, передачу сообщения и расшифровку на стороне получателя. Проанализирована криптографическая стойкость протокола несколькими методами. Построен график зависимости криптографической стойкости протокола от характеристики конечного поля, над которым строится эллиптическая кривая. Написана программа на языке C++ в среде программирования Visual C++ 2010 с поддержкой библиотеки больших чисел GMP, производящая шифрование и дешифрование сообщения в соответствии с построенным протоколом. Разработанная программа является инструментом, позволяющим передавать и получать сообщения с достаточной степенью криптографической стойкости и приемлемой скоростью.
Бесплатно

Распараллеливание рекуррентных циклов с предварительным вычислением суперпозиций
Статья научная
Как правило, именно циклы являются участками кода, вычисление которых занимает много времени. Поэтому, именно на них направляются особые усилия при ускорении программ, в частности, через распараллеливание. В статье описывается алгоритм распараллеливания циклов, вычисляющих элементы рекуррентно заданной последовательности. Рекуррентные циклы, рассматриваемые в статье, непосредственно распараллелены быть не могут. С помощью вспомогательных преобразований иногда их можно привести к циклам, допускающим параллельное выполнение. Ранее автором статьи был опубликован другой алгоритм распараллеливания циклов, вычисляющих элементы рекурсивно заданной последовательности. В современных процессорах время выполнения арифметических операций оказывается на порядок меньше, чем считывание аргументов этих операций из оперативной памяти. В данной статье приводятся оценки сложности по обращению к памяти, для описываемого алгоритма. Представленный в статье параллельный алгоритм оказывается более эффективным по обращениям к памяти, чем алгоритм, описанный автором ранее.
Бесплатно

Статья научная
В статье рассмотрен подход к масштабированию сервиса морфологического разбора слов естественного языка при обработке различных коллекций документов на русском языке. Выполнен обзор и критический анализ существующих решений. Сформированы требования к инструментальной среде словарного морфологического анализатора. Распределенная архитектура Web-сервиса морфологического анализа, предназначенного для обработки крупных коллекций документов на русском языке, представлена в виде структурной модели. Данная архитектура реализована в виде прототипа системы на языке программирования Ruby. Приведена структура используемого морфологического словаря в виде реляционной схемы. Испытания данного метода в распределенной вычислительной среде показали линейную масштабируемость предлагаемого решения. Конфигурация эксперимента включает систему генерации нагрузки в виде HTTP-запросов, систему балансировки нагрузки на рабочие узлы распределенной системы, серверы приложений с функционирующим анализатором и базу данных морфологического словаря, а также кэширующий узел для снижения издержек при выполнении запросов к словарю. Применение данного подхода позволяет получить линейный рост производительности в распределенных системах автоматической обработки больших объемов текста.
Бесплатно

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

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