En la literatura sobre métodos numéricos, muchos trabajos de investigación consisten en una descripción de una nueva variación algorítmica, seguida de algunos problemas de prueba que comparan el nuevo método con uno o dos métodos existentes. Esto dificulta determinar Cómo funcionaría el nuevo...
Estoy buscando un software de generación de mallas que es gratis y de código abierto, proporciona una interfaz de secuencia de comandos sana para la especificación de dominio, funciona para geometrías complejas, puede generar mallas 2D y 3D, ¿Que opciones
Siempre calculo las cosas a mano, pero ahora mis camaradas se están poniendo desagradables y hacen muchos ejercicios repetitivos que implican simplemente conectar cosas como la expresión anterior. Estoy particularmente interesado en software de código abierto como Python o R para simplificar este...
Supongamos que tengo alguna función y quiero encontrar tal que . Puede ser que utilice el método de Newton-Raphson. Pero esto requiere que sé que el derivado de la función . Una expresión analítica para puede no estar disponible. Por ejemplo, puede definirse mediante un código informático...
Necesito escribir algunos datos de un cálculo, que serán leídos más tarde por Paraview (archivo .vtu o vtk). Cuando se trata del tamaño del archivo, ¿debo elegir el formato ASCII o el formato
Acabo de instalar una tarjeta gráfica Nvidia GT660 en mi escritorio y, después de una lucha, me las arreglo para interactuar con R. He estado jugando con varios paquetes de R que las GPU utilizando, especialmente gputools, y yo estaba comparando el tiempo empleado por mi GPU y CPU para realizar...
Muchos métodos numéricos para PDE hiperbólicas se basan en el uso de solucionadores de Riemann. Tales solucionadores son esenciales para capturar con precisión las ondas de choque. Hay una gama de tales solucionadores disponibles para los sistemas mejor estudiados (por ejemplo, solucionadores...
Los métodos numéricos para resolver PDEs (u ODE) se dividen en dos grandes categorías: métodos explícitos e implícitos. Los métodos implícitos permiten pasos de tiempo estables más grandes pero requieren más trabajo por paso. Para las PDE hiperbólicas, la sabiduría común es que los métodos...
Dada una evolución PDE tut= A u + B utut=UNtu+situu_t = Au + Bu donde son operadores diferenciales (posiblemente no lineales) que no conmutan, un enfoque numérico común es alternar entre resolverA , BUN,siA,B tut= A utut=UNtuu_t = Au y tut= B u .tut=situ.u_t = Bu. La implementación más simple...
Estoy en un campo que no necesariamente hace una gran cantidad de trabajo de HPC, y cuando lo encuentra, a menudo es el resultado de investigadores de otros campos que exploran nuevas aplicaciones para sus métodos y similares. Principalmente, lo que esto significa es que nunca se introduce...
Esto está dirigido principalmente a PDE elípticas sobre dominios convexos, para que pueda obtener una buena visión general de los dos
Considere con casi singular, lo que significa que hay un valor propio de que es muy pequeño. El criterio habitual de detención de un método iterativo se basa en el valor residual y considera que las iteraciones pueden detenerse cuando con n el número de iteración. Pero en el caso que estamos...
Estoy haciendo una tarea de clasificación de texto con R, y obtengo una matriz de términos de documentos con un tamaño de 22490 por 120,000 (solo 4 millones de entradas distintas de cero, menos del 1% de entradas). Ahora quiero reducir la dimensionalidad utilizando PCA (Análisis de componentes...
Lo que más me disgusta de MPI es tratar con tipos de datos (es decir, mapas / máscaras de datos) porque no encajan muy bien con C ++ orientado a objetos. boost::mpisolo es compatible con MPI 1.1, sin embargo, desde su sitio web: boost :: mpi es una interfaz amigable de C ++ para la interfaz...
Estoy usando Valgrind + Callgrind para perfilar un solucionador que he escrito. Como dice el manual de usuario de Valgrind, he compilado mi código con las opciones de depuración para el compilador: "Sin información de depuración, las mejores herramientas de Valgrind serán capaces de adivinar a...
Cuando se calcula la factorización QR en la práctica, uno usa las reflexiones de Householder para poner a cero la parte inferior de una matriz. Sé que para calcular valores propios de matrices simétricas, lo mejor que puedes hacer con las reflexiones de Householder es llevarlo a una forma...
Me gustaría saber si hay una manera rápida de calcular la distancia euclidiana de dos vectores en Octave. Parece que no hay una función especial para eso, así que ¿debería usar la fórmula con
Probablemente sea una pregunta a nivel de estudiante, pero no puedo dejarla clara para mí. ¿Por qué es más preciso usar cuadrículas no uniformes en los métodos numéricos? Estoy pensando en el contexto de algún método de diferencia finita para el PDE de la forma . Y suponga que estoy interesado en...
Me gustaría ser independiente del software comercial para mi trabajo científico. Encuentro insatisfactorio una dependencia de paquetes comerciales como Matlab y sus cajas de herramientas, porque no sé si tendré acceso a Matlab en el futuro y porque no me gusta el idioma. Por lo tanto, estoy...
Recientemente he estado comparando diferentes solucionadores no lineales de scipy y quedé particularmente impresionado con el ejemplo de Newton-Krylov en el Scipy Cookbook en el que resuelven una ecuación de ecuación diferencial de segundo orden con un término de reacción no lineal en...