Wpis z mikrobloga

@zarev:
from SDToolbox import *
import csv
import numpy as np
import pandas as pd
import math
import cantera as ct

P1 = 101325
T = [293,313,333,353,373]
mech = 'gri30.cti'
q = 'H2:2.0 O2:1.0 N2:3.76'
fname = 'wyniki dla gri30'

#for j in range(0, len(T)):
gas = Solution(mech)
gas.TPX = T[0], P1, q
r = gas.density
b = 10000; j = gas.nspecies
out = cvoutput(b, j)
out = explosion(gas, fname,
@Drahearys: Wygląda na błąd w bibliotece SDToolbox w pliku CV.py. Przekazujesz do funkcji dane, które ani razu nie wywołują któregoś z warunków temp_grad[j] >= (0.5* MAX) albo temp_grad[j] <= (0.5* MAX) i nie przypisują niczego do tstep2 i/lub tstep1.
@Drahearys: jak już przedprzedmówca zdążył zauważyć, "wcięcia driven" bez wcięć to tak słabo raczej, no ale nawet jeśli, analizując kolejne poprawne permutacje póki co zawsze widzę możliwy scenariusz że któraś ze zmiennych tstep1 lub tstep2 jest niezainicjalizowana dla wartości halfTflag1 = halfTflag2 = 0 i ważne jest jak mocno wcięta jest sporna linia (ta gdzie jest błąd)
@morsik: @DyrektorWykopu: @luzny_lori:
pastebin
mam nadzieje że jak teraz rzucicie okiem to moze cos sie wyjasni, dodam że poprzednim razem kiedy uzywałem funkcji explosion ( w innym projekcie) inputowe dane odnosnie gazu były inne tj. liczone przez inną funkcję(postschock, a wczeniej cjspeed).
Także pytanie mam takie.
Czy moze być tak że całość zadziała tylko jeżeli ta poprzednia funkcja (postshock) dostarczy dane funkcji explosion. Chodzi o to że w poprzenim
@Drahearys: Są dwie możliwości - albo funkcja explosion jest po prostu błędnie napisana i wykłada się akurat na danych, które do niej przekazujesz, albo autor założył, że funkcja będzie stosowana tylko dla pewnego zakresu danych, a Ty przekazujesz jej dane spoza tego zakresu. Za pierwszą opcją przemawia to, że tamta pętla z warunkami jest po prostu bardzo brzydko napisana i aż prosi się o takie błędy. Musisz sam sprawdzić, czy dla