(I need to finish an introduction. Read paper in misc
.)
from weleber_tobler import kineticfield
kineticfield.compute([(54.7, 1.0), (51, 14.9), ... ])
The result is a dictionary that contains the solid angle and the percent of a sphere it occupies:
{
'percent_of_sphere': 12.368853569901212,
'steradians': 1.5543159803411815
}
No idea why the computed steradians are "close enough" but don't match those in the paper for the given input sets (see test_weleber_tobler.py
). The equation used seems right:
- Here's a nice introduction to Visual Fields.
- Used Chipmunk Basic for the original code.
- Used this manual as reference for GW-BASIC.
This line means "go to line 400"; it doesn't reassign the array value (d'oh...):
IF RD(I)<90 THEN 400
This means "use double-precision":
RD(I)=90#