ajusta opções para solucionadores de EDO
odeoptions()
Esta função exibe interativamente um comando que deveria ser
executado para se ajustar diversas opções para solucionadores de EDOs. A
variável global %ODEOPTIONS ajusta as opções.
CUIDADO: a função ode verifica se esta variável
existe e, neste caso, a utiliza. Para se utilizar valores padrões, deve-se
limpar esta variável. Note que odeoptions não cria esta
variável. Para criá-la, você deve executar a linha de comando exibida por
odeoptions.
A variável %ODEOPTIONS é um vetor com os
seguintes elementos:
[itask, tcrit, h0, hmax, hmin, jactyp, mxstep, maxordn, maxords, ixpr, ml, mu]
O valor padrão é:
[1, 0, 0, %inf, 0, 2, 500, 12, 5, 0, -1, -1]
O significado de cada elemento é descrito abaixo
itask 1 : computação normal nos tempos
especificados 2 : computação nos tempos de malha (dados na primeira linha
da saída de ode) 3 : um passo em um ponto interno da
malha e retorna 4 : computação normal sem overshooting ("desvalorização
excessiva") tcrit 5 : um passo, sem passagem de
tcrit, e retorna
tcrit assume itask igual a 4
ou 5, descrito acima
h0 primeiro passo tentado
hmax tamanho máximo do passo
hmin tamanho mínimo do passo
jactype 0 : iterações funcionais, nenhum
jacobiano usado ("adams" ou "stiff"
apenas) 1 : jacobiano completo fornecido pelo usuário 2 : jacobiano
completo internamente gerado 3 : jacobiano diagonal internamente gerado
("adams" ou "stiff" apenas) 4 :
jacobiano com bandas fornecido pelo usuário (ver ml e
mu abaixo) 5 : jacobiano com bandas internamente gerado
(ver ml e mu abaixo)
maxordn máxima ordem não-stiff permitida, no
máximo 12
maxords máxima ordem stiff permitida, no máximo
5
ixpr nível de impressão, 0 ou 1
ml,mu se
jactype é igual a 4 ou 5, ml e
mu são as semi-bandas inferior e superior do jacobiano
com bandas: a banda são os i,j's com i-ml <= j <= ny-1. Se
jactype é igual a 4 a função jacobiana deve retornar
uma matriz J que é ml+mu+1 x ny (onde ny=dim de y em ydot=f(t,y)) tal que
a coluna 1 de J é feita de mu zeros seguidos por df1/dy1, df2/dy1,
df3/dy1, ... (1+ml entradas não-nulas possivelmente) a coluna 2 é feita de
mu-1 zeros seguidos por df1/dx2, df2/dx2, etc