Tengo una pregunta simple que es realmente difícil para Google (además del canónico Lo que todo informático debe saber sobre el papel de aritmética de punto flotante ). ¿Cuándo deberían usarse funciones como log1po en expm1lugar de logy exp? ¿Cuándo no deberían ser utilizados? ¿Cómo difieren las...
¿Alguien me puede recomendar un buen sistema de gestión de flujo de trabajo (WMS), preferiblemente en Python? Hasta ahora he estado usando GNU Make, pero introduce una capa de complejidad que quiero evitar. Un buen WMS debe tener las siguientes características: integrarse fácilmente con...
Los ingenieros a menudo insisten en usar métodos localmente conservadores como el volumen finito, la diferencia finita conservadora o los métodos discontinuos de Galerkin para resolver PDEs. ¿Qué puede salir mal cuando se usa un método que no es localmente conservador? Bien, entonces la...
Antecedentes: creo que podría querer portar algún código que calcule productos de vectores exponenciales de matriz utilizando un método de subespacio Krylov de MATLAB a Python. (Específicamente, la función expmvp de Jitse Niesen , que utiliza un algoritmo descrito en este documento ). Sin embargo,...
Como alguien que tiene una licenciatura en física, me escandalizó un poco cuando comencé a trabajar con simulaciones moleculares. Fue un poco sorprendente descubrir que incluso las simulaciones más detalladas y computacionalmente caras no pueden reproducir cuantitativamente el comportamiento...
El problema Actualmente estoy trabajando en una simulación de elementos finitos Navier Stokes y me gustaría investigar los efectos de una variedad de parámetros. Algunos parámetros se especifican en un archivo de entrada o mediante las opciones de una línea de comando; otros parámetros se...
Recientemente me encontré con una formulación del meta-fenómeno : " dos es fácil, tres es difícil " (redactado de esta manera por Federico Poloni), que se puede describir de la siguiente manera: Cuando se formula un cierto problema para dos entidades, es relativamente fácil de resolver; sin...
Si tengo una matriz invertible cuadrada y tomo su determinante, y encuentro que , ¿implica esto que la matriz está mal condicionada?det(A)≈0det(A)≈0\det(A) \approx 0 ¿Lo inverso también es cierto? ¿Una matriz mal acondicionada tiene un determinante casi cero? Aquí hay algo que probé en Octave: a...
El Journal of Computational Physics ha sido una salida importante para la ciencia computacional en el pasado, y he publicado allí antes. Para el beneficio de aquellos (como yo) que firmaron el boicot de Elsevier , ¿qué revistas que no sean de Elsevier serían lugares apropiados para publicar...
¿Quieres mejorar esta publicación? Proporcione respuestas detalladas a esta pregunta, incluidas citas y una explicación de por qué su respuesta es correcta. Las respuestas sin suficiente detalle pueden ser editadas o eliminadas. La versión 3.0 del estándar MPI eliminó...
Sé y he usado f2py2e para envolver un código Fortran 77 antiguo, pero entiendo que no funciona con el código Fortran 95 más nuevo. He investigado qué debo usar, y he encontrado fwrap y G3 f2py, ninguno de los cuales parece dar una explicación de su estado actual o cómo usarlos (más allá del uso...
Quiero modificar una matriz de transición cuadrada densa en el lugar cambiando el orden de varias de sus filas y columnas, usando la biblioteca numpy de python. Matemáticamente, esto corresponde a la multiplicación previa de la matriz por la matriz de permutación P y la multiplicación posterior por...
EDITAR: Estoy probando si algún valor propio tiene una magnitud de uno o más. Necesito encontrar el valor propio absoluto más grande de una gran matriz dispersa, no simétrica. He estado usando la eigen()función de R , que usa el algoritmo QR de EISPACK o LAPACK para encontrar todos los valores...
Esta pregunta se migró de Stack Overflow porque se puede responder en Computational Science Stack Exchange. Migrado hace 7 años . ¿Existen problemas / algoritmos famosos en la informática científica que no puedan acelerarse mediante paralelización? Al leer libros sobre...
En estadística y sus diversas aplicaciones, a menudo calculamos la matriz de covarianza , que es positiva definida (en los casos considerados) y simétrica, para varios usos. A veces, necesitamos el inverso de esta matriz para varios cálculos (formas cuadráticas con este inverso como la (única)...
¿Qué puede salir mal cuando se utilizan los métodos de Krylov preconditonado de KSP ( paquete de solucionador lineal de PETSc ) para resolver un sistema lineal disperso como los obtenidos al discretizar y linealizar ecuaciones diferenciales parciales? ¿Qué pasos puedo tomar para determinar qué...
¿Existen herramientas diseñadas específicamente para comprimir datos científicos de coma flotante? Si una función es fluida, obviamente existe una gran correlación entre los números que representan esa función, por lo que los datos deberían comprimirse bien. Sin embargo, los datos binarios de...
No estoy muy familiarizado con los esquemas comunes de discretización para PDEs. Sé que Crank-Nicolson es un esquema popular para discretizar la ecuación de difusión. ¿También es una buena opción para el término de advección? Estoy interesado en resolver la ecuación Reacción-Difusión-Advección...
Uno de los principales problemas con los que tenemos que lidiar en las simulaciones moleculares es el cálculo de las fuerzas dependientes de la distancia. Si podemos restringir las funciones de fuerza y distancia para tener potencias iguales de la distancia de separación , entonces podemos...
Los desarrolladores de software tienen la opción de elegir una licencia apropiada de acuerdo con los objetivos del trabajo. ¿Alguien puede dar algunas recomendaciones / experiencias sobre qué licencia elegir para el software? ¿Cuáles son las ventajas y desventajas de "regalar" todo el trabajo...