testing docs

c_src/_web.c

@@ -437,7 +437,6 @@ static PyObject *web_lightcurve(PyObject *self, PyObject *args)
 
     PyObject *t_array = PyArray_FROM_OTF(t_obj, NPY_DOUBLE, NPY_IN_ARRAY);
 
-
     /* If that didn't work, throw an exception. */
     if (t_array == NULL) {
         Py_XDECREF(t_array);

c_src/web.c

@@ -25,6 +25,7 @@ double *lightcurve(int n_layers, int n_points, double *t, double tc, double per,
     double *output = malloc(sizeof(double) * n_points);
 
     // generate the planet grid
+
     double **planet = generate_planet(n_layers);
 
     double *transit_coords = separation_of_centers(tc,tc,per,a,inc,ecc,omega,a_rs,r2);
@@ -50,14 +51,13 @@ double *lightcurve(int n_layers, int n_points, double *t, double tc, double per,
 
         phase = calc_phase(t[n],tc,per);
 
-
         // make correction for finite light travel speed
 
         phase_z = old_coords[3];
         phase_dz = transit_z-phase_z;
         phase_dt = (phase_dz/c)/(3600.0*24.0);
 
-        double *substellar = calc_substellar(phase,coords);
+        double *substellar = calc_substellar(phase,old_coords);
 
         lambda0 = substellar[0];
         phi0 = substellar[1];

docs/index.rst

@@ -13,6 +13,7 @@ Contents:
 
    intro
    installation
+   quickstart
    ack
 
 Release Notes

docs/quickstart.rst

@@ -0,0 +1,21 @@
+
+Quickstart
+============
+Making a lightcurve
+---
+This is likely the basic task that you need SPIDERMAN to perform, to, for example, form the basis for a likelihood function being fed into an mcmc. Here is how to return a simulated lightcurve from SPIDERMAN, the fast way.
+
+.. figure:: images/f1.png
+    :width: 200px
+    :align: center
+    :height: 100px
+    :alt: alternate text
+    :figclass: align-center
+
+    figure are like images but with a caption
+
+    and whatever else youwish to add
+
+    .. code-block:: python
+
+        import image

spiderman/params.py

@@ -1,6 +1,7 @@
 import numpy as np
 import spiderman as sp
 import spiderman._web as _web
+import spiderman.plot as splt
 
 class ModelParams(object):
 
@@ -76,3 +77,13 @@ def calc_substellar(self,t):
 		substellar = _web.calc_substellar(self.phase,np.array(coords))
 		self.lambda0 = substellar[0]
 		self.phi0 = substellar[1]
+
+	def plot_system(self,t,ax=False,min_temp=False,max_temp=False,temp_map=False,min_bright=0):
+		return splt.plot_system(self,t,ax=ax,min_temp=min_temp,max_temp=max_temp,temp_map=temp_map,min_bright=min_bright)
+
+	def plot_planet(self,t,ax=False,min_temp=False,max_temp=False,temp_map=False,min_bright=0,scale_planet=1.0,planet_cen=[0.0,0.0]):
+		return splt.plot_planet(self,t,ax=ax,min_temp=min_temp,max_temp=max_temp,temp_map=temp_map,min_bright=min_bright,scale_planet=scale_planet,planet_cen=planet_cen)
+
+	def lightcurve(self,t):
+		brightness_params = self.format_bright_params()
+		return _web.lightcurve(self.n_layers,t,self.t0,self.per,self.a_abs,self.inc,self.ecc,self.w,self.a,self.rp,self.p_u1,self.p_u2,self.brightness_type,brightness_params)

spiderman/plot.py

@@ -48,6 +48,7 @@ def plot_system(spider_params,t,ax=False,min_temp=False,max_temp=False,temp_map=
 
 	ax.set_xlim(star_offset_pix+40*p_imrat,star_offset_pix+-40*p_imrat)
 	ax.set_ylim(star_offset_pix+-10*p_imrat,star_offset_pix+10*p_imrat)
+	return ax