10.01.09
Птица
Магнитный вихрь, если его сверху стукнуть полем, ведет себя как птица
Прочтите эту запись до конца »
08.08.09
Физика-лирика
Ферромагнитных свойств ясна причина —
Непарный электрон в них виноват:
Все атомы по направленью спина,
Глдядящего вперед или назад,
Построены, как войско на парад.
Во внешнем поле, разрушая стены,
Сливаются соседние домены!
Так создает гармонию Вселенной
Ничтожных сил суммарный результат.
Дж. Апдайк, Танцы твердых тел.
Ж-л «Химия и жизнь», № 11, 1969.
06.21.09
Ландавшиц
Из предисловия к четвертому изданию третьего тома (Квантовая механика) в 1989 г.:
Феноменальное владение Львом Давидовичем аппаратом теоретической физики позволяло ему сплошь и рядом обходиться без обращения к оригинальным работам для воспроизведения тех или иных результатов своим путем.
Е. М. Лившиц
Будем честными, с подобным феноменальным владением примерно половину курса можно было бы заменить фразой «очевидно что»
06.05.09
Система управления заданиями PBS
Часто для управления выполнением пользовательских задач на кластере используют систему PBS. Она очень проста в использовании, но разрозненные в сети сведения о том, как нужно с ней работать, хотелось бы собрать в одном месте, что я и попытаюсь ниже сделать.
Наиболее часто используемыми программами из системы PBS являются qsub, qdel и qstat. Первая ставит задание на очередь, вторая удаляет его оттуда, а третья позволяет следить за состоянием задания. Прочтите эту запись до конца »
05.27.09
Использование OpenMP: старт
Распараллеливание задач, подчас, бывает крайне необходимой штукой. Однако в самом начале изучения одной из библиотек для этого — OpenMP — сталкиваешься со страшной проблемой: нифига не паралеллится. Итак, раз и на всегда что нужно делать:
#include <stdio.h>
#include <stdlib.h>
#include <omp.h>
int main(int argc, char ** argv)
{
#pragma omp parallel
{
printf("Hello, world!\n");
}
return 0;
}
Теперь компилируем:
alexander@A64x2 ~/Temp $ gcc -fopenmp -c test.c
alexander@A64x2 ~/Temp $ gcc -o test test.o -lgomp
Здесь ключевыми являются флаги -fopenmp и -lgomp, использование которых, похоже, очевидно всему миру кроме меня… А результат выполнения у меня был таков:
alexander@A64x2 ~/Temp $ ./test
Hello, world!
Hello, world!
05.08.09
Метод Рунге-Кутта с адаптивным шагом (Runge-Kutta-Fehlberg)
Метод интегрирования Рунге-Кутта с адаптивным шагом является одним из наиболее часто используемых методов численного решения дифференциальных уравнений. Его реализаций существует довольно много, но наибольшее распространение получил метод Рунге-Кутта четвертого порядка. Для достижения заданной точности применяется изменение шага интегрирования. Можно изменять по-простому, деля шаг на два, и интегрируя дифференциальне уравнение сразу с двумя шагами по сетке, а можно поступить немного иначе: интегрировать дважды с одним шагом, но разным порядком метода. Получается быстрее. Далее приведены три матлабловских скрипта, решающие уравнение
с помощью интегрирования 4-м и 5-м порядками метода. Прочтите эту запись до конца »
02.18.09
wxplot2d & wxplot2d в wxmaxima
Если не работают функции wxplot2d и wxplot3d в wxmaxima, следует перекомпилировать gnuplot с флагом gd
11.25.08
МНК и его погрешность
В эксперименте были намеряны значения , которым соответствуют
. После получения апроксимирующего полинома (в нашем случае прямой
), нужно узнать погрешности коэффициентов и самой аппроксимации. Прочтите эту запись до конца »
11.17.08
Решение нелинейных минимаксных задач с помощью r-алгоритмов
Цитата из книги Н. З. Шора «Методы минимизации недифференцируемых функций и их приложения», К., «Наукова думка» 1979.
Глава4, §5. Использование r-алгоритмов для решения нелинейных минимаксных задач.
Рассмотрим вопросы минимизации функций вида , где
— гладкие функции,
. Если
— выпуклые функции,
, то в качестве обобщенного градиента можно брать субградиент
, где
— индекс функции
для которой
.
В связи с важностью минимаксных задач большое внимание в настоящее время уделяется разработке методов их решения. Ряд таких методов для определенных классов задач предложен в работах В. М. Демьянова, Б. Н. Пшеничного, Е. Г. Евтушенко и др.
В принципе задача отыскания минимакса может быть сведена к решению следующей задачи выпуклого программирования: найти при ограничениях
,
. Прочтите эту запись до конца »
10.15.08
Таблицы спектров атомов
NIST Atomic Spectra Database
По всем элементам, интервалам частот и еще многим не известным мне параметрам.