Ciencia computacional

Estoy realizando una búsqueda de línea como parte de un algoritmo cuasi-Newton BFGS. En un paso de la búsqueda de línea, uso una interpolación cúbica para acercarme al minimizador local. Sea ser la función de interés. Quiero encontrar una tal que .x ∗ f ′ ( x ∗ ) ≈ …

9 optimization  numerical-analysis  interpolation 

En el informe técnico sobre Galahad [1], los autores afirman, en el contexto de problemas generales de programación no lineal, En nuestra opinión, nunca había habido muchas dudas de que los métodos SQP [programación cuadrática secuencial] tendrían más éxito [que los métodos lagrangianos aumentados] a largo plazo ... ¿Cuál podría …

9 nonlinear-programming 

Muchas bibliotecas FEM adaptativas utilizan estructuras de datos de malla más avanzadas para manejar la adición / eliminación de nodos, bordes, triángulos, tetraedros, etc. Por ejemplo, la biblioteca p4est usa estructuras de datos octree para el refinamiento de malla adaptable; a menudo no encontraría octrees utilizados para los cálculos en …

9 linear-algebra  finite-element  sparse-matrix 

Estoy buscando una solución explícita rápida (¿me atrevo a decir que es óptima?) El problema real lineal de 3x3, , . Ax=bAx=b\mathbf{A}\mathbf{x} = \mathbf{b}A∈R3×3,b∈R3A∈R3×3,b∈R3\mathbf{A} \in \mathbf{R}^{3 \times 3}, \mathbf{b} \in \mathbf{R}^{3} La matriz es general, pero está cerca de la matriz de identidad con un número de condición cercano a …

9 linear-solver  reference-request 

Siempre he escuchado que la paralelización fácil era una de las ventajas de los métodos DG, pero realmente no veo por qué alguna de esas razones no se aplica también a Galerkin continuo.

9 finite-element  parallel-computing  discontinuous-galerkin 

Tengo un gran problema de valor propio cúbico: (A0+λA1+λ2A2+λ3A3)x=0.(UNA0 0+λUNA1+λ2UNA2+λ3UNA3)X=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Podría resolver esto convirtiendo a un problema de valor propio lineal, pero daría como resultado un sistema tan grande:32323^2 ⎡⎣⎢- A0 00 00 00 0yo0 00 00 0yo⎤⎦⎥⎡⎣⎢Xyz⎤⎦⎥= λ ⎡⎣⎢UNA1yo0 0UNA20 0yoUNA30 …

9 c++  eigensystem  eigenvalues 

Tengo un proyecto donde necesito usar campos cuadráticos Específicamente números de la forma a + b - 3---√una+si-3a + b \sqrt{-3} cona , b ∈ Quna,si∈Qa,b \in \mathbb{Q}. Por ejemplo, aquí están los números primos en enteros de Eisenstein : No quiero usar salvia. Me gustaría escribir mi propio tipo …

9 linear-algebra  precision 

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Lo anterior es mi código C …

9 c++  performance 

Estas preguntas pueden estar un poco fuera de tema en comp-sci. si es necesario, sugiera dónde encaja. La pregunta es sobre cómo administrar todas las ejecuciones de simulación de manera eficiente. digamos, por ejemplo, una simulación requiere la fijación de 2 parámetros que deben definirse en un cierto rango de …

9 monte-carlo  data-management 

¿Cuáles son las principales diferencias entre FEM y XFEM? ¿Cuándo deberíamos (no) usar XFEM intead de FEM y viceversa? En otras palabras, cuando encuentro un nuevo problema, ¿cómo puedo saber cuál de ellos debo usar?

9 finite-element  numerical-analysis  adaptive-mesh-refinement  numerics 

¿Hay alguna regla general simple para decir si vale la pena hacer SOR en lugar de Gauss-Seidel? (y posible forma de estimar el parámetro de realxation )ωω\omega Me refiero solo a mirar la matriz , o al conocimiento de un problema particular que representa la matriz. Estaba leyendo la respuesta …

9 linear-solver  iterative-method  heuristics 

Me pregunto si hay alguna lista de correo muy buena o grupos de Google para Computational Science, donde podamos discutir preguntas en lugar de solo hacer y responder preguntas. De hecho, estoy más interesado en la computación paralela y las soluciones numéricas de PDE. Pero no sé qué y cómo …

9 parallel-computing  reference-request 

He estado tratando de depurar este error los últimos días, me preguntaba si alguien tiene consejos sobre cómo proceder. Estoy resolviendo la ecuación de Poisson para una distribución de carga escalonada (un problema común en la electrostática / física de semiconductores) en una malla de volumen finito no uniforme donde …

9 discretization  finite-volume  poisson 

Considere el problema de optimización estrictamente convexo sin restriccionesDeje que x_ \ text {opt} denote sus mínimos únicos y x_0 sea ​​una aproximación inicial dada a x_ \ text {opt}. Llamaremos a un vector x una \ epsilon- solución cercana de \ mathcal {O} if \ begin {ecation} \ frac …

9 convex-optimization  complexity 

Estoy escribiendo un pequeño paquete MATLAB que resolverá ciertos problemas numéricos de clase. Hay 3 etapas del algoritmo y el usuario tiene 5 opciones para cada etapa. He implementado todo el problema usando funciones y 3 casos de cambio (uno para cada etapa de algoritmo). Funciona bien, pero estoy considerando …

9 matlab