Wann ist es vorteilhaft, Integrale numerisch zu iterieren?

8

Wenn es ein -dimensionales Integral der Form [ 0 , 1 ] n + 1 f ( x , y ) gibt(n+1) normalerweise würde man dies unter Verwendung einer mehrdimensionalen Integrationsbibliothek über die gesamte Domäne [ 0 , 1 ] n + 1 bewerten.

[0,1]n+1f(x,y)dnxdy,
[0,1]n+1

Aber gibt es einige Bedingungen, unter denen es sinnvoll sein könnte, das Integral über separat unter Verwendung einer eindimensionalen Quadratur auszuführen und dann die mehrdimensionale Integrationsbibliothek zu verwenden, um den Integranden über die anderen n Koordinaten auszuwerten ? [ 0 , 1 ] n g ( x )yn

[0,1]ng(x)dnx,g(x)=01f(x,y)dy.

Dies kann beispielsweise sinnvoll sein, wenn als Funktion von y besonders glatt ist, nicht jedoch x . Aber wie glatt muss es in diesem Fall genau sein? Meine Vermutung war, dass es fast nie Sinn macht, weil viel zu viele der 1-D-Quadratur-Bewertungspunkte "verschwendet" würden, aber ich bin mir nicht sicher, ob dies immer zutrifft. Wird dies durch das Design der hochdimensionalen Integrationsmethoden garantiert?fyx

In meinem Fall ist Black-Box, aber stückweise glatt in y , und hat eine unbekannte Menge Knicke und Sprünge in x an unbekannten Orten, und n ist ziemlich hoch ( n 4 ), so dass das Integral über x hat mit etwas speziell für viele Dimensionen gemacht werden. Das Integral über y kann mit etwas Regelmäßigem wie gemacht werden . In diesem Beispiel ist die Funktion in y so flüssig, dass sie fast zu funktionieren scheint, aber die wiederholte Integration ist 30-mal langsamer, sodass ich mich frage, ob der Ansatz falsch ist.fyxnn4xyquadgky

Wenn Sie wissen, wo dies bereits in der Literatur diskutiert wird, wäre dies ebenfalls hilfreich.

[0,1]nex1x2xndnx=F({1,,1}n{2,,2}n|1).
n(n1)x1g(x2:n)=(ea1)/aa=x2xn

(n=5)N0.00244N10.00167N14g1.5

1.51.5g=(1ea)/a1.5

xy

Kirill
quelle
2
fy
2
y(n+1)
y1n(n+1)
Ah, tut mir leid. Ich drehte es in meinem Kopf herum. Ich kann nicht sagen, dass ich hier eine gute Faustregel kenne, also muss jemand anderes diese nehmen.
Tyler Olsen

Antworten:

4

Klarstellung: Meine Antwort wurde speziell für adaptive Integrationsroutinen mit deterministischer Fehlerkontrolle wie dieser geschrieben . Es wird für spärliche Gitter- und Monte-Carlo-basierte Integrationsroutinen, deren Fehlerkontrolle nicht auf die unten beschriebene Weise durchgeführt wird, umstritten.

Ein erheblicher Kostenfaktor für automatische Tensorprodukt-basierte Black-Box-Integrationsroutinen ist die Fehlerkontrolle unter zwei Gesichtspunkten

  1. Verschwendete Funktionsbewertungen. Alle adaptiven Integrationen funktionieren, indem der Integrand und der Fehler unter Verwendung einer Regel niedriger Ordnung oder einer gröberen Partitionierung geschätzt und mit Regeln höherer Ordnung oder feineren Partitionen wiederholt werden, bis die Fehleranforderungen erfüllt sind. Durch verschachtelte Integrationsregeln können einige der in früheren Schritten geleisteten Arbeiten recycelt werden, häufig jedoch nicht alle.
  2. Um Funktionsbewertungen zu erhalten, werden in adaptiven Integrationscodes häufig stark verschachtelte Regeln wie Gauss-Kronrod oder Newton-Cotes verwendet. Verschachtelte Quadraturregeln sind suboptimale Quadraturregeln, da sie für eine bestimmte Funktionsklasse über eine feste Quadraturordnung erheblich ungenauer sind als optimale Regeln (z. B. Gauss-Legendre und Clenshaw Curtis). Mit anderen Worten, verschachtelte Regeln nutzen Funktionsbewertungen weniger effizient.

yyf(x,y)yr

|k=1rwkf(x,yk)[0,1]f(x,y)dy|ϵfor all x[0,1]n,
f(x,y)yxrxϵ

g(x)ng(x)

f(x,y)

Um ein Anwendungsbeispiel zu nennen: Dieses genaue Problem trat für mich bei der Bewertung von Volumen-zu-Volumen-Singularintegralen in diesem Artikel auf , und meine Behandlung ähnelt der oben vorgeschlagenen. Als Faustregel ist es immer ratsam, so viele Dimensionen wie möglich mithilfe analytischer Argumente zu entfernen, bevor das Problem durch eine Black-Box-Integrationsroutine geleitet wird.

Richard Zhang
quelle
y
Mein Argument ist einfach, dass, wenn Sie eine adaptive Methode mit deterministischer Fehlerkontrolle annehmen, die Auswertung einiger Dimensionen in geschlossener Form (oder halb geschlossener Form) unvermeidliche Schritte eliminiert, die sonst numerisch ausgeführt würden. Ihr (ausgezeichnetes) Beispiel ist jedoch eines, bei dem eine deterministische adaptive Methode überhaupt nicht angewendet werden würde.
Richard Zhang
Angenommen, Sie haben eine deterministische Standardanpassungsmethode verwendet, z. B. ab-initio.mit.edu/wiki/index.php/Cubature . Dann wäre ich sehr, sehr überrascht, wenn Sie durch Ausschneiden nicht eine Beschleunigung des gesamten Faktors erhalten eine der Dimensionen semi-analytisch.
Richard Zhang
Wenn Sie Cubature genau so verwenden, wie Sie es vorschlagen (was ein durchaus vernünftiger Vorschlag ist), habe ich in erster Linie die "30-mal langsamere" Zahl in meiner Frage erhalten, also war ich überrascht (daher die Frage). Ich habe das Monte-Carlo-Beispiel nur als etwas leicht zu analysierendes gemeint, ich interessiere mich eigentlich mehr für deterministische Methoden.
Kirill
n+1(n+1)