Computational Science

Pytania i odpowiedzi dla naukowców używających komputerów do rozwiązywania problemów naukowych

1
Czysto obrotowe dopasowanie najmniejszych kwadratów
Czy ktoś mógłby polecić metodę dla następującego problemu najmniejszych kwadratów: znajdź który minimalizuje: , gdzie jest jednostką (obrót) matryca.N ∑ i = 0 ( R x i - b i ) 2 → minR ∈ R3 × 3R∈R3×3R \in \mathbb{R}^{3 \times 3}∑i = 0N.( R xja- bja)2)→ min∑i=0N(Rxi−bi)2→min\sum\limits_{i=0}^N (Rx_i - …
3
Pytanie o przybliżenie metodą najmniejszych kwadratów
Idę na kurs obliczeń naukowych i właśnie przekroczyliśmy przybliżenie metodą najmniejszych kwadratów. Moje pytanie dotyczy konkretnie przybliżenia za pomocą wielomianów. Rozumiem, że jeśli masz n + 1 punktów danych, możesz znaleźć unikalny wielomian stopnia n, który opisuje wszystkie te punkty. Ale rozumiem też, dlaczego nie zawsze jest to idealne. Przy …
2
Jak wykreślić powierzchnię wykresu 4D?
Próbuję wykreślić funkcję fali dla cząstki w polu 3D. Wymaga to ode mnie wykreślenia 4 zmiennych: osi x, y, z oraz funkcji gęstości prawdopodobieństwa. Funkcja gęstości prawdopodobieństwa to: abs((np.sin((p*np.pi*X)/a))*(np.sin((q*np.pi*Y)/b))*(np.sin((r*np.pi*Z)/c)))**2 Używam np.arange()dla X, Y i Z. Przeczytałem, że aby to zrobić, musisz narysować powierzchnię wykresu 4D. Oto jak powinno to wyglądać:
6
Czy istnieje implementacja BLAS na poziomie odniesienia w C / C ++?
Implementacja BLAS w netlib jest doskonałym źródłem informacji, ponieważ jest w większości niezoptymalizowana i dobrze udokumentowana (np. Zgemm ). Jest jednak w Fortran 77, co czyni go nieco niedostępnym dla osób z bardziej nowoczesnym wykształceniem programistycznym. Czy istnieje implementacja BLAS na poziomie odniesienia, taka jak netlib, w C / C …
11 blas  education 
1
Numeryczne metody odwracania transformat całkowych?
Próbuję odwrócić liczbowo następującą transformację całkowitą: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)d xF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Więc dla danego F(y)F(y)F(y) muszę aproksymować f( x )f(x)f(x) gdzie: fa( x )f(x)f(x) ifa(y)F(y)F(y) są rzeczywiste i dodatnie(są to ciągłe rozkłady prawdopodobieństwa) x , yx,yx,y są rzeczywiste i dodatnie(są wielkościami) W tej chwili mam bardzo niechlujną …
2
Solidna mechanika ze skończonymi różnicami: jak obsługiwać „węzły narożne”?
Mam pytanie dotyczące kodowania warunków brzegowych dla mechaniki bryły (elastyczność liniowa). W szczególnym przypadku muszę użyć różnic skończonych (3D). Jestem bardzo nowy w tym temacie, więc być może niektóre z poniższych pytań mogą być bardzo podstawowe. Aby doprowadzić do mojego konkretnego problemu, przede wszystkim chcę pokazać, co już zaimplementowałem (aby …
3
Czy mając trójridiagonalny układ liniowy SPD, możemy wstępnie obliczyć, aby dowolne trzy wskaźniki można było połączyć w czasie O (1)?
Rozważmy symetryczny dodatnio określona tridiagonal system liniowy , gdzie i . Biorąc pod uwagę trzy wskaźniki , jeśli przyjmiemy tylko rzędy równań ściśle między i hold, możemy wyeliminować zmienne pośrednie, aby uzyskać równanie w postaci gdzie . To równanie odnosi wartość do niezależnie od wpływu „zewnętrznego” (powiedzmy, jeśli wprowadzono ograniczenie …
2
Struktury danych dla skończonego kodu woluminu: tablice vs klasy
Muszę napisać skończony kod głośności dla magnetohydrodynamiki (MHD). Pisałem wcześniej kod numeryczny, ale nie w tej skali. Chciałem tylko zapytać, który będzie dobrym wyborem, używając struktury danych (podejście zorientowane obiektowo) z klasami lub po prostu używając wielu tablic dla różnych właściwości, pod względem szybkości, skalowalności itp. Planuję napisać kod w …
3
Siatkowanie danych powierzchni 3D w pythonie
Mam zestaw danych 3-wymiarowych punktów, dla których chciałbym zbudować siatkę, używając Pythona. Całe oprogramowanie, które widziałem, wymaga podania krawędzi. Czy w Pythonie jest program, który jako dane wejściowe przyjmuje zestaw punktów w 3D i generuje siatkę? Jeśli to możliwe, chciałbym, aby siatka była jednolita.
1
Czy ktoś używa metod szacowania oprogramowania w swoich badaniach w dziedzinie obliczeń?
W pracy zasadniczo działam jako niezależny konsultant. W przypadku kierownictwa i klientów muszę oszacować czas potrzebny na opracowanie oprogramowania w ramach moich badań w dziedzinie obliczeń. Jednak moje prognozy czasu są zwykle wyłączone. Wiem, że istnieją metody szacowania czasu potrzebnego na opracowanie oprogramowania. Czy metody te pozwolą mi dokładnie oszacować …
11 software 
3
Najlepsze metody zarządzania siatką w równoległych obliczeniach elementów skończonych?
Obecnie opracowuję metodę dekompozycji domen dla rozwiązania problemu rozpraszania. Zasadniczo rozwiązuję iteracyjnie system BVP firmy Helmholtz. Dyskretyzuję równania metodą elementów skończonych na siatkach trójkątnych lub czworościennych. Rozwijam kod w kierunku mojej pracy doktorskiej. Zdaję sobie sprawę z niektórych istniejących bibliotek elementów skończonych, takich jak deal.ii lub DUNE i chociaż uważam, …
3
Jak należy traktować niestałe współczynniki za pomocą schematu podmuchu wiatru o skończonej objętości pierwszego rzędu?
Począwszy od równania doradczego w formie konserwatorskiej. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x gdzie a(x)a(x)a(x) jest prędkością zależną od przestrzeni, uuu jest stężeniem gatunku, który jest zachowany. Dyskretyzacja strumienia (gdzie strumień f=a(x)uf=a(x)uf=a(x)u jest zdefiniowany na krawędziach komórek między punktami siatki) daje, ut=1h(fjot−12)-fjot+12))ut=1h(fajot-12)-fajot+12)) u_t = \frac{1}{h}\left( f_{j-{\frac{1}{2}}} - f_{j+{\frac{1}{2}}} \right) Stosując podmuch wiatru …
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.