LINUX.ORG.RU

История изменений

Исправление thunar, (текущая версия) :

Вот жеж треш! Походу, я первый раз где то опечатался и от того оно у меня и не сходилось, сейчас заново переписал в стиле: 

phi0 = phi[ix,iy]
if   iy == 0 or iy == ny-1:
	phi[ix,iy] = (0.25*(phi[ix-1,0] + phi[ix+1,0] + phi[ix,ny-2] + phi[ix,1])+ q[ix,0])*w_relax + (1.-w_relax)*phi0
else:
	phi[ix,iy] = (0.25*(phi[ix-1,iy] + phi[ix+1,iy] + phi[ix,iy-1] + phi[ix,iy+1])+ q[ix,iy])*w_relax + (1.-w_relax)*phi0
и заработало... правда сходится медленно по сравнению с обычными границами.

Исправление thunar, :

Вот жеж треш, у меня походу я первый раз где то опечатался и от того оно у меня и не сходилось, сейчас заново переписал в стиле: 

phi0 = phi[ix,iy]
if   iy == 0 or iy == ny-1:
	phi[ix,iy] = (0.25*(phi[ix-1,0] + phi[ix+1,0] + phi[ix,ny-2] + phi[ix,1])+ q[ix,0])*w_relax + (1.-w_relax)*phi0
else:
	phi[ix,iy] = (0.25*(phi[ix-1,iy] + phi[ix+1,iy] + phi[ix,iy-1] + phi[ix,iy+1])+ q[ix,iy])*w_relax + (1.-w_relax)*phi0
и заработало... правда сходится медленно по сравнению с обычными границами.

Исходная версия thunar, :

Вот жеж треш, у меня походу я первый раз где то опечатался и от того оно у меня и не сходилось, сейчас заново переписал в стиле: 

					if   iy == 0 or iy == ny-1:
						phi[ix,iy] = (0.25*(phi[ix-1,0] + phi[ix+1,0] + phi[ix,ny-2] + phi[ix,1])+ q[ix,0])*w_relax + (1.-w_relax)*phi0
					else:
						phi[ix,iy] = (0.25*(phi[ix-1,iy] + phi[ix+1,iy] + phi[ix,iy-1] + phi[ix,iy+1])+ q[ix,iy])*w_relax + (1.-w_relax)*phi0
и заработало... правда сходится медленно по сравнению с обычными границами.