Recentemente, comparei diferentes solucionadores não lineares do scipy e fiquei particularmente impressionado com o exemplo de Newton-Krylov no Scipy Cookbook, no qual eles resolvem uma equação de equação diferencial de segunda ordem com o termo de reação não linear em cerca de 20 linhas de...
Eu implementei um solucionador de Euler para trás no python 3 (usando numpy). Para minha própria conveniência e como exercício, também escrevi uma pequena função que calcula uma aproximação de diferença finita do gradiente para que nem sempre seja necessário determinar analiticamente o jacobiano...
Estou tentando descobrir se existe uma maneira mais rápida de calcular todos os autovalores e autovetores de uma matriz de adjacência muito grande e esparsa do que usar scipy.sparse.linalg.eigsh Até onde eu sei, esse método usa apenas a escassez e atributos de simetria da matriz. Uma matriz de...
Eu tenho escrito uma caixa de ferramentas do sistema de controle a partir do zero e puramente em Python3 (plug descarado:) harold. Nas minhas pesquisas anteriores, sempre tenho reclamações sobre o solucionador de Riccati care.mpor motivos técnicos / irrelevantes. Por isso, tenho escrito meu...
Estou procurando o caminho mínimo entre os mínimos de uma superfície potencial que já é conhecida em uma grade. (fonte: http://www.math.nus.edu.sg/~matrw/string/ ) Qualquer ponto no caminho é o mínimo potencial em todas as direções perpendicular ao caminho. Existe algum método SciPy ou...