Обчислювальна наука

Тому я обговорював, чи варто мені турбуватися над вивченням Python. Якщо говорити з моїми професорами, Матлаб, здається, є загальною мовою, що використовується у прикладній математиці / обчислювальній науці, що стосується наукових шкіл; тоді як у галузі , мої професори (особливо ті, хто працював у галузі) сказали, що навчання с ++ …

9 matlab  python  c++ 

Я використовую панель інструментів Matlab pde для вирішення певного еліптичного рівняння в 2D. Розв’язання чудово, хоча мені потрібно побудувати його по заданій лінії, тобто вирізати плоский зріз з 3D-сітки, що представляє рішення. Я не можу з’ясувати спосіб, який інтелектуально включає функції інструментарію (тобто не залучає інтерполяцію низького рівня на трикутну …

9 pde  matlab 

Для своєї дослідницької роботи я розробив невелику бібліотеку C ++, спрямовану на полегшення зв'язку між обчислювальними кодами C ++ та Octave / Matlab (коли останні використовуються для цілей після обробки). Я хотів би випустити таку бібліотеку під GPL на деяких безкоштовних веб-сайтах, таких як Sourceforge або Github. Чесно кажучи, мені …

9 matlab  libraries  c++ 

Я хочу вирішити Ku=bKu=bK u = b де KKKмоя матриця жорсткості. Однак деякі обмеження можуть бути відсутніми, тому деякий жорсткий рух тіла може залишатися в системі (через власне значення нуля). Оскільки я використовую CG для вирішення лінійної системи, це не прийнятно, оскільки іноді CG не сходить на напівпозитивних проблемах (але …

9 finite-element  iterative-method  conjugate-gradient 

З визначення числа умови виходить, що для його обчислення потрібна інверсія матриці, мені цікаво, чи для загальної квадратної матриці (або краще, якщо симетричний позитивний певний) можливо використовувати деяку декомпозицію матриці для обчислення числа умови в швидший шлях.

9 linear-algebra  matrix  condition-number 

Чи є кращий спосіб, як здійснити швидку (приблизну) оцінку інтерполяційного полінома Чебишева на рівномірній сітці (з урахуванням значень функції у вузлах Чебишева)? Моя проблема полягає в тому, що інтерполяція стає повільною, коли ступінь інтерполяційного полінома зростає. На думку мені прийшли такі ідеї: Спробуйте адаптувати нерівномірні методи FFT (NFFT) Використовуйте FFT …

9 interpolation  fourier-analysis  fftw 

Я намагаюся реалізувати метод Бройдена-Флетчера-Голдфарба-Шенно, щоб знайти мінімум функції. Мені потрібні два початкові здогадки & та початкове наближення матриці . Єдині вимоги, які я знаходжу до - це те, що якщо гессіан є симетричним позитивним певним, то також повинен бути . Дивлячись на wikipedia, я бачу, що типовим початковим наближенням …

9 optimization  algorithms 

Проведення прямого чисельного моделювання ізотропної турбулентності псевдоспектральним методом Фур'є (Orzag & Patterson, PRL, 1972) за допомогою FFT. Для ознайомлення з методом, який широко використовується в турбулентній спільноті, ви можете ознайомитись з цим курсом: http://www.math.ualberta.ca/~bowman/m655/lab3d.pdf Використовуючи т. Зв 2 / 32/32/3 правила дезаліювання, які складаються з виконання у^( к , т …

9 fluid-dynamics  simulation 

У мене є проблема, коли мені потрібно знайти всі позитивні (як у власній значенні позитивні) власні пари невеликої (зазвичай менше 60x60) несиметричної матриці. Я можу припинити обчислення, коли власне значення менше певного порогу. Я знаю, що власні значення справжні. Будь-які пропозиції щодо алгоритмів я міг би використати, щоб спробувати витіснити …

9 performance  eigensystem  gpu 

Більшість моїх програмувань - це одноразові дослідницькі коди на C для мого власного використання. Я ніколи не поширював жодного коду, окрім близьких співробітників. Я розробив алгоритм, який публікую у науковому журналі. Я хочу надати вихідний код та, можливо, виконуваний код в онлайн-додатку до статті. Колега попросив мене зробити узагальнення до …

9 linear-solver  libraries  c++  c 

Я хочу мінімізувати складну цільову функцію, і я не впевнений, чи вона опукла. Чи є приємний алгоритм, який намагається довести, що він не опуклий? Звичайно, алгоритм міг цього не довести, і в такому випадку я б не знав, випуклий він чи ні, і це нормально; Я просто хочу спробувати виключити …

9 convex-optimization 

У мене складається враження, що з дуже різних ресурсів та переговорів з дослідженнями зростає попит на високоточні обчислення в числових часткових диференціальних рівняннях. Тут висока точність означає більшу точність, ніж просто стандартна 64-бітна подвійна точність. Цікаво про сучасність цієї теми. Для порівняння, в числовому PDE є комуніти, які конкретно націлені, …

9 reference-request  precision 

Мені хотілося б знати, чи існує алгоритм, який задає набір точок o, а кут обчислює опуклий корпус, якщо кут дорівнює α=0α=0\alpha = 0 і дав α>0α>0\alpha > 0 обчислює конверт, який більш уважно слідує "периметру". І якщо є визначення непересічного периметра набору точок, в цьому випадку отриманий многокутник, коли αα\alpha …

9 algorithms  computational-geometry 

Відповідь дає програмне забезпечення для обчислення умовного показника Ляпунова (CLE) для зв'язаних осциляторів в синхронізації хаосу. Однак, це важко прослідкувати, і графічного виводу сюжету немає (і він є в С, більш складний). Хтось знає, як змінити панель інструментів LET, що дуже добре для нерозв’язаної системи, але я не розумію, як …

9 matlab  ode  simulation  quantum-mechanics 

Переглядаючи кілька підручників, я помітив, що проблема спочатку обмеження дужки під час пошуку рядків має тенденцію (принаймні, в моїх студентських текстах). Чи є налагоджені методи чи найкраща практика для подібного типу проблем, чи рішення зазвичай залежать від додатків? Хтось може порекомендувати деякі посилання на тему?

9 optimization  reference-request