IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 840


Ignore:
Timestamp:
Jun 2, 2004, 5:13:23 PM (22 years ago)
Author:
eugene
Message:

added the old SRS text after document end

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/doc/design/hardware.tex

    r669 r840  
    806806
    807807\end{document}
     808
     809
     810\subsection{Computer Hardware}
     811
     812\subsubsection{Overview}
     813
     814This section discusses the Pan-STARRS Image Processing Pipeline (IPP)
     815PS-1 hardware requirements.  The hardware requirements addressed in
     816this section consist of:
     817
     818\begin{itemize}
     819\item Total Disk Volume
     820\item Total Processing Power
     821\item Sustained Switch Bandwidth
     822\item Sustained Node Network I/O
     823\item Sustained Disk I/O
     824\end{itemize}
     825
     826Even without the complete IPP design, it is possible to identify the
     827major drivers on the hardware requirements.  The total disk volume
     828requirements are dominated by the need to store raw images for a
     829certain period, the need to store calibration images for a longer
     830period, and the need to store the static sky images.  Of the various
     831analysis stages, Phase 2 and Phase 4 present the most significant
     832demands in terms of data I/O throughput on the network.  Phase 2 and
     833Phase 4 also present the most significant CPU demands.  In this
     834discusion, Phase 2 refers to the per-OTA image pre-processing in which
     835the instrumental signature is removed and a first pass object
     836detection is performed.  Phase 4 refers to the multiple OTA
     837combination in which the pre-processed images are merged and combined,
     838in both addition and subtraction, with the static sky image, and up to
     839three object detection passes are performed.
     840
     841This document does not address the hardware requirements implied by
     842Phase 1 or 3, nor the load required by the calibration or reference
     843catalog creation stages.  In the first instance, the operations are
     844only performed on the metadata and are extremely minimal both in terms
     845of data I/O and computation requirements.  In the second case, the
     846processing is less time critical than the per-image processing and is
     847performed only infrequently (once per night to once per week, month or
     848year).  \tbd{The software implementation for metadata storage (RDBMS,
     849FITS tables, etc) will have a very large impact and will be evaluated
     850along with the needed hardware at a later date.}
     851
     852We will address the various hardware requirements by referring to an
     853assumed data processing and data organization scenario.  The
     854organization of the data and certain aspects of the data processing
     855scheme have very large implications for the hardware requirements.  In
     856this analysis, we assume that data types are chosen to minimize the
     857data volume and that the data is organized to minimize the I/O
     858bandwidth needs, as defined below.  We address the data requirements
     859of the single-telescope Pan-STARRS-1 scenario based on the Design
     860Reference Mission \tbd{REF}.
     861
     862\subsubsection{Data Organization}
     863
     864The IPP hardware system must provide both data storage and
     865computational resources.  The IPP requires relatively large amounts of
     866data storage space, primarily for the image data.  Image data is
     867organized in two categories.  First, there is the per-OTA data -- data
     868associated with specific OTAs, including the raw images, the
     869calibration images, and temporary processed images at various stages.
     870Second, there is the data associated with the static sky imagery,
     871which is in turn organized into smaller sky-cell units.  The first
     872assumption we make is that the hardware is organized into nodes which
     873provide both data storage and computational resources.  The second
     874assumption we make is that the data storage nodes are divided into two
     875classes: those which deal with the per-OTA data and those that provide
     876the static sky storage.  In addition, we assume that the computational
     877tasks related to Phase 2 take place on the per-OTA storage nodes and
     878the Phase 4 computation takes place on the static sky storage nodes.
     879
     880Figure~\ref{hardware} shows our basic concept for the hardware
     881organization for the IPP.  This diagram shows the two types of compute
     882nodes: OTA-level processing and storage nodes (dominated by Phase 2)
     883and static sky processing and storage nodes (mostly Phase 4).  Also
     884shown are two switches used in this configuration; although it is
     885currently possible to buy a single switch with sufficient number of
     886ports, this organization represents a minimal configuration for the
     887PS-1 IPP hardware.  In such a case, the interswitch communication must
     888also meet the required throughput needs.  We discuss the hardware
     889requirements in the assumption that such an organization will be
     890necessary.
     891
     892The way in which the images are distributed among the storage and
     893compute nodes will largely determine the I/O bandwidth requirements.
     894For data bandwidth requirements calculations, it is necessary to make
     895some assumptions about the data organization.  We make the assumption
     896that the OTA data is optimally distributed to the OTA nodes such that
     897the OTA processing is always on a machine with local OTA data.  This
     898implies that all OTA data from a specific OTA are targetted to a
     899specific machine.  (see below for discussion of data duplication).
     900
     901A second factor which will have a significant impact on the I/O
     902requirements is the image storage format for the processed and
     903calibration images.  We have two basic choices: 32 bit floating point
     904format or 16 bit integer format with appropriate scaling.  In the
     905former case, additional dynamic range is retained, while in the latter
     906case, we reduce the data volume by a factor of 2.  Since the science
     907requirements for PS-1 do not specify a need for dynamic range greater
     908than 16 bits, we assume all images are stored as 16 bit data.
     909
     910A third determining factor is the number of calibration images needed
     911by the processing system.  Since the complete analysis is not yet
     912defined, this number is difficult to ascertain.  However, we can make
     913a reasonable guess at the total number for scaling purposes.  We
     914assume that each frame requires a total of 4 calibration frames on
     915average
     916
     917\begin{table}[b]
     918\begin{center}
     919\caption{Data Storage Requirements \label{storage}}
     920\begin{tabular}{lrrrr}
     921\hline
     922\hline
     923Raw data           & 200 TB \\
     924static sky         & 256 TB \\
     925calibration frames &   5 TB \\
     926metadata db        & 0.3 TB \\
     927object db          &   4 TB \\
     928\hline
     929total              & 116 TB \\
     930\hline
     931\end{tabular}
     932\end{center}
     933\end{table}
     934
     935\subsubsection{Data Storage Requirements}
     936
     937The Pan-STARRS IPP data storage requirements may be divided into five
     938principal areas: raw image data, static sky image data, master
     939calibration images, the metadata database, and the object database.
     940We discuss each of these data items and their impact on the data
     941storage requirements for the IPP for PS-1.  Table~\ref{storage}
     942summarizes the data storage requirements in the different scenarios.
     943
     944\paragraph{Raw Data Storage}
     945
     946There are two basic image types which will be acquired: night-time
     947science images and calibration images.  The night-time science images
     948consist of 1Gpix per image, or 2GB in raw format.  At nominal cadence,
     949the PS-1 telescope can obtain images at a sustained rate of 1 image
     950per 30 seconds for the entire night of 10 hours (36000 seconds).  A
     951total of 100 calibration images per night would be a substantial
     952overestimate of the typical expectation.  Combining these numbers, we
     953can expect to receive a total of 1300 images, or 2.6 TB of data per
     954night.  The total data storage requirements for the raw data are
     955governed by the number of nights' worth of data we are required to
     956keep online.  \tbd{for the first year, we are required to keep all
     957images from the AP and IPV surveys.  This amounts to a total of 200
     958TB of data}.
     959
     960\paragraph{Static Sky Data Storage}
     961
     962The static sky is represented by images with 0.2 arcsec per pixel.
     963There will be one summed image and one weight image for each of the
     964\tbd{6} filters, each stored with 16 bits of resolution, for a total
     965of 24 bytes per sky pixel.  At this resolution, there are 324 Mpix per
     966square degree, and we will observe a potential total area of 30,000
     967square degrees.  Allowing for 10\% overage for overlapping tiling, we
     968require a total of 10.7 Tpix to cover the sky once, or a total of
     969$\sim 256$ TB to maintain a single image of the static sky in all 6
     970filters.
     971
     972\paragraph{Calibration Frame Storage}
     973
     974The possible required calibration frames consist of the bias, dark,
     975and mask images, along with one flat, one flat-correction, and
     976multiple sky/fringe library frames per filter.  In fact, not all types
     977are needed at all stages.  It is very likely that we will not require
     978bias or dark images, and mask images may be represented by a single
     979byte per pixel.  Nonetheless, it is necessary for us to generate and
     980store all master calibration frames at least until we prove that they
     981are not needed.  We assume a total of 21 calibration images are
     982necessary (one flat, fringe, and sky per filter, along with a bias,
     983dark, and mask).  If we intend to keep all master calibration frames
     984for the project lifetime, and generate a new master on a weekly basis
     985(a reasonable time-scale), then we can expect to require a total of 5
     986TB of calibration image by the end of the 2 years of PS-1.  We note
     987that this is likely to be a drastic overestimate as we are unlikely to
     988need to regenerate all master calibration frames on a weekly
     989time-scale.
     990
     991\paragraph{Metadata Database Storage}
     992
     993The metadata data storage requirements are driven by the need to store
     994the data for the project lifetime.  There are two types of metadata
     995generated at the summit: data associated with images and environmental
     996data.  The environmental data consists of measurements on a regular
     997cadence, roughly 1 per minute, of a variety of parameters.  We suggest
     998an expected of 1kB per entry, for a total of 1 GB over the two-year
     999term of PS-1.  The additional systems, such as the DIMM, SkyProbe, NIR
     1000Sky Camera, and the LRProbe will have higher data requirements, but
     1001should be considered as separate, self-contained systems.  Their data
     1002products are distilled to a limited number of parameters per minute
     1003which are included in the 1kB given above.  Furthermore, items such as
     1004guide-star history, if saved, will be saved with the image data and
     1005represents only a small fraction of the total image data volume.  Some
     1006subset of the telescope diagnosic information may be a high volume
     1007data product as well, but only retained by the telescope control
     1008system for the purpose of diagnostic studies.  Such data will be
     1009excluded from this analysis.
     1010
     1011The image metadata consists of values associated with the FPA (1), the
     1012OTAs (64), and the Cells (4096).  Aside from the guide star history,
     1013the total data requirements for each of these entries will be scaled
     1014by the number of bytes required for the metadata from each data level.
     1015Clearly, if the Cell entry is allowed to be large, it will dominate
     1016the total Metadata data volume.  We suggest an expected number of 64
     1017bytes per Cell, 256 B per OTA, and 1k per FPA, yielding a total
     1018metadata volume per exposure of roughly 0.3 MB, completely dominated
     1019by the Cell metadata.  With the exposure rates above, we find a total
     1020of metadata volume of 0.3 TB over the two-year term of PS-1.
     1021
     1022\paragraph{Object Database Storage}
     1023
     1024The hardware requirements for the IPP object database are rather
     1025flexible: the total volume depends critically on the depth to which
     1026the object detection analyses are performed (and thus the total number
     1027of object detections) and the number of object parameters which are
     1028measured.  We can make very rough estimates that the total number of
     1029detections over the 2 year lifetime of the project may be in the
     1030vicinity of $10^{11}$.  We can conservatively estimate the number of
     1031bytes needed to represent each detection as 128 B, resulting in a
     1032total data storage for the object detections of 12 TB.  However, this
     1033number depends strongly on the timescale for which the IPP is required
     1034to maintain all object detections, and may potentially be
     1035significantly reduced.
     1036
     1037\subsubsection{CPU Requirements}
     1038
     1039Phase 2 and Phase 4 dominate the processing requirement, primarily
     1040because they must keep up with the image delivery rate of 1 per 30
     1041seconds.  We have performed benchmarks of a demonstration version for
     1042both the Phase 2 and Phase 4 analyses.
     1043
     1044For the Phase 2, a substantial fraction of the processing time is
     1045consumed by the need to perform FFTs on the images in order to
     1046convolve them with the guide-star kernel, and in the smoothing used
     1047for the object detection process.  Additional processing time is
     1048needed by the object detection, deblending, and analysis.  Experiments
     1049with the FFTW package show that FFTs may be performed on Intel
     1050processors at rates of approximately 0.25 GHz-sec / Mpix for data sets
     1051of order 1 Megapixel.  The FFTs required for the Phase 2 analysis are
     1052performed on the 512$^2$ pixel cells, so these numbers may roughly be
     1053scaled linearly to determine the total time required for OTA
     1054processing.  A single FFT on a full OTA, with 64 Cells, therefore
     1055requires roughly 4 GHz-sec.  For the full Phase 2 analysis, there are
     1056roughly 4 single direction FFTs required excluding those associated
     1057with object detection; thus the total processing time for these FFTs
     1058is approximately 16 GHz-sec.  The addtional analysis steps, excluding
     1059object detection and characterization, account for a small fraction of
     1060this compute time, which we estimate at 10\%.  The object detection
     1061stage depends somewhat on the depth to which the analysis is
     1062performed, and the number of measurements made per object.  Typical
     1063analysis performed by the Sextractor routine, which performs a
     1064substantial number of per-object analyses, requires 27 GHz-sec for a
     1065full OTA, including the FFTs used for smoothing.  We can therefore
     1066assume a total of 50 GHz-sec per OTA for the Phase 2 processing.  This
     1067converts to a total of 12800 GHz-sec for a complete major frame.
     1068
     1069For Phase 4, the main computational tasks are combining the multiple
     1070images, with cosmic-ray rejection, and performing the object detection
     1071tasks.  Nick Kaiser has done tests of the Phase 4 image combine and
     1072rejection stages, and finds a total processing time of roughly 96
     1073GHz-sec for a full stack of 4 OTA images.  If we add in an additional
     107434 GHz-sec for detailed object detection and image differencing, we
     1075find a conservative estimage of 130 GHz-sec for a 4-image OTA stack,
     1076equivalent to 7800 GHz-sec for a major frame.
     1077
     1078For PS-1, the typical time for a major frame is $4 \times 30$ seconds.
     1079Some reduction in the load may be gained by reducing the complexity
     1080and depth of analysis for PS-1.  Depending on the details and depth of
     1081the analysis, we may reduce the computational load by a factor of 2.
     1082
     1083\begin{table}
     1084\begin{center}
     1085\caption{Data I/O (MB per OTA or Sky-cell) \label{scenarios}}
     1086\begin{tabular}{lrrrr}
     1087\hline
     1088\hline
     1089{\em Phase 2 input}                                \\
     1090from summit    &                 $2 \times 32$ MB  \\
     1091input image    &                       {\bf 32 MB} \\
     1092calibration    &            {\bf 4 $\times$ 32 MB} \\
     1093mask image     &                       {\bf  8 MB} \\
     1094\hline
     1095network I/O:   &                            64 MB  \\
     1096disk I/O:      &                           176 MB  \\
     1097               &                                   \\
     1098{\em Phase 2 output}                               \\
     1099output image   &                      {\bf  32 MB} \\
     1100output mask    &                      {\bf   8 MB} \\
     1101image to P4    &               $1.5 \times 32$ MB  \\
     1102mask to P4     &               $1.5 \times  8$ MB  \\
     1103\hline
     1104network I/O:   &                            60 MB  \\
     1105disk I/O:      &                            40 MB  \\
     1106               &                                   \\
     1107{\em Phase 4}  &                                   \\
     1108input images   &      $1.5 \times 4 \times 32$ MB  \\
     1109input masks    &      $1.5 \times 4 \times  8$ MB  \\
     1110static sky     &                            32 MB  \\
     1111static weight  &                            32 MB  \\
     1112\hline
     1113input:         &                           304 MB  \\
     1114output:        &                            96 MB  \\
     1115\hline
     1116\multicolumn{5}{l}{\em Bold-faced entries are access to local-disk} \\
     1117\multicolumn{5}{l}{\em parenthesised disk I/O numbers are parallel with the network I/O} \\
     1118\end{tabular}
     1119\end{center}
     1120\end{table}
     1121
     1122\subsubsection{Per-Node I/O Requirements}
     1123
     1124Data I/O per node is defined as the number of bytes per second passed
     1125through the node's network adapter.  The data throughput for each node
     1126depends strongly on the how the data is organized and processed.  In
     1127this section, we identify the data which is passed between nodes for
     1128the two stages of the science analysis process.  Table~\ref{scenarios}
     1129lists the per-node data I/O for the analysis stages.
     1130
     1131For PS-1, there are 120 seconds of compute time allowed for each of
     1132the Phase 2 and Phase 4 analyses for the collection of four images
     1133which makes up a cannonical major frame.  We use the data I/O volumes
     1134and some assumptions about expected network and disk bandwidth to
     1135estimate the I/O and processing timeline for the four scenarios. From
     1136this analysis, we can judge the total CPU requirements in terms of
     1137GHz, not just GHz-sec.  We have assumed that GigE network adapters are
     1138capable of delivering data at 50MB/sec sustained and that a disk RAID
     1139can deliver sustained 100 MB/sec reads and writes.  These numbers are
     1140conservative estimates based on recent tests discussed below.  Using
     1141these assumptions, Table~\ref{throughput} lists the time allocations
     1142for the processing stages.
     1143
     1144\paragraph{Phase 2 Node I/O Requirements}
     1145
     1146In the assumed data distribution scenario, there is a single CPU
     1147allocated to each OTA in the OTA farm and a single CPU for each Sky
     1148cell process.  In addition, all data for the specified OTA are stored
     1149on local disks attached to the same computer as the CPU, with the
     1150result that all Phase 2 I/O is made to a local disk.  For each science
     1151OTA image which is observed, each OTA node will read from the network
     1152a total of 2 raw images (one for the original image, one for a backup
     1153copy) and write an average of roughly 1.5 processed images and masks
     1154to the Phase 4 machines for a total of 124 MB of network I/O.  During
     1155the processing stage, the OTA node will read from disk a total of 176
     1156MB (4 calibration frames at 32 MB each, one 16 MB mask, and one raw
     1157science image at 32 MB) and write a total of 40 MB (one processed
     1158image at 32 MB and one mask at 8 MB).  Given the assumptions for the
     1159network and disk bandwidths (50 MB/s and 100 MB/s respectively), the
     1160data volumes imply a total I/O period of 4.6 seconds.  In this
     1161instance, the network I/O is presumed to be sequential with the disk
     1162I/O.
     1163
     1164\paragraph{Phase 4 Node I/O Requirements}
     1165
     1166Although it is easy to arrange the OTA data in such a way that the
     1167majority of I/O is performed locally, it is not as easy to arrange
     1168this for the Static Sky data used by the Phase 4 analysis.  We
     1169therefore make the assumption that the Phase 4 analysis will require
     1170all input OTA data to be loaded across the network, as well as all
     1171Static Sky data.  This is somewhat of an overestimate as some of the
     1172Static Sky data will be processed by machines with the data stored
     1173locally, and clever Static-Sky data organization schemes can enhance
     1174this chance. 
     1175
     1176In the Phase 4 analysis, the images from the 4 separate telescopes are
     1177combined into a single image, confronted with the appropriate segment
     1178of the static sky, with output difference image and updated static sky
     1179image.  If we restrict input access to the individual OTA cells, the
     1180maximum read overhead is 50\% (need to read a 10x10 set of cells for
     1181an 8x8 input image).  If the processing is performed on Static Sky
     1182segments equivalent in size to the OTAs, the total volume of input
     1183data per node is 304 MB (192 MB of processed science image, 48 MB of
     1184input mask, 32 MB of static sky image and 32 MB of static sky weight
     1185map) while the output data is 96 MB (32 MB static sky, 32 MB weight
     1186map, and 32 MB difference image).  Thus, we require a total of 400 MB
     1187network I/O, which implies an I/O period of 8 seconds.
     1188
     1189\begin{table}
     1190\begin{center}
     1191\caption{Data Throughput \label{throughput}}
     1192\begin{tabular}{lrrrr}
     1193\hline
     1194\hline
     1195Phase 2 per-node network I/O       & 2.2 s           \\
     1196Phase 2 per-node disk I/O (read)   & 1.3 s           \\
     1197Phase 2 per-node disk I/O (write)  & 1.2 s           \\       
     1198Phase 2 CPU total                  &  25 s : 128 GHz \\
     1199Phase 4 per-node I/O               &   8 s           \\
     1200Phase 4 CPU total                  & 112 s : 70 GHz  \\
     1201Phase 2 switch load                & 264 MB/s \\
     1202Phase 4 switch load                & 215 MB/s \\
     1203Phase 2 to Phase 4 switch load     & 160 MB/s \\
     1204Summit to Phase 2 switch load      &  70 MB/s \\
     1205\hline
     1206\end{tabular}
     1207\end{center}
     1208\end{table}
     1209
     1210\subsubsection{Switch I/O Requirements}
     1211
     1212The switch I/O requirements are defined by the total number of bytes
     1213per second serviced by the two switches in the system. 
     1214
     1215The Phase 2 network I/O is 124 MB per OTA.  With 64 OTAs per image,
     1216and 30 seconds average between images, this implies a total of 264
     1217MB/s switch bandwidth.  The Phase 4 network I/O is 400 MB per sky
     1218cell.  With 64 cells and 120 seconds between major frames, this is an
     1219average switch bandwidth of 215 MB/s switch bandwidth.  The total
     1220switch-to-switch load is 304 MB per OTA, with an average timescale of
     1221120 seconds.  With 64 OTAs, this corresponds to 160 MB/s.  The
     1222summit-to-Phase 2 switch load is 70 MB/s.
     1223
     1224\begin{table}
     1225\begin{center}
     1226\caption{Hardware Throughput Tests \label{existing-hardware}}
     1227\begin{tabular}{lrrrr}
     1228\hline
     1229\hline
     1230Test        & where \& when     & model                & result                             \\
     1231\hline
     1232node I/O    & CFHT 11/2002      & Intel 1000 Gigabit   & 35 - 40 MB/s sustained             \\
     1233node I/O    & CFHT 2/2004       & Intel 1000 Gigabit   & 65 - 70 MB/s sustained             \\
     1234RAID write  & CFHT 2/2004       & 3ware RAID cntl + IDE & 110 MB/s sustained                 \\
     1235Switch Load & VeriTest          & Cisco                & 3 GB/s (for 32 ports)              \\
     1236\hline
     1237\end{tabular}
     1238\end{center}
     1239\end{table}
     1240
     1241\subsubsection{Existing Hardware Throughput}
     1242
     1243We have collected a few representative tests of various pieces of
     1244modern hardware to give a reference for the throughput capabilities.
     1245A number of hardware configurations have been tested at CFHT for the
     1246Elixir project, and we include here their recent reported hardware
     1247RAID-5 I/O speeds and GigE card speeds.  We also have included data
     1248from VeriTest studies of Cisco switch throughput, commissioned by
     1249Cisco for a 32 port GigE switch.  These tests are summarized in
     1250Table~\ref{existing-hardware}.
     1251
Note: See TracChangeset for help on using the changeset viewer.