«quadrature» 태그된 질문

나는 1967 년 Bulirsch & Stoer (Numerische Mathematik, 9, 271-278)가 발표 한 ALGOL 프로그램의 C ++ 적응 인 수치 적분 (quadrature)을위한 작은 서브 루틴을 가지고 있습니다. 더 현대적인 (적응) 알고리즘으로 업그레이드하고 그러한 (무료) C ++ 라이브러리가 있는지 궁금합니다. 나는 GSL (C)로 보였지만 끔찍한 API가 제공됩니다 (숫자가 좋을 수도 있지만). 다른 …

10 c++  quadrature 

하자 삼각형하고하자 에 원활하게 기능 할 수 .TTTfffTTT 중간 점 구적법 . 여기서 은 의 중간 점입니다 .∫fdx≈|T|⋅f(xM)∫fdx≈|T|⋅f(xM)\int f dx \approx |T|\cdot f(x_M)xMxMx_MTTT 심플 렉스에 고차 공식을 제공 할 수 있습니까?

10 quadrature 

제목에서 알 수 있듯이 삼각형에서 컴팩트하게 지원되는 함수 (Wendland의 quintic polynomial)의 적분을 계산하려고합니다. 함수의 중심은 3 차원 공간 어딘가에 있습니다. I는 임의의 기능이 통합되지만 작은 삼각형 ( R은 전자 &lt; ( R은 D I U S / 4 ) (2)area&lt;(radius/4)22area&lt;(radius/4)22area < \frac{(radius/4)^2}{2} ). 현재 Dunavant, 1985 (p = 19)에서 설명한 …

10 quadrature 

한동안 업데이트되지 않았고 액세스가 제한되어있는 직교 규칙의 백과 사전 이 하나 이상 있습니다. 이 출처는 몇 가지 고전적이고 현대적인 출처를 말하며 일반적으로 잘 정리되어 있습니다. 그러나 순전히 이론적 인 접근 방식에서 직교 규칙의 구성에 접근하므로 유한 요소 계산을위한보다 실용적인 방법을 놓칩니다. 직교 법 규칙에 대한 다 분야 학문적 개요가 있습니까? …

10 libraries  finite-element  quadrature 

진동 적분에 대한 대부분의 방법이 나는 형태의 적분과 거래에 대해 알고 IS 큽니다.∫f(x)eiωxdx∫f(x)eiωxdx \int f(x)e^{i\omega x}\,dx ωω\omega I는 형태의 중요한있는 경우 뿌리 만 약 알려진 진동 함수가 있지만 점근선 형태의 일종 는 주파수가 완전히 다른 (및 선형 적으로 독립적 임) 알려진 것으로, 이 적분을 어떻게 평가할 수 있습니까?∫f(x)g1(x)⋯gn(x)dx,∫f(x)g1(x)⋯gn(x)dx, \int f(x)g_1(x)\cdots …

9 quadrature 

예를 들어, 일부 도메인에서 의 -norm 을 수치 적으로 계산하고 싶습니다. 0을 포함하고 Gauss quadrature를 시도했지만 실패했습니다. 구 좌표를 사용하여 단위 공 의 실제 -norm과는 거리가 멀었습니다.이 작업을 수행하는 좋은 방법이 있습니까? 이 문제는 종종 재진입 코너가있는 도메인의 유한 요소 컴퓨팅 장난감 문제에서 나타납니다. 감사.엘2엘2L^2u =1(엑스2+와이2+지2)1 / 3유=1(엑스2+와이2+지2)1/삼\displaystyle u = …

9 finite-element  quadrature 

통합하려고합니다 ∫10t2n+2exp(αr0t)dt∫01t2n+2exp⁡(αr0t)dt\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt 이것은 간단한 변환입니다 ∫∞1x2nexp(−αr0x)dx∫1∞x2nexp⁡(−αr0x)dx\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx 사용 t=1xt=1xt = \frac1{x}부적절한 적분을 수치 적으로 근사하기가 어렵 기 때문입니다. 그러나 이로 인해 새로운 정수를 거의 0으로 평가하는 문제가 발생합니다. 간격이 길이가 1이므로 적절한 구적 노드 수를 얻는 것이 매우 쉽습니다 (따라서 비슷한dtdtdt 아주 작게 만들 수 있습니다), …

9 quadrature  accuracy