вопрос по sparse solver

[Борман 2 375]

В разборки. Такие, какие Ньютон прописал, теормех попистайте. И когда только безграмотные дебилы закончатся?

Ваше желание оказалось пророческим.

[Борман 2 375]

Зачем оскорбляешь?

Т.е. амплитуда груза, колеблющегося на пружине не равна F/(N-Mw^2)?

У него по пьяни бывает. Отправлен на спецолимпиаду. Вернется не ранее чем через месяц.

[Цветочек 15]

Я бы не отказалась от титула "Знаток", если это понравится Модератору

[Борман 2 375]

Для этого решателя можно вывести расширенный протокол. Смотрите bcsoption.

[Борман 2 375]

Вот тест.

ut_fun_printf ("1") ;

double a = 0. ;

double b = 1. ;

for (int j = 0 ; j < 8 ; j++) {

for (i = 1 ; i <= 1000000000 ; i++) {

a += b ;

}

}

sprintf (str, "a = %e", a) ;

ut_fun_printf (str) ;

Получаем

a = 8.000000e+009

На компе 3.6 Ггц он идет 8 сек.

А вы утверждаете, что на компе 2.8 Ггц он идет 1 сек.

Чтобы ускорить эту операцию необходимо:

1. Если все пользовательские процессы являются дочерними процесса explorer.exe то необходимо измнить класс приоритета этого процесса на 32 (real time).

2. Повысить приоритет процесса, вызывающий поток в котором запущен цикл, на наивысший 32.

3. Изменить класс приоритета потока с циклом на Critical Priority.

Так что все возможно. Если на одноядерной тачке считает Ансис, то он вешает ее полностью. Не исключено, что Ансис повышает поиоритеты нужных процессов до очень высоких.

Кстати, можно попробовать через WinAPI повысить приоритет Ансиса до максимума, и посмотреть что получится.

[Борман 2 375]

будет иметь 1,2 миллиарда ненулевых элементов.

В то же время Ansys в процессе решения этой же матрицы на ж.д. в файл LN09 пишет всего 9ГБ.

Для хранения sparse-матрицы надо хранить сами значения (double - 8 байт) и целый индекс (long, ну пусть 4 байта). Итого:

1.2e9 * 12 / 1024^3 = 13.4 Гб. Неплохо сходится.

[Борман 2 375]

Алгоритм упоряд???вания также управляется bcsoption

Посмотрите здесь <noindex>http://fsapr2000.ru/index.php?act=Attach&a...st&id=42802</noindex> глава 7 - первое недокументированное поле этой команды.

Само рисование матриц уже реализовано в последних версиях...

[KostyaPC 0]

ISPA, <noindex>http://en.wikipedia.org/wiki/SIMD</noindex> может объяснить, почему количество GFLOPS не равно частоте процессора Изменено 6 марта 2012 пользователем KostyaPC

[KostyaPC 0]

Объясните.

Какое это имеет отношение к полной факторизации разреженных матриц.

Особенно комплксных.

отношение это имеет такое, что если проводятся однотипные операции сложения/умножения/..., то использование расширенного набора комнад процессора позволяет выполнять их параллельно. ansys использует <noindex>http://software.intel.com/en-us/articles/intel-mkl/</noindex> . соответственно там все максимально оптимизировано для работы с конкретным процессором, критические секции кода написаны на ассемблере и по максимуму используют расширенный набор команд.

Какие команды использовать?

читать мануалы к процессорам, писать с использованием ассемблерных вставок

[KostyaPC 0]

ISPA, вы впервые услышали про расширенный набор инструкций и утверждаете, что это необоснованная реклама. удачи вам, ispa.

[Борман 2 375]

Что же Вы в такую чепуху верите.

<noindex>http://en.wikipedia.org/wiki/Instructions_per_cycle</noindex>

[Борман 2 375]

ИСПА - лучшая программа всех времен и народов.

Все лапшеухие студенты пользуются ансис, а все профессионалы - ИСПА.

Все остальные численные методы должны быть верифицированы с обязательным сравнением с ИСПА.

УРА ИСПА! Этой эталонной программе для КЭМ расчетов.

Честь и хвала ее создателям!

Суть передана правильно.

[Борман 2 375]

ИСПА, какие еще мифы об АНСИС вы знаете ?

suvorov.as

В САЕ есть ряд тем, обсуждение которых обязательно заканчивается выводом, который вы и озвучили. Об объективности этого вывода речь, конечно, не идет.

[Борман 2 375]

А я словами Леонида Соловьева ("Повесть о Ходже Насреддине"):

- Тот, кто носит медный щит, тот имеет медный лоб. На месте сокола сидит филин. О джины, вы ищете там, где не прятали, поцелуйте за это под хвост моего ишака!

[Борман 2 375]

Когда АНСИС пишет скорость 8 Гфл это рекламный ход.

Так и вы напишите, ИСПА.

[soklakov 1 475]

Вот народ то на форуме над этим спором смеется наверное.

Вот уж да уж.

[a_schelyaev 366]

Я позволю себе постебаться немного

Осторожнее, он это крайне болезненно воспринимает.

[Цветочек 15]

какой у вас тут накал, я то думала уже потухло

[Цветочек 15]

Я бы не смогла переписывать то, что не писала

[Цветочек 15]

Я пишу о себе, ибо ничего ни о ком не знаю, и не осуждаю никого и ничего.