Adding paper, paper build workflow.

.github/workflows/draft-pdf.yml

@@ -0,0 +1,24 @@
+on: [push]
+
+jobs:
+  paper:
+    runs-on: ubuntu-latest
+    name: Paper Draft
+    steps:
+      - name: Checkout
+        uses: actions/checkout@v4
+      - name: Build draft PDF
+        uses: openjournals/openjournals-draft-action@master
+        with:
+          journal: joss
+          # This should be the path to the paper within your repo.
+          paper-path: paper.md
+      - name: Upload
+        uses: actions/upload-artifact@v1
+        with:
+          name: paper
+          # This is the output path where Pandoc will write the compiled
+          # PDF. Note, this should be the same directory as the input
+          # paper.md
+          path: paper.pdf
+

paper.bib

@@ -0,0 +1,59 @@
+@article{Treeby2010,
+  title = {k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields},
+  volume = {15},
+  ISSN = {1083-3668},
+  url = {http://dx.doi.org/10.1117/1.3360308},
+  DOI = {10.1117/1.3360308},
+  number = {2},
+  journal = {Journal of Biomedical Optics},
+  publisher = {SPIE-Intl Soc Optical Eng},
+  author = {Treeby,  Bradley E. and Cox,  B. T.},
+  year = {2010},
+  pages = {021314}
+}
+
+@INPROCEEDINGS{Garcia2021,
+  author={Garcia, Damien},
+  booktitle={2021 IEEE International Ultrasonics Symposium (IUS)}, 
+  title={Make the most of MUST, an open-source Matlab UltraSound Toolbox}, 
+  year={2021},
+  volume={},
+  number={},
+  pages={1-4},
+  keywords={Analytical models;Adaptation models;Ultrasonic imaging;Transducers;Array signal processing;Speckle;Acoustics;Ultrasonic imaging;Medical simulation;Signal processing;Image processing;Open access},
+  doi={10.1109/IUS52206.2021.9593605}}
+
+@INPROCEEDINGS{USTB2017,
+  author={Rodriguez-Molares, Alfonso and Rindal, Ole Marius Hoel and Bernard, Olivier and Nair, Arun and Lediju Bell, Muyinatu A. and Liebgott, Hervé and Austeng, Andreas and L⊘vstakken, Lasse},
+  booktitle={2017 IEEE International Ultrasonics Symposium (IUS)}, 
+  title={The UltraSound ToolBox}, 
+  year={2017},
+  volume={},
+  number={},
+  pages={1-4},
+  keywords={Ultrasonic imaging;Pipelines;Imaging;Array signal processing;Program processors;Signal processing algorithms;MATLAB;Beamforming;signal processing;comparison;toolbox},
+  doi={10.1109/ULTSYM.2017.8092389}}
+
+@inproceedings{heinisch2018MatCL,
+  title={{MatCL}: {A} new easy-to use {OpenCL} toolbox for {MathWorks} {Matlab}},
+  author={Heinisch, Philip and Ostaszewski, Katharina},
+  year={2018},
+  pages={8:1--8:1},
+  booktitle={Proceedings of the International Workshop on OpenCL},
+  series={IWOCL '18, May 2018, Oxford (United Kingdom)},
+  publisher={ACM},
+  address={New York, NY, USA},
+  note={\url{https://github.com/IANW-Projects/MatCL}},
+  doi={10.1145/3204919.3204927}
+}
+
+@article{Jensen1996,
+author = {Jensen, Jørgen},
+year = {1996},
+month = {01},
+pages = {351-352},
+title = {FIELD: A program for simulating ultrasound systems},
+volume = {34},
+journal = {Medical and Biological Engineering and Computing}
+}
+

paper.md

@@ -0,0 +1,66 @@
+---
+Title: 'QUPS: A MATLAB Toolbox for Ultrasound Imaging Prototyping'
+tags:
+  - MATLAB
+  - ultrasound
+  - beamforming
+  - imaging
+  - prototyping
+  - 'signal processing'
+  - CUDA
+  - OpenCL
+
+authors:
+  - name:
+    orcid: 0000-0002-0720-6955
+    first: Thurston 
+    surname: Brevett
+    affiliation: 1
+affiliations:
+  - name: Stanford University
+    index: 1
+bibliography: paper.bib
+---
+
+# Summary:
+
+Medical ultrasound encompasses a rapidly developing field of research with many potential applications from diagnostics to therapeutics. Public data initiatives
+and recent industry trends have reduced the cost and complexity to engage in research via publicly available data, on-board data acquisition functions, and commercial off the shelf (COTS) imaging systems. However, terminology, data formatting standards, and computing practices tend to be highly fragmented throughout the ultrasound research field, leading to a high barrier of entry that hinders communication, code reuse, and shareability. 
+
+`QUPS` aims to facilitate ultrasound research by offering a flexible, performant, modular, and intuitive toolbox that can be used by experienced researchers as well as curious students for data acquisition, generation, representation, and processing. It pursues this goal by offering tools designed to reduce the burden of common computational tasks involved in the scientific process by standardizing data formats and offering hardware accelerated functionality.
+
+# Statement of need:
+
+The lack of data and terminology standardization in ultrasound technology research is an easily overlooked impediment to developing research, as it encourages reinventing the wheel and can easily lead to misinterpretation.
+At its lowest level, `QUPS` provides memory efficient definition classes capable of representing a wide variety of common imaging systems and enforces standardized units, coordinate systems, and nomenclature. 
+All classes have data visualization support, making it far easier to share and validate data amongst collaborators and with the public. 
+
+`QUPS` provides a variety of compute and memory optimized implementations of common signal processing methods (e.g. beamforming, resampling, or convolution), which take advantage of native parallelization tools within MATLAB or hardware acceleration via CUDA or OpenCL when feasible by leveraging `MatCL`[@heinisch2018MatCL]. `QUPS` exploits data parallellism primarily using a single program multiple data (SPMD) approach, and optimistically dispatches to the most efficient implementation available on a user's computer. For access to cutting edge research, it can also export or import to `USTB`[@USTB2017], an Ultrasound toolbox which contains a library of advanced algorithms from the literature.
+
+To facilitate data acquisition, `QUPS` provides tools to import data from [Verasonics Vantage](https://verasonics.com/vantage-systems/) ultrasound research scanners, which uses a flexible yet highly complex data standard. Similarly, to facilitate data generation, `QUPS` also interfaces with multiple popular free ultrasound field simulators including `FieldII`[@Jensen1996], `MUST`[@Garcia2021], and `k-Wave`[@Treeby2010]. Simulations can be batched on a compute cluster via MATLAB's parallel server toolbox.
+
+Alternatively, when a coarse but fast approximation is preferable, `QUPS` offers a hardware accelerated simulator with minimal complexity. This method uses a green's function to simulate ideal isotropic point sources (transmitters), point receivers, and point scatterers with radial propagation loss using the equation
+[$$ x_{mn}(t) = \sum_s \frac{v(t - c_0^{-1} \left( ||\mathbf{p}_s - \mathbf{q}_m|| + ||\mathbf{p}_s - \mathbf{q}_n|| \right) )}{||\mathbf{p}_s - \mathbf{q}_m|| \cdot ||\mathbf{p}_s - \mathbf{q}_n||} $$]{label="eq:greens"}
+where $v(\cdot)$ is the transmit waveform, $m$, $n$, and $s$ are the transmit, receive, and scatterer indices, $\mathbf{q}_m$ and $\mathbf{q}_n$ $\mathbf{p}_s$ is the location of scatterer $s$, and $x_{mn}(\cdot)$ is the received echo waveform as a function of time.
+\autoref{eq:greens} is embarassingly parallel over $m$ and $n$, and can be further accelerated with a branch-and-bound approach used to minimize unnecessary computation. Both of these aspects are exploited in native MATLAB, as well as in CUDA and OpenCL kernels.
+
+Broadcasting conventions are employed throughout to reduce the memory footprint while offering significant flexibility. 
+For example, for delay-and-sum beamforming, one of the most common operations in ultrasound imaging, `QUPS` provides a GPU-enabled method employing the equation
+$$ b_{ijkf} = \sum_m \sum_n \alpha_{ijknm} \cdot x_{mnf}(||\mathbf{p}_{ijk} - \mathbf{q}_n|| + (-1)^{s(\mathbf{p}_{ijk}, \mathbf{q}_m)} \cdot ||\mathbf{p}_{ijk} - \mathbf{q}_m|| - t_m) $$
+where $m$, $n$, and $f$ are the transmit, receive and frame indices, $i$, $j$, and $k$ are the voxel (or pixel) indices, $x_{mnf}(\cdot)$ is the voltage signal (trace) as a function of time, $t_m$ is the start time for transmit $m$, $\mathbf{p}_{ijk}$ are the voxel positions in 3D cartesian coordinates, $\mathbf{q}_m$ and $\mathbf{q}_n$ are the virtual transmit and receiver positions, and $s(\cdot,\cdot)$ is a threshold function based on the relation between the virtual transmit position and the imaging voxel position. 
+Rather than require the user to specify the five dimensional tensor $\alpha_{ijknm}$, any subset of these dimensions may be singular to efficiently form e.g. a transmit by lateral weighting (apodization) scheme or a receiver by depth weighting scheme.
+Permutations of the input data e.g. $x_{nfm}$ or $\mathbf{p}_{kji}$ are handled via properties specifying the indices rather than by enforcing a convention, which may lead to costly memory operations.
+
+Engaging in research often involves a plethora of small choices in order to configure software to suit the needs of a particular scientific inquiry. While exposing the user is necessary for some fields of research, this forces all users to consider each aspect of the methodology. 
+In most cases, the vast majority of these choices are either ultimately inconsequential or there exists a clear optimal or at a minumum, a clear _acceptable_ choice.
+For example, for finite-difference time-domain (FDTD) simulators, knowledge of a transducer's frequency bandwidth is sufficient to determine a reasonably sized temporal sampling interval that satisfies the Courant-Friedrichs-Lewy (CFL) stability criterion.
+For research regarding the stability and accuracy of such simulators, the sampling interval is a critical choice. In many other cases however, this aspect is irrelevant as long as the sampling interval is _sufficiently_ small.
+Requiring a user to be sufficiently well-versed creates a barrier to entry with little benefit.
+
+Throughout the toolbox, `QUPS` provides reasonable defaults that abstract and simplify the number and complexity of such choices so that researchers can focus on the big picture questions.
+When these choices do become critical as is the case for researchers and experts, such users can engage with the toolbox at a lower level, or simply use the available functionality as a template to be modified to suit their needs.
+
+# Acknowledgements
+
+# References
+