All Samples(1787) | Call(1617) | Derive(0) | Import(170)
radians(x) Convert angle x from degrees to radians.
src/p/y/pyx-HEAD/trunk/pyx/examples/drawing2/insert.py pyx(Download)
from math import sin, cos, radians from pyx import * angle = 10 factor = 1.0 / (cos(radians(angle)) + sin(radians(angle))) cc = canvas.canvas()
src/p/y/pyx-HEAD/pyx/examples/drawing2/insert.py pyx(Download)
from math import sin, cos, radians from pyx import * angle = 10 factor = 1.0 / (cos(radians(angle)) + sin(radians(angle))) cc = canvas.canvas()
src/o/w/owyl-0.3/examples/boids.py owyl(Download)
import os import random from math import radians, degrees, sin, cos, pi, atan2 pi_2 = pi*2.0 pi_1_2 = pi/2.0 pi_1_4 = pi/4.0
def move(self, **kwargs):
"""Move the actor forward perpetually.
"""
bb = kwargs['blackboard']
while True:
dt = bb['dt']
r = radians(getR(self)) # rotation
dx = gx-self.x
dy = gy-self.y
seek_heading = self.getFacing(dx, dy)
my_heading = radians(self.rotation)
rsize = degrees(self.findRotationDelta(my_heading, seek_heading))
rate = kwargs['rate']
while True:
dt = bb['dt'] or 0.01
n_heading = radians(self.findAverageHeading(*self.neighbors))
if n_heading is None:
yield None
continue
my_heading = radians(self.rotation)
heading_away_from_neighbors = self.getFacing(dx, dy)
flee_heading = heading_away_from_neighbors
my_heading = radians(self.rotation)
rsize = degrees(self.findRotationDelta(my_heading, flee_heading))
dx = np_x-self.x
dy = np_y-self.y
seek_heading = self.getFacing(dx, dy)
my_heading = radians(self.rotation)
# Find the rotation delta
rsize = degrees(self.findRotationDelta(my_heading, seek_heading))
src/r/a/Rabbyt-0.8.3/examples/driving.py Rabbyt(Download)
from __future__ import division import pygame import rabbyt from math import cos, sin, radians import random
a = [0.0,0.0]
if self.boost_endtime > rabbyt.get_time():
f = 3*(self.boost_endtime - rabbyt.get_time())/self.boost_length
a[0] += cos(radians(self.boost_rot))*f
a[1] += sin(radians(self.boost_rot))*f
self.create_boost_particle()
if self.accelerating:
a[0] += cos(radians(self.rot))*.9
a[1] += sin(radians(self.rot))*.9
src/o/w/owyl-0.3/examples/testboids.py owyl(Download)
# -*- coding: utf-8 -*- """ """ import unittest from math import radians, pi
b.rotation = 0
pos = (0, 5)
expect = radians(0)
rr = b.getFacing(*pos)
self.assertEqual(rr, expect,
"%s: got %s, expected %s" % (pos, rr, expect))
pos = (5, 0)
expect = radians(90)
"%s: got %s, expected %s" % (pos, rr, expect))
pos = (-5, 0)
expect = radians(-90)
rr = b.getFacing(*pos)
self.assertEqual(rr, expect,
"%s: got %s, expected %s" % (pos, rr, expect))
pos = (0, -5)
expect = radians(180)
)
for current, match, delta in current_match_delta:
current = radians(current)
match = radians(match)
delta = radians(delta)
dr = b.findRotationDelta(current, match)
src/r/a/Rabbyt-0.8.3/examples/pymunk_integration.py Rabbyt(Download)
import pygame import rabbyt from math import cos, sin, radians, degrees, pi import random import os.path
src/n/e/netpylab-HEAD/netpylab/paths.py netpylab(Download)
import time from math import log, cos, sin, radians, degrees, atan2, tan, pi, sqrt from xml.dom import minidom import bisect pow2_25 = 2**25
def osm2map(self):
x = (self.lon + 180.0) / 360.0
y = (1 - log(tan(radians(self.lat)) + 1/cos(radians(self.lat))) / pi) / 2
return(x, y)
def geo2map(self):
return (2*self.xo-1)*pow2_25, -(2*self.yo-1)*pow2_25
def get_distance(self, other):
radius = 6371000
deltaLat = other.lat-self.lat
deltaLong = other.lon-self.lon
a = (sin(radians(deltaLat)/2))**2
b = cos(radians(other.lat))*cos(radians(self.lat))*(sin(radians(deltaLong)/2))**2
def get_direction(self, other):
deltaLong = other.lon - self.lon
a = sin(radians(deltaLong))*cos(radians(self.lat))
b = cos(radians(other.lat))*sin(radians(self.lat))
c = sin(radians(other.lat))*cos(radians(self.lat))*cos(radians(deltaLong))
direction = degrees(atan2(a, b-c))
return direction #maybe a sign problem ????
src/a/g/agtl-0.7.1.0-freerunner0/advancedcaching/astral.py agtl(Download)
# Shortened for AGTL by Daniel Fett import datetime from math import cos, sin, tan, acos, asin, atan2, floor, radians, degrees SUN_POSITION_CACHE_DURATION = 3600 # seconds
if hourangle < -180:
hourangle = hourangle + 360.0
harad = radians(hourangle)
csz = sin(radians(latitude)) * sin(radians(solarDec)) + \
cos(radians(latitude)) * cos(radians(solarDec)) * cos(harad)
zenith = degrees(acos(csz))
azDenom = (cos(radians(latitude)) * sin(radians(zenith)))
if (abs(azDenom) > 0.001):
azRad = ((sin(radians(latitude)) * cos(radians(zenith))) - sin(radians(solarDec))) / azDenom
if hourangle < -180:
hourangle = hourangle + 360.0
harad = radians(hourangle)
csz = sin(radians(latitude)) * sin(radians(solarDec)) + \
cos(radians(latitude)) * cos(radians(solarDec)) * cos(harad)
zenith = degrees(acos(csz))
azDenom = (cos(radians(latitude)) * sin(radians(zenith)))
if (abs(azDenom) > 0.001):
azRad = ((sin(radians(latitude)) * cos(radians(zenith))) - sin(radians(solarDec))) / azDenom
if exoatmElevation > 85.0:
refractionCorrection = 0.0
else:
te = tan(radians(exoatmElevation))
if exoatmElevation > 5.0:
refractionCorrection = 58.1 / te - 0.07 / (te * te * te) + 0.000086 / (te * te * te * te * te)
elif exoatmElevation > -0.575:
def _obliquity_correction(self, juliancentury):
e0 = self._mean_obliquity_of_ecliptic(juliancentury)
omega = 125.04 - 1934.136 * juliancentury
return e0 + 0.00256 * cos(radians(omega))
def _geom_mean_long_sun(self, juliancentury):
def _eq_of_time(self, juliancentury):
epsilon = self._obliquity_correction(juliancentury)
l0 = self._geom_mean_long_sun(juliancentury)
e = self._eccentricity_earth_orbit(juliancentury)
m = self._geom_mean_anomaly_sun(juliancentury)
y = tan(radians(epsilon) / 2.0)
y = y * y
sin2l0 = sin(2.0 * radians(l0))
sinm = sin(radians(m))
cos2l0 = cos(2.0 * radians(l0))
sin4l0 = sin(4.0 * radians(l0))
sin2m = sin(2.0 * radians(m))
def _sun_eq_of_center(self, juliancentury):
m = self._geom_mean_anomaly_sun(juliancentury)
mrad = radians(m)
sinm = sin(mrad)
sin2m = sin(mrad + mrad)
sin3m = sin(mrad + mrad + mrad)
def _sun_apparent_long(self, juliancentury):
O = self._sun_true_long(juliancentury)
omega = 125.04 - 1934.136 * juliancentury
return O - 0.00569 - 0.00478 * sin(radians(omega))
def _sun_declination(self, juliancentury):
e = self._obliquity_correction(juliancentury)
lambd = self._sun_apparent_long(juliancentury)
sint = sin(radians(e)) * sin(radians(lambd))
def _hour_angle(self, latitude, solar_dec, solar_depression):
latRad = radians(latitude)
sdRad = radians(solar_dec)
HA = (acos(cos(radians(90 + solar_depression)) / (cos(latRad) * cos(sdRad)) - tan(latRad) * tan(sdRad)))
return HA
def _sun_rad_vector(self, juliancentury):
v = self._sun_true_anomoly(juliancentury)
e = self._eccentricity_earth_orbit(juliancentury)
return (1.000001018 * (1 - e * e)) / (1 + e * cos(radians(v)))
def _sun_rt_ascension(self, juliancentury):
e = self._obliquity_correction(juliancentury)
lambd = self._sun_apparent_long(juliancentury)
tananum = (cos(radians(e)) * sin(radians(lambd)))
tanadenom = (cos(radians(lambd)))
src/w/e/weewx-HEAD/trunk/experimental/astral.py weewx(Download)
import datetime import time from math import cos, sin, tan, acos, asin, atan2, floor, radians, degrees __version__ = "0.3+" __author__ = "Simon Kennedy <python@sffjunkie.co.uk>"
hourangle = ha + eqtime/4.0
#print "hour angle = ", hourangle
harad = radians(hourangle)
csz = sin(radians(latitude)) * sin(radians(solarDec)) + \
cos(radians(latitude)) * cos(radians(solarDec)) * cos(harad)
zenith = degrees(acos(csz))
azDenom = (cos(radians(latitude)) * sin(radians(zenith)))
if (abs(azDenom) > 0.001):
azRad = ((sin(radians(latitude)) * cos(radians(zenith))) - sin(radians(solarDec))) / azDenom
if exoatmElevation > 85.0:
refractionCorrection = 0.0
else:
te = tan(radians(exoatmElevation))
if exoatmElevation > 5.0:
refractionCorrection = 58.1 / te - 0.07 / (te * te * te) + 0.000086 / (te * te * te * te * te)
elif exoatmElevation > -0.575:
def _obliquity_correction(juliancentury):
e0 = _mean_obliquity_of_ecliptic(juliancentury)
omega = 125.04 - 1934.136 * juliancentury
return e0 + 0.00256 * cos(radians(omega))
def _geom_mean_long_sun(juliancentury):
def _eq_of_time(juliancentury):
epsilon = _obliquity_correction(juliancentury)
l0 = _geom_mean_long_sun(juliancentury)
e = _eccentricity_earth_orbit(juliancentury)
m = _geom_mean_anomaly_sun(juliancentury)
y = tan(radians(epsilon) / 2.0)
y = y * y
sin2l0 = sin(2.0 * radians(l0))
sinm = sin(radians(m))
cos2l0 = cos(2.0 * radians(l0))
sin4l0 = sin(4.0 * radians(l0))
sin2m = sin(2.0 * radians(m))
def _sun_eq_of_center(juliancentury):
m = _geom_mean_anomaly_sun(juliancentury)
mrad = radians(m)
sinm = sin(mrad)
sin2m = sin(mrad + mrad)
sin3m = sin(mrad + mrad + mrad)
def _sun_apparent_long(juliancentury):
O = _sun_true_long(juliancentury)
omega = 125.04 - 1934.136 * juliancentury
return O - 0.00569 - 0.00478 * sin(radians(omega))
def _sun_declination(juliancentury):
e = _obliquity_correction(juliancentury)
lambd = _sun_apparent_long(juliancentury)
sint = sin(radians(e)) * sin(radians(lambd))
def _hour_angle(latitude, solar_dec, solar_depression):
latRad = radians(latitude)
sdRad = radians(solar_dec)
HA = (acos(cos(radians(90 + solar_depression)) / (cos(latRad) * cos(sdRad)) - tan(latRad) * tan(sdRad)))
return HA
def _sun_rad_vector(juliancentury):
v = _sun_true_anomoly(juliancentury)
e = _eccentricity_earth_orbit(juliancentury)
return (1.000001018 * (1 - e * e)) / (1 + e * cos(radians(v)))
def _sun_rt_ascension(juliancentury):
e = _obliquity_correction(juliancentury)
lambd = _sun_apparent_long(juliancentury)
tananum = (cos(radians(e)) * sin(radians(lambd)))
tanadenom = (cos(radians(lambd)))
src/a/s/astral-0.3/src/astral.py astral(Download)
"""
import datetime
from math import cos, sin, tan, acos, asin, atan2, floor, radians, degrees
try:
import pytz
if hourangle < -180:
hourangle = hourangle + 360.0
harad = radians(hourangle)
csz = sin(radians(latitude)) * sin(radians(solarDec)) + \
cos(radians(latitude)) * cos(radians(solarDec)) * cos(harad)
zenith = degrees(acos(csz))
azDenom = (cos(radians(latitude)) * sin(radians(zenith)))
if (abs(azDenom) > 0.001):
azRad = ((sin(radians(latitude)) * cos(radians(zenith))) - sin(radians(solarDec))) / azDenom
if hourangle < -180:
hourangle = hourangle + 360.0
harad = radians(hourangle)
csz = sin(radians(latitude)) * sin(radians(solarDec)) + \
cos(radians(latitude)) * cos(radians(solarDec)) * cos(harad)
zenith = degrees(acos(csz))
azDenom = (cos(radians(latitude)) * sin(radians(zenith)))
if (abs(azDenom) > 0.001):
azRad = ((sin(radians(latitude)) * cos(radians(zenith))) - sin(radians(solarDec))) / azDenom
if exoatmElevation > 85.0:
refractionCorrection = 0.0
else:
te = tan(radians(exoatmElevation))
if exoatmElevation > 5.0:
refractionCorrection = 58.1 / te - 0.07 / (te * te * te) + 0.000086 / (te * te * te * te * te)
elif exoatmElevation > -0.575:
def _obliquity_correction(self, juliancentury):
e0 = self._mean_obliquity_of_ecliptic(juliancentury)
omega = 125.04 - 1934.136 * juliancentury
return e0 + 0.00256 * cos(radians(omega))
def _geom_mean_long_sun(self, juliancentury):
def _eq_of_time(self, juliancentury):
epsilon = self._obliquity_correction(juliancentury)
l0 = self._geom_mean_long_sun(juliancentury)
e = self._eccentricity_earth_orbit(juliancentury)
m = self._geom_mean_anomaly_sun(juliancentury)
y = tan(radians(epsilon) / 2.0)
y = y * y
sin2l0 = sin(2.0 * radians(l0))
sinm = sin(radians(m))
cos2l0 = cos(2.0 * radians(l0))
sin4l0 = sin(4.0 * radians(l0))
sin2m = sin(2.0 * radians(m))
def _sun_eq_of_center(self, juliancentury):
m = self._geom_mean_anomaly_sun(juliancentury)
mrad = radians(m)
sinm = sin(mrad)
sin2m = sin(mrad + mrad)
sin3m = sin(mrad + mrad + mrad)
def _sun_apparent_long(self, juliancentury):
O = self._sun_true_long(juliancentury)
omega = 125.04 - 1934.136 * juliancentury
return O - 0.00569 - 0.00478 * sin(radians(omega))
def _sun_declination(self, juliancentury):
e = self._obliquity_correction(juliancentury)
lambd = self._sun_apparent_long(juliancentury)
sint = sin(radians(e)) * sin(radians(lambd))
def _hour_angle(self, latitude, solar_dec, solar_depression):
latRad = radians(latitude)
sdRad = radians(solar_dec)
HA = (acos(cos(radians(90 + solar_depression)) / (cos(latRad) * cos(sdRad)) - tan(latRad) * tan(sdRad)))
return HA
def _sun_rad_vector(self, juliancentury):
v = self._sun_true_anomoly(juliancentury)
e = self._eccentricity_earth_orbit(juliancentury)
return (1.000001018 * (1 - e * e)) / (1 + e * cos(radians(v)))
def _sun_rt_ascension(self, juliancentury):
e = self._obliquity_correction(juliancentury)
lambd = self._sun_apparent_long(juliancentury)
tananum = (cos(radians(e)) * sin(radians(lambd)))
tanadenom = (cos(radians(lambd)))
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 Next