Al evaluar los armónicos cilíndricos, uno debe evaluar las funciones trigonométricas y , potencialmente para enteros grandes y . ¿Cuál es la mejor manera de hacer esto en código C? Actualmente, solo evalúo en el ángulo , pero sospecharía que las bibliotecas estándar pierden precisión con argumentos...
Tengo el siguiente programa lineal: dondex∈Rn,1Txdenota la suma de las entradas dex, yaes conocido y tiene entradas distintas estrictamente positivas.MaximizarSujeto aunaTXXmin≤ x ≤ xmax1Tx = 1MaximizeaTxSubject toxmin≤x≤xmax1Tx=1 \begin{array}{cc} \text{Maximize} & a^T x \\ \text{Subject to} &...
Tengo necesidad de calcular las funciones: y donde y y a menudo es muy pequeño ( ). ¿Hay alguna forma general de generar algoritmos altamente precisos para funciones "especiales" como estas?F( x ) = pecado- 1XXF(X)=pecado-1XX f(x) = \frac{\sin^{-1}x}{x}sol( x ) = pecadoa
Me gustaría saber cuál es más ventajoso cuando se trata de resolver ecuaciones hiperbólicas no lineales, elementos finitos o métodos de diferencias finitas. ¿Qué método será mejor para capturar choques? ¿Es posible proporcionar una respuesta detallada / referencias? Además, quiero resolver...
Pronto implementaré el método nodal discontinuo de Galerkin, y habiendo hecho esto antes de conocer las matrices de indexación básicas que tendré que calcular, dada una malla y datos polinomiales. El problema que encontré en el código anterior fueron errores sutiles que cometí al calcular cosas...
Necesito implementar una inversión de matriz densa en paralelo para un lenguaje que estoy usando que parece no tener una biblioteca existente para esto (específicamente IDL usando IDL Bridge para pasar mensajes). Estoy familiarizado con los métodos de programación paralela a través de la...
Me he encontrado con una ecuación rígida para resolver el problema circular restringido de tres cuerpos. [Un objeto se mueve considerando el efecto de las fuerzas gravitacionales causadas por dos fuentes gravitacionales fijadas en un espacio 2D]. Las ecuaciones son
Necesito aprender y utilizar un paquete de solución automatizada de volumen finito para un proyecto en el que estoy trabajando y lo he reducido a estos dos paquetes. Me preguntaba si alguien tiene experiencia con ambos paquetes y quizás podría comentar una comparación. Considerando el objetivo...
Sea un proceso Ito d X t = a ( X t , t ) d t + b ( X t , t ) d W t donde W t es un proceso Wiener.XtXtX_treXt= a ( Xt, t ) dt + b ( Xt, t ) dWtreXt=una(Xt,t)ret+si(Xt,t)reWt dX_t=a(X_t,t)dt + b(X_t,t)dW_t WtWtW_t Milstein propone una aproximación numérica de la solución de estas...
Estoy leyendo un libro sobre métodos numéricos y el cuadrado de la norma discreta L2L2L^2 se define como ||x||22=h∑1Nx2i||x||22=h∑1Nxi2||x||^2_2=h\sum_1^Nx^2_i Cada punto obtiene un "peso", que es hhh , por lo tanto, es como un promedio sobre los cuadrados de los valores en todos los puntos. De...
Supongamos que tengo el sistema parabólico con condiciones de contorno de Dirichlet y condición inicial u = g ,ut=∇⋅(k(x)∇u)+f,(x,t)∈Ω×Iut=∇⋅(k(x)∇u)+f,(x,t)∈Ω×Iu_t=\nabla\cdot(k(x)\nabla u)+f,\quad (x,t)\in\Omega\times Iu ( x , t ) = h ,u=g,x∈∂Ωu=g,x∈∂Ωu=g, \quad...
Soy un estudiante de doctorado en Informática Científica y en los últimos meses, pasé una buena cantidad de tiempo aprendiendo Python y C ++ de la manera correcta. Siento que he aprendido bien C ++ y puedo usar Python para hacer lo que quiero si mantengo un buen libro de referencia. También...
Publiqué esta pregunta anteriormente en stackoverflow, donde se cerró como fuera de tema. Espero que sobreviva aquí. En nuestro gimnasio de escalada, las rutas deben volver a establecerse de vez en cuando. Se aplican las siguientes reglas: Tenemos presas de escalada con diferentes colores en...
Hola, si alguien puede ofrecer información sobre cómo resolver mi problema, ¡sería genial! Estoy buscando calcular el área de intersección entre un volumen 3D y un plano 2D. Volumen 3D: definido por 6 puntos (siempre será una cuña 3D que se asemeja a un triángulo grueso) Plano 2D: definido por 3...
Estoy buscando encontrar una solución al siguiente problema, pero tengo problemas para formularlo con sensatez y luego encontrar un algoritmo adecuado para resolverlo. Considere una lista de artículos colocados en una bolsa de compras: 1, 2, 3, 4 ... Cada artículo puede ser parte de una o más...
¿Existe un algoritmo para reorganizar una matriz en forma de diagonal de bloque, dado que la matriz es de naturaleza diagonal de bloque pero aleatoria con una elección de base imprudente? En particular, ¿hay algún módulo de Python para
Supongamos que queremos resolver una ley de conservación hiperbólica . Realmente me gusta usar Lax-Wendroff, que dicetut+ f( u )X= 0ut+f(u)x=0u_t+f(u)_x=0 tun + 1j= unortej- Δ tΔ x( g( unortej + 1, unortej) - g( unortej, unortej - 1) )ujn+1=ujn−ΔtΔx(g(uj+1n,ujn)−g(ujn,uj−1n))u_j^{n+1} = u_j^n...
¿Cuál es una buena manera de extraer la diagonal de una matriz simétrica que ya es casi diagonal, pero donde no tiene los elementos de la matriz (solo la capacidad de aplicarla a los vectores)? Otras restricciones son, (1) aplicar la matriz n-por-n n veces para construir explícitamente la...
Tengo un código de modelado basado en el filtro Kalman que he desarrollado para una aplicación de mapeo ionosférico regional casi en tiempo real. El código asimila datos de diferentes sensores en un mapa (descrito por un conjunto de funciones básicas) utilizando un filtro Kalman. Estoy tratando de...
No estoy realmente seguro de cómo explicar este problema claramente, así que tengan paciencia conmigo. Tengo una base de 3 vectores de unidades ortonormales y una posición, una matriz de transformación 4x4 estándar en gráficos de computadora. También tengo varios puntos (compensaciones) en ese...