Ciencia computacional

Estoy buscando un método o algoritmo para producir un valor que describa cuán diferentes son dos mallas geométricamente pero que tienen topologías diferentes. Un ejemplo serían algunos datos CAD a los que se les ha aplicado una malla FE y la parte posterior a su fabricación y escaneo en una …

8 optimization  finite-element 

Para un protector de pantalla tonto que estoy tratando de desarrollar, me gustaría generar aleatoriamente una matriz 2D de vectores 2D sin divergencia, y luego usarla para generar un gráfico de convolución integral de línea. He escuchado que una forma de hacer esto es generar ruido aleatorio y luego proyectar …

8 fourier-analysis 

Ahora que Formulize / Eureqa comenzó a cobrar $ 2500 al año por usarlo y haber estropeado la versión de prueba, ¿alguien sabe de algún reemplazo que pueda hacer cosas similares como encontrar una ecuación dados los puntos de datos? Era gratis, pero ahora está fuera del rango de precios …

8 machine-learning  regression  symbolic-computation 

dada una matriz simétrica , y una matriz arbitraria X ∈ R n × n , y un vector v ∈ R n × 1 , ¿es posible calcular la siguiente expresión en el tiempo O ( n 2 ) ?Y∈Rn×nY∈Rn×nY \in \mathbb{R}^{n \times n}X∈Rn×nX∈Rn×nX \in \mathbb{R}^{n \times n}v∈Rn×1v∈Rn×1v \in …

8 optimization  algorithms  performance  matrix  complexity 

Al escribir el código MPI 3.0, ¿debería pasar argcy argva la MPI_Initllamada o no, y por qué? EDITAR: Dado que la respuesta de Bill Barth planteó algunas preguntas, quiero hacer algunos comentarios: Pasar argc/ argvno es necesario desde MPI 1.1. La pregunta es específicamente sobre por qué debería / no …

8 c++  c  mpi 

Durante mi deambulación en Mathematica.se , gradualmente me di cuenta de que un cierto tipo de problema de resolución de ecuaciones diferenciales nos "preocupa" todo el tiempo, es decir, el problema del valor límite (BVP) de las ecuaciones diferenciales ordinarias no lineales (EDO). El método de disparo , que es …

8 ode  nonlinear-equations 

Quiero publicar el paquete de software que he escrito para mi trabajo de posgrado junto con un documento que describe el paquete. Hasta ahora he estado considerando licencias liberales como BSD. Sin embargo, ahora sé que publicar su código puede implicar un acuerdo con los términos de licencia de la …

8 publications 

Estoy entrando en mi tercer año de mi programa de doctorado, y hasta ahora mi código de investigación (problemas numéricos inversos / imagen médica / procesamiento de imágenes / etc.) consiste principalmente en scripts y funciones MATLAB desorganizados, con algunos C ++ y CUDA MEX'd rutinas lanzadas para un poco …

8 open-source 

He estado leyendo muchos documentos sobre ADMM últimamente, y también intenté resolver varios problemas al usarlo, en todos los cuales fue muy efectivo. A diferencia de otros métodos de optimización, no puedo tener una buena intuición de cómo y por qué este método es tan efectivo (por supuesto, he visto …

8 optimization  admm 

Estoy preguntando sobre la solución de la ecuación de Helmholtz en un dominio periódico con velocidad de onda constante por partes en diferentes regímenes de frecuencia. Un enfoque posible para resolver este problema es escribir ecuaciones integrales en las superficies límite en términos de la función del sistema de Green. …

8 integral-equations  helmholtz-equation 

Codifiqué un problema de Stokes Flow usando elementos finitos y estoy en el proceso de verificar que funcione. Simplemente no estoy seguro de qué tasa de convergencia debería esperar ya que refino globalmente la malla. Sé que para problemas escalares usando funciones de base lineal esperaría convergencia de orden ( …

8 pde  finite-element  convergence 

Necesito calcular la integral yo= ∫R0 0F( r ) Jnorte( zn mrR) rdrI=∫0Rf(r)Jn(znmrR)rdrI = \int_0^R f(r)J_n\left(\frac{z_{nm}r}{R}\right)rdr donde es el n t h para funciones de Bessel del primer tipo, z n m es su m t h cero y f ( r ) es una función real que es algo …

8 quadrature  integration  special-functions  numerics 

Antecedentes: Actualmente estoy ejecutando una gran cantidad de experimentos de variación de parámetros. Se están ejecutando en Python 2.6+, usando numpy. Estos experimentos demorarán aproximadamente 2 semanas en ejecutarse. Aproximadamente estoy variando 3 parámetros (variables independientes) en un rango de valores. Estoy arreglando 6 variables independientes adicionales (por ahora) estoy …

8 software  visualization  data-analysis 

Hice una pregunta similar en physics.stackexchange , siendo ignorante sobre este sitio web. Básicamente estoy buscando una manera eficiente de implementar campos gravitacionales. Tengo un gran espacio 2D, con miles de objetos en él. Luego necesito simular cómo estos objetos se ven afectados por la gravedad del otro. Pensé que …

8 optimization  simulation  approximation 

Cuando aprendí sobre SOR, se dio principalmente como uno de los primeros ejemplos de métodos iterativos, y luego los métodos iterativos que terminaría usando serían los métodos del subespacio de Krylov. ¿Alguno de los métodos iterativos como Gauss-Seidel y SOR se ha usado alguna vez en la práctica? ¿Conoces algún …

8 linear-algebra