Changeset 356 for trunk/doc/pslib/psLibADD.tex
- Timestamp:
- Mar 31, 2004, 10:29:31 PM (22 years ago)
- File:
-
- 1 edited
-
trunk/doc/pslib/psLibADD.tex (modified) (39 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/doc/pslib/psLibADD.tex
r347 r356 1 %%% $Id: psLibADD.tex,v 1. 9 2004-04-01 04:40:14eugene Exp $1 %%% $Id: psLibADD.tex,v 1.10 2004-04-01 08:29:31 eugene Exp $ 2 2 \documentclass[panstarrs]{panstarrs} 3 3 … … 18 18 \maketitle 19 19 20 21 20 % -- Revision History -- 22 21 % provide explicit values for the old versions … … 28 27 \RevisionsEnd 29 28 30 \pagebreak 29 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 30 31 \DocumentsInternal 32 PSCD-430-xxx & PS-1 Design Reference Mission \\ \hline 33 PSCD-430-004 & Pan-STARRS IPP C Code Conventions \\ \hline 34 PSCD-430-005 & Pan-STARRS IPP SRS \\ \hline 35 PSCD-430-006 & Pan-STARRS IPP ADD \\ \hline 36 PSCD-430-008 & Pan-STARRS IPP Architecture SDR \\ 37 \DocumentsExternal 38 Posix Standard & Open Group Based Specifications Issue 6, IEEE Std 1003.1, 2003 \\ 39 \DocumentsEnd 40 31 41 \tableofcontents 32 33 42 \pagebreak 34 43 \pagenumbering{arabic} 35 44 36 \section{Pan STARRS Library PSLib}45 \section{Pan-STARRS Library PSLib} 37 46 38 47 \subsection{Math Utilities} 39 48 40 \subs ection{Sorting}49 \subsubsection{Sorting} 41 50 42 51 A variety of sorting algorithms exist, with a wide range in speed for … … 61 70 function \code{psSortIndex} shall return an integer index to the 62 71 sequence of the input array without overwriting the input array. 63 Given the following line of code, \code{out = psSortIndex (NULL, 64 &in);}, the elements of the array \code{out} are in the sequence 65 \code{in.arr[out->arr[0]]} to \code{in.arr[out->arr[in.n - 1]]}. 66 67 \subsection{Smoothing: Boxcar and Gaussian} 72 Given the following line of code: 73 \begin{verbatim} 74 out = psSortIndex (NULL,&in);} 75 \end{verbatim} 76 the elements of the array \code{out} are in the sequence 77 \code{in.arr[out->arr[0]]} to \code{in.arr[out->arr[in.n - 1]]}. 78 79 \subsubsection{Smoothing: Boxcar and Gaussian} 68 80 \label{smooth} 69 81 … … 92 104 \end{equation} 93 105 94 \subs ection{Statistics}106 \subsubsection{Statistics} 95 107 96 108 The general statistics function \code{psStats} performs a variety of … … 105 117 sample and robust estimators. 106 118 107 \ subsubsection{Sample Statistics}119 \paragraph{Sample Statistics} 108 120 109 121 We define the following statistical terms, assuming there is a set of 110 122 data elements $x_i$. 111 123 112 \ paragraph{Mean}124 \subparagraph{Mean} 113 125 114 126 The mean is defined as: … … 117 129 \end{equation} 118 130 119 \ paragraph{Median}131 \subparagraph{Median} 120 132 121 133 The median is defined as the value for which 50\% of the data values … … 127 139 robust median is quicker and more accurate. 128 140 129 \ paragraph{Upper and Lower Quartiles}141 \subparagraph{Upper and Lower Quartiles} 130 142 131 143 The upper and lower quartiles ($U_{\frac{1}{4}}$ and … … 141 153 are quicker and more accurate. 142 154 143 \ paragraph{Standard Deviation}155 \subparagraph{Standard Deviation} 144 156 145 157 The standard deviation of the sample is given by: … … 156 168 \end{equation} 157 169 158 \ subsubsection{Clipped Statistics}170 \paragraph{Clipped Statistics} 159 171 160 172 The clipped statistics are used to determine the mean and standard … … 186 198 \end{enumerate} 187 199 188 \ subsubsection{Robust Statistics}200 \paragraph{Robust Statistics} 189 201 190 202 The robust version of the statistics provides estimators of basic … … 232 244 233 245 234 \subs ection{Matrix Operations}246 \subsubsection{Matrix Operations} 235 247 236 248 In this section, we define the linear algebra operations performed on … … 251 263 \code{gsl_linalg_LU_decomp}. 252 264 253 \ subsubsection{LU Decomposition}265 \paragraph{LU Decomposition} 254 266 \label{LUdecomp} 255 267 … … 292 304 \end{equation} 293 305 294 \ subsubsection{Calculate a matrix determinant}306 \paragraph{Calculate a matrix determinant} 295 307 296 308 The determinant $D$ of a matrix $a_{ij}$ is calculated from the … … 308 320 shall be used. 309 321 310 \ subsubsection{Solving a Linear Equation}322 \paragraph{Solving a Linear Equation} 311 323 312 324 The LU decomposition of a matrix may be used to solve the … … 328 340 \end{eqnarray} 329 341 330 \ subsubsection{Invert a matrix}342 \paragraph{Invert a matrix} 331 343 332 344 Inversion of a matrix using the LU decomposition is performed by … … 337 349 operation shall be implemented using the GSL function \code{gsl_linalg_LU_invert}. 338 350 339 \ subsubsection{Perform matrix addition, subtraction and multiplication}351 \paragraph{Perform matrix addition, subtraction and multiplication} 340 352 341 353 Matrix binary arithmetic operations differ from image binary … … 361 373 \times$. 362 374 363 \ subsubsection{Transpose a matrix}375 \paragraph{Transpose a matrix} 364 376 365 377 The transpose of a matrix is simply the reorganization of the matrix … … 374 386 where $M_{ij}$ is the matrix to be transposed. 375 387 376 \ subsubsection{Convert a matrix to a vector}388 \paragraph{Convert a matrix to a vector} 377 389 378 390 Matrix-to-vector conversion is only defined for a matrix that has a … … 383 395 matrix is converted to a \code{PS_DIMEN_TRANV}-type vector. 384 396 385 \subs ection{Fitting}386 387 \ subsubsection{Chi-squared}397 \subsubsection{Fitting} 398 399 \paragraph{Chi-squared} 388 400 \label{chisq} 389 401 … … 396 408 \end{equation} 397 409 398 \ subsubsection{General Polynomial Fitting}410 \paragraph{General Polynomial Fitting} 399 411 400 412 Given a set of data values $y_i$ with errors $\sigma_i$, related to … … 421 433 (section~\ref{LUdecomp}). 422 434 423 \ subsubsection{Non-linear Fitting: Levenberg-Marquardt Method}435 \paragraph{Non-linear Fitting: Levenberg-Marquardt Method} 424 436 425 437 \TBD{describe LMM for psMinimize and psMinimizeChi2} … … 497 509 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 498 510 499 \subs ection{Polynomials}511 \subsubsection{Polynomials} 500 512 501 513 We will employ Chebyshev polynomials (NR \S 5.8) to approximate functions: … … 534 546 $-1 < x < 1$. 535 547 536 \ subsubsection{Multi-dimensional polynomials}548 \paragraph{Multi-dimensional polynomials} 537 549 538 550 Multi-dimensional polynomials shall be composed of multiplications of … … 540 552 the appropriate rank. 541 553 542 543 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 544 545 \subsection{(Fast) Fourier Transforms} 554 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 555 556 \subsubsection{(Fast) Fourier Transforms} 546 557 547 558 (Fast) Fourier Transforms (FFTs) shall be implemented using the … … 549 560 library}. 550 561 551 \ subsubsection{FFTW Plans}562 \paragraph{FFTW Plans} 552 563 553 564 FFTW requires the user to create a ``plan'' for each transform size, … … 566 577 saved at the conclusion (\code{psFFTDone()}). 567 578 568 \ subsubsection{Function mapping}579 \paragraph{Function mapping} 569 580 570 581 The forward and reverse transforms call the corresponding … … 583 594 place to avoid the need to pad the input array to hold the output. 584 595 585 \ subsubsection{More Complicated Functions}596 \paragraph{More Complicated Functions} 586 597 587 598 The \code{psFFTFilter()} and \code{psFFTFilterComplex()} functions provide … … 665 676 \end{eqnarray} 666 677 667 \ subsubsubsection{Gnomonic}678 \paragraph{Gnomonic} 668 679 669 680 The Gnomonic projection (``TAN'') is a zenithal projection. … … 674 685 \end{eqnarray} 675 686 676 \ subsubsubsection{Orthographic}687 \paragraph{Orthographic} 677 688 678 689 The Orthographic projection (``SIN'') is a zenithal projection. … … 683 694 \end{eqnarray} 684 695 685 \ subsubsubsection{Cartesian}696 \paragraph{Cartesian} 686 697 687 698 The Cartesian projection (``CAR'') is a very simple cylindrical projection. … … 692 703 \end{eqnarray} 693 704 694 \ subsubsubsection{Mercator}705 \paragraph{Mercator} 695 706 696 707 The Mercator projection (``MER'') is a cylindrical projection. … … 702 713 \end{eqnarray} 703 714 704 \ subsubsubsection{Hammer-Aitoff}715 \paragraph{Hammer-Aitoff} 705 716 706 717 The Hammer-Aitoff projection is a general projection, and is defined: … … 722 733 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 723 734 724 \subs ection{Tangent Plane to Sky}735 \subsubsection{Tangent Plane to Sky} 725 736 726 737 Mappings between the tangent plane and the sky will be implemented … … 755 766 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 756 767 757 \subs ection{The One-to-Many Problem with Mosaic Cameras}768 \subsubsection{The One-to-Many Problem with Mosaic Cameras} 758 769 759 770 The \PS{} focal plane consists of several chips, so we will often want … … 777 788 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 778 789 779 \subs ection{General Astronomy Functions}790 \subsubsection{General Astronomy Functions} 780 791 781 792 The airmass is calculated using the SLALIB function \code{sla_AIRMAS}. … … 807 818 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 808 819 809 \subs ection{Positions of Major Solar System Objects}820 \subsubsection{Positions of Major Solar System Objects} 810 821 811 822 The SLALIB function \code{SLA_RDPLAN} returns the apparent position of … … 818 829 \rightarrow \alpha + 12 {\rm hrs}$ and $\delta \rightarrow -\delta$). 819 830 820 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 821 822 \subsection{Offsets} 823 824 825 826 827 828 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 829 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 830 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 831 831 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 832 833 \subsubsection{Offsets} 834 835 836 837 838 839 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 840 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 841 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 842 843 844 \end{document} 832 845 833 846 \section{Modules} … … 866 879 \section{Architectual Components} 867 880 868 \end{document}
Note:
See TracChangeset
for help on using the changeset viewer.
