Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matlab data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-E-J-117: Verification Method for the Detection Accuracy of Automatic Winston Lutz Test

Purpose: Winston Lutz test (WLT) has been a standard QA procedure performed prior to SRS treatment, to verify the mechanical iso-center setup accuracy upon different Gantry/Couch movements. Several detection algorithms exist,for analyzing the ball-radiation field alignment automatically. However, the accuracy of these algorithms have not been fully addressed. Here, we reveal the possible errors arise from each step in WLT, and verify the software detection accuracy with the Rectilinear Phantom Pointer (RLPP), a tool commonly used for aligning treatment plan coordinate with mechanical iso-center. Methods: WLT was performed with the radio-opaque ball mounted on a MIS and irradiated onto EDR2 films. The films were scanned and processed with an in-house Matlab program for automatic iso-center detection. Tests were also performed to identify the errors arise from setup, film development and scanning process. The radioopaque ball was then mounted onto the RLPP, and offset laterally and longitudinally in 7 known positions ( 0, ±0.2, ±0.5, ±0.8 mm) manually for irradiations. The gantry and couch was set to zero degree for all irradiation. The same scanned images were processed repeatedly to check the repeatability of the software. Results: Miminal discrepancies (mean=0.05mm) were detected with 2 films overlapped and irradiated but developed separately. This reveals the error arise from film processor and scanner alone. Maximum setup errors were found to be around 0.2mm, by analyzing data collected from 10 irradiations over 2 months. For the known shift introduced using the RLPP, the results agree with the manual offset, and fit linearly (R{sup 2}>0.99) when plotted relative to the first ball with zero shift. Conclusion: We systematically reveal the possible errors arise from each step in WLT, and introduce a simple method to verify the detection accuracy of our in-house software using a clinically available tool.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Modeling and Control System Design for an Integrated Solar Generation and Energy Storage System with a Ride-Through Capability

This paper presents a generic approach for PV panel modeling. Data for this modeling can be easily obtained from manufacturer datasheet, which provides a convenient way for the researchers and engineers to investigate the PV integration issues. A two-stage power conversion system (PCS) is adopted in this paper for the PV generation system and a Battery Energy Storage System (BESS) can be connected to the dc-link through a bi-directional dc/dc converter. In this way, the BESS can provide some ancillary services which may be required in the high penetration PV generation scenario. In this paper, the fault ride-through (FRT) capability is specifically focused. The integrated BESS and PV generation system together with the associated control systems is modeled in PSCAD and Matlab platforms and the effectiveness of the controller is validated by the simulation results.

14 SOLAR ENERGY↗
DOE OSTI2023

Soil temperature and soil moisture raw data, permafrost table depths, and accompanying environmental variable data, Kenai Wildlife Refuge, 2019-2022

Data package purpose: This data package was created to contain all data used in an upcoming article, "Canopy Cover and Microtopography Control Precipitation-Enhanced Thaw of Ecosystem-Protected Permafrost." In review.This data package includes: Raw output from 19 distributed temperature profilers with a thermistor every 10 cm along a 160 cm length at a measurement interval of 15 minutes (.CSV). Raw output from two soil moisture and temperature profilers (90 cm length and 120 cm length) that took composite soil moisture readings every 15 cm along the sensor length at a measurement interval of 30 minutes (.CSV). Permafrost depths were measured annually in mid-September at DTP sensor locations (.CSV) and along an across-site transect (.CSV). Environmental variables (snow depth, canopy closure, moss depth, and elevation) for all sensor locations. Real-time kinetic (RTK) GPS points showing site microtopography (.CSV).Analysis software: Our analysis was done in Matlab. File types can be used with any software.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

TH-C-BRD-05: Reducing Proton Beam Range Uncertainty with Patient-Specific CT HU to RSP Calibrations Based On Single-Detector Proton Radiography

Purpose: Beam range uncertainty in proton treatment comes primarily from converting the patient's X-ray CT (xCT) dataset to relative stopping power (RSP). Current practices use a single curve for this conversion, produced by a stoichiometric calibration based on tissue composition data for average, healthy, adult humans, but not for the individual in question. Proton radiographs produce water-equivalent path length (WEPL) maps, dependent on the RSP of tissues within the specific patient. This work investigates the use of such WEPL maps to optimize patient-specific calibration curves for reducing beam range uncertainty. Methods: The optimization procedure works on the principle of minimizing the difference between the known WEPL map, obtained from a proton radiograph, and a digitally-reconstructed WEPL map (DRWM) through an RSP dataset, by altering the calibration curve that is used to convert the xCT into an RSP dataset. DRWMs were produced with Plastimatch, an in-house developed software, and an optimization procedure was implemented in Matlab. Tests were made on a range of systems including simulated datasets with computed WEPL maps and phantoms (anthropomorphic and real biological tissue) with WEPL maps measured by single detector proton radiography. Results: For the simulated datasets, the optimizer showed excellent results. It was able to either completely eradicate or significantly reduce the root-mean-square-error (RMSE) in the WEPL for the homogeneous phantoms (to zero for individual materials or from 1.5% to 0.2% for the simultaneous optimization of multiple materials). For the heterogeneous phantom the RMSE was reduced from 1.9% to 0.3%. Conclusion: An optimization procedure has been designed to produce patient-specific calibration curves. Test results on a range of systems with different complexities and sizes have been promising for accurate beam range control in patients. This project was funded equally by the Engineering and Physical Sciences Research Council (UK) and Ion Beam Applications (Louvain-La-Neuve, Belgium)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-D-201-06: Random Walk Algorithm Seed Localization Parameters in Lung Positron Emission Tomography (PET) Images

Purpose: The objective of this study was to find the best seed localization parameters in random walk algorithm application to lung tumor delineation in Positron Emission Tomography (PET) images. Methods: PET images suffer from statistical noise and therefore tumor delineation in these images is a challenging task. Random walk algorithm, a graph based image segmentation technique, has reliable image noise robustness. Also its fast computation and fast editing characteristics make it powerful for clinical purposes. We implemented the random walk algorithm using MATLAB codes. The validation and verification of the algorithm have been done by 4D-NCAT phantom with spherical lung lesions in different diameters from 20 to 90 mm (with incremental steps of 10 mm) and different tumor to background ratios of 4:1 and 8:1. STIR (Software for Tomographic Image Reconstruction) has been applied to reconstruct the phantom PET images with different pixel sizes of 2×2×2 and 4×4×4 mm{sup 3}. For seed localization, we selected pixels with different maximum Standardized Uptake Value (SUVmax) percentages, at least (70%, 80%, 90% and 100%) SUVmax for foreground seeds and up to (20% to 55%, 5% increment) SUVmax for background seeds. Also, for investigation of algorithm performance on clinical data, 19 patients with lung tumor were studied. The resulted contours from algorithm have been compared with nuclear medicine expert manual contouring as ground truth. Results: Phantom and clinical lesion segmentation have shown that the best segmentation results obtained by selecting the pixels with at least 70% SUVmax as foreground seeds and pixels up to 30% SUVmax as background seeds respectively. The mean Dice Similarity Coefficient of 94% ± 5% (83% ± 6%) and mean Hausdorff Distance of 1 (2) pixels have been obtained for phantom (clinical) study. Conclusion: The accurate results of random walk algorithm in PET image segmentation assure its application for radiation treatment planning and diagnosis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

NWEI Azura RTI 1/20th Model Validation Wave Tank Test Data

Data from the 1/20th wave tank test of the RTI model. NWEI has licensed intellectual property from RTI, and modified the PTO and retested the 1/20th RTI model that was tested as part of the Wave Energy Prize. The goal of the test was to validate NWEI's simulation models of the model. The test occurred at the University of Maine in Orono.

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

Neutron Detector Signal Processing to Calculate the Effective Neutron Multiplication Factor of Subcritical Assemblies

This report describes different methodologies to calculate the effective neutron multiplication factor of subcritical assemblies by processing the neutron detector signals using MATLAB scripts. The subcritical assembly can be driven either by a spontaneous fission neutron source (e.g. californium) or by a neutron source generated from the interactions of accelerated particles with target materials. In the latter case, when the particle accelerator operates in a pulsed mode, the signals are typically stored into two files. One file contains the time when neutron reactions occur and the other contains the times when the neutron pulses start. In both files, the time is given by an integer representing the number of time bins since the start of the counting. These signal files are used to construct the neutron count distribution from a single neutron pulse. The built-in functions of MATLAB are used to calculate the effective neutron multiplication factor through the application of the prompt decay fitting or the area method to the neutron count distribution. If the subcritical assembly is driven by a spontaneous fission neutron source, then the effective multiplication factor can be evaluated either using the prompt neutron decay constant obtained from Rossi or Feynman distributions or the Modified Source Multiplication (MSM) method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

SU-E-J-240: The Impact On Clinical Dose-Distributions When Using MR-Images Registered with Stereotactic CT-Images in Gamma Knife Radiosurgery

Purpose: Most Gamma knife treatments are based solely on MR-images. However, for fractionated treatments and to implement TPS dose calculations that require electron densities, CT image data is essential. The purpose of this work is to assess the dosimetric effects of using MR-images registered with stereotactic CT-images in Gamma knife treatments. Methods: Twelve patients treated for vestibular schwannoma with Gamma Knife Perfexion (Elekta Instruments, Sweden) were selected for this study. The prescribed doses (12 Gy to periphery) were delivered based on the conventional approach of using stereotactic MR-images only. These plans were imported into stereotactic CT-images (by registering MR-images with stereotactic CT-images using the Leksell gamma plan registration software). The dose plans, for each patient, are identical in both cases except for potential rotations and translations resulting from the registration. The impact of the registrations was assessed by an algorithm written in Matlab. The algorithm compares the dose-distributions voxel-by-voxel between the two plans, calculates the full dose coverage of the target (treated in the conventional approach) achieved by the CT-based plan, and calculates the minimum dose delivered to the target (treated in the conventional approach) achieved by the CT-based plan. Results: The mean dose difference between the plans was 0.2 Gy to 0.4 Gy (max 4.5 Gy) whereas between 89% and 97% of the target (treated in the conventional approach) received the prescribed dose, by the CT-plan. The minimum dose to the target (treated in the conventional approach) given by the CT-based plan was between 7.9 Gy and 10.7 Gy (compared to 12 Gy in the conventional treatment). Conclusion: The impact of using MR-images registered with stereotactic CT-images has successfully been compared to conventionally delivered dose plans showing significant differences between the two. Although CTimages have been implemented clinically; the effect of the registration has not been fully investigated.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-241: Creation of Ventilation CT From Daily 4D CTs Or 4D Conebeam CTs Acquired During IGRT for Thoracic Cancers

Purpose: To develop a method to create ventilation CTs from daily 4D CTs or 4D KV conebeam CTs (4DCBCT) acquired during image-guided radiation therapy (IGRT) for thoracic tumors, and to explore the potential for using the ventilation CTs as a means for early detection of lung injury during radiation treatment. Methods: 4DCT acquired using an in-room CT (CTVision, Siemens) and 4DCBCT acquired using the X-ray Volume Imaging (XVI) system (Infinity, Elekta) for representative lung cancer patients were analyzed. These 4D data sets were sorted into 10 phase images. A newly-available deformable image registration tool (ADMIRE, Elekta) is used to deform the phase images at the end of exhale (EE) to the phase images at the end of inhale (EI). The lung volumes at EI and EE were carefully contoured using an intensity-based auto-contour tool and then manually edited. The ventilation images were calculated from the variations of CT numbers of those voxels masked by the lung contour at EI between the registered phase images. The deformable image registration is also performed between the daily 4D images and planning 4DCT, and the resulting deformable field vector (DFV) is used to deform the planning doses to the daily images by an in-house Matlab program. Results: The ventilation images were successfully created. The tide volumes calculated using the ventilation images agree with those measured through volume difference of contours at EE and EI, indicating the accuracy of ventilation images. The association between the delivered doses and the change of lung ventilation from the daily ventilation CTs is identified. Conclusions: A method to create the ventilation CT using daily 4DCTs or 4D KV conebeam CTs was developed and demonstrated.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Integrate Latimer Controls' Solution into RTAC (CRADA Final Report, CRD-23-24672)

Latimer Controls, Inc. was awarded two vouchers under the Department of Energy's American-Made Solar Prize Round 6 to conduct collaborative research at a national laboratory. The National Renewable Energy Laboratory (NREL) was selected as a partner to assist Latimer Controls in the performance evaluation of its photovoltaic (PV) control software. This collaboration focuses on developing a hardware-in-the-loop (HIL) testbed at NREL, which will be used to test and validate the Latimer PV control technology in a realistic yet de-risked environment. Both Latimer and NREL teams will work together to analyze the collected test data, derive insights, and disseminate the scientific findings. Recent studies underscore the potential of solar energy as a zero-marginal-cost and zero-emission flexibility resource within the bulk power system, particularly when integrated with advanced control systems. To enhance the performance of such systems, Latimer Controls has developed leading-edge technologies, including machine learning (ML) algorithms and hierarchical inverter set-point allocation methods. These innovations are designed to estimate the operational headroom of large PV plants for grid integration and control. However, comprehensive validation under real-world conditions remains necessary. To address this gap, the concurrent CRADA project proposes the real-world application and validation of the Latimer Control solution within a HIL environment. Initially, the Latimer algorithm was developed and tested within MATLAB Simulink, a platform suitable for research-level simulations and iterative development. However, transitioning this technology to a real solar site as an industry-ready solution necessitates implementation in a format compatible with widely used solar power plant controllers. In this additional CRADA work, the MATLAB Simulink-based logic will be translated into Structured Text, a programming language compliant with IEC 61131 standards, which is commonly used for custom logic implementations in industry-leading programmable logic controllers (PLCs), such as the Schweitzer SEL real-time automation controller (RTAC). This transition will facilitate the deployment of the Latimer Control solution in real-world solar power plants, thereby advancing the technology towards commercialization.

14 SOLAR ENERGY↗
DOE OSTI2012

A novel conformity index for intensity modulated radiation therapy plan evaluation

Purpose: Intensity modulated radiation therapy (IMRT) has gained popularity in the treatment of cancers. Manual evaluation of IMRT plans for head-and-neck cancers has been especially challenging necessitating efficient and objective assessment tools. In this work, the authors address this issue by developing a personalized conformity index (CI) for comparison of IMRT plans for head-and-neck cancers and evaluating its plan quality discerning power in comparison with other widely used CIs. Methods: A two-dimensional CI with dose and distance incorporated (CI{sub DD}) was developed using the MATLAB program language, to quantify the planning target volume (PTV) coverage. Valuable information contained in the digital imaging and communication in medicine (DICOM) RT objects were harvested for computation of each of the CI{sub DD} components. Apart from the dose penalty factor, a distance-based exponential function was employed by varying the penalty weight associated with the location of cold spots within the PTV. With the goal of deriving a customized penalty factor, the distances between individual pixel and its nearest PTV boundary was found. Using the exponential function, the impact of distance penalty was substantially larger for cold spots closer to the PTV centroid but petered out quickly wherever they were situated in the vicinity of PTV border. In order to evaluate the CI{sub DD} scoring system, three CT image data sets of nasopharyngeal carcinoma (NPC) patients were collected. Ten IMRT plans with degrading qualities were generated from each dataset and were ranked based on CI{sub DD} and other existing indices. The coefficient of variance was calculated for each dataset to compare the degree of variation. Results: The CI{sub DD} scoring system that considered spatial importance of each voxel within the PTV was successfully developed. The results demonstrated that the CI{sub DD} including four discrete factors could provide accurate rankings of plan quality by examining the relative importance of each cold spot within the PTVs. Apart from the dose penalty factor, a distance-based exponential function was employed taking the specific tumor geometry into account. Compared with other commonly used CIs, the CI{sub DD} resulted in the largest coefficient of variance among the ten IMRT plans for each dataset, indicating that its discerning power was the best among the CIs being compared. Conclusions: The CI{sub DD} scoring system was successfully developed to incorporate patient-specific spatial dose information and provide a geometry-based physical index for comparison of IMRT plans for head-and-neck cancers. By taking individual tumor geometry into account, the superiority of CI{sub DD} in plan discerning power was demonstrated. The use of CI{sub DD} could provide an effective means of benchmarking performance, reducing treatment plan variability, and advancing the quality of current IMRT planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-BRA-02: Development of a Learning Based Block Matching Algorithm for Ultrasound Tracking in Radiotherapy

Purpose: To develop an ultrasound learning-based tracking algorithm with the potential to provide real-time motion traces of anatomy-based fiducials that may aid in the effective delivery of external beam radiation. Methods: The algorithm was developed in Matlab R2015a and consists of two main stages: reference frame selection, and localized block matching. Immediately following frame acquisition, a normalized cross-correlation (NCC) similarity metric is used to determine a reference frame most similar to the current frame from a series of training set images that were acquired during a pretreatment scan. Segmented features in the reference frame provide the basis for the localized block matching to determine the feature locations in the current frame. The boundary points of the reference frame segmentation are used as the initial locations for the block matching and NCC is used to find the most similar block in the current frame. The best matched block locations in the current frame comprise the updated feature boundary. The algorithm was tested using five features from two sets of ultrasound patient data obtained from MICCAI 2014 CLUST. Due to the lack of a training set associated with the image sequences, the first 200 frames of the image sets were considered a valid training set for preliminary testing, and tracking was performed over the remaining frames. Results: Tracking of the five vessel features resulted in an average tracking error of 1.21 mm relative to predefined annotations. The average analysis rate was 15.7 FPS with analysis for one of the two patients reaching real-time speeds. Computations were performed on an i5-3230M at 2.60 GHz. Conclusion: Preliminary tests show tracking errors comparable with similar algorithms at close to real-time speeds. Extension of the work onto a GPU platform has the potential to achieve real-time performance, making tracking for therapy applications a feasible option. This work is partially funded by NIH grant R01CA190298.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Bayesian multiple-instance motif discovery with BAMBI: inference of recombinase and transcription factor binding sites

Finding conserved motifs in genomic sequences represents one of essential bioinformatic problems. However, achieving high discovery performance without imposing substantial auxiliary constraints on possible motif features remains a key algorithmic challenge. This work describes BAMBI—a sequential Monte Carlo motifidentification algorithm, which is based on a position weight matrix model that does not require additional constraints and is able to estimate such motif properties as length, logo, number of instances and their locations solely on the basis of primary nucleotide sequence data. Furthermore, should biologically meaningful information about motif attributes be available, BAMBI takes advantage of this knowledge to further refine the discovery results. In practical applications, we show that the proposed approach can be used to find sites of such diverse DNA-binding molecules as the cAMP receptor protein (CRP) and Din-family site-specific serine recombinases. Results obtained by BAMBI in these and other settings demonstrate better statistical performance than any of the four widely-used profile-based motif discovery methods: MEME, BioProspector with BioOptimizer, SeSiMCMC and Motif Sampler as measured by the nucleotide-level correlation coefficient. Additionally, in the case of Din-family recombinase target site discovery, the BAMBI-inferred motif is found to be the only one functionally accurate from the underlying biochemical mechanism standpoint. C++and Matlab code is available at http://www.ee.columbia.edu/guido/ BAMBI or http://genomics.lbl.gov/BAMBI/.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Integration of a DER Management System in Riverside. Final report

The tasks in this project covered various aspects, including algorithm development, algorithm integration into a commercial Active Network Management (ANM) platform, hardware-in-the-loop (HIL) testing in an industry-standard testing platform, pilot demonstration in Riverside, California, and also cost and benefit analysis. The DERMS platform in this project can host different algorithms developed on different platforms (e.g., MATLAB and Python) and it can interact with different hardware devices (e.g., different PV inverters, battery inverters, and different sensors). The DER control solution are based on an advanced model-free, layered, and clustered DER control paradigm. At the core of the DER control algorithms was the concept of Extremum Seeking (ES), which is a model-free probing-based control technique. The ES-based control algorithms were tested on major real-world inverters; both individually and in a cluster. It was shown that even legacy equipment (or when paired with a few additional advanced equipment) can support such advanced control. The monitoring algorithms utilize a heterogeneous set of legacy and advanced sensor measurements, such as behind-the-meter DER sensors, distribution-level Phase Measurement Units, distribution-substation Supervisory Control and Data Acquisition (SCADA), and line current sensors, with their limited availability; in order to infer practical network conditions. Sensor data are utilized to achieve resource forecasting, phase identification, and distribution system state estimation. The technology that was developed and demonstrated in this project could be transformational to utilities, including the smaller municipal utilities such as in Riverside, which may not have the resources to deploy advanced distribution system and DERMS solutions in order to support high penetration of solar power integration. This project created a real-world prototype to provide utilities with an assessment of smart grid monitoring and control technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2013

The need for application-based adaptation of deformable image registration

Purpose: To utilize a deformable phantom to objectively evaluate the accuracy of 11 different deformable image registration (DIR) algorithms. Methods: The phantom represents an axial plane of the pelvic anatomy. Urethane plastic serves as the bony anatomy and urethane rubber with three levels of Hounsfield units (HU) is used to represent fat and organs, including the prostate. A plastic insert is placed into the phantom to simulate bladder filling. Nonradiopaque markers reside on the phantom surface. Optical camera images of these markers are used to measure the positions and determine the deformation from the bladder insert. Eleven different DIR algorithms are applied to the full and empty-bladder computed tomography images of the phantom (fixed and moving volumes, respectively) to calculate the deformation. The algorithms include those from MIM Software (MIM) and Velocity Medical Solutions (VEL) and nine different implementations from the deformable image registration and adaptive radiotherapy toolbox for Matlab. These algorithms warp one image to make it similar to another, but must utilize a method for regularization to avoid physically unrealistic deformation scenarios. The mean absolute difference (MAD) between the HUs at the marker locations on one image and the calculated location on the other serves as a metric to evaluate the balance between image similarity and regularization. To demonstrate the effect of regularization on registration accuracy, an additional beta version of MIM was created with a variable smoothness factor that controls the emphasis of the algorithm on regularization. The distance to agreement between the measured and calculated marker deformations is used to compare the overall spatial accuracy of the DIR algorithms. This overall spatial accuracy is also utilized to evaluate the phantom geometry and the ability of the phantom soft-tissue heterogeneity to represent patient data. To evaluate the ability of the DIR algorithms to accurately transfer anatomical contours, the rectum is delineated on both the fixed and moving images. A Dice similarity coefficient is then calculated between the contour on the fixed image and that transferred, via the calculated deformation, from the moving to the fixed image. Results: The phantom possesses sufficient soft-tissue heterogeneity to act as a proxy for patient data. Large discrepancies appear between the algorithms and the measured ground-truth deformation. VEL yields the smallest mean spatial error and a Dice coefficient of 0.90. MIM produces the lowest MAD value and the highest Dice coefficient of 0.96, but creates the largest spatial errors. Increasing the MIM smoothness factor above the default value improves the overall spatial accuracy, but the factor associated with the lowest mean error decreases the Dice coefficient to 0.85. Conclusions: Different applications of DIR require disparate balances between image similarity and regularization. A DIR algorithm that is optimized only for its ability to transfer anatomical contours will yield large deformation errors in homogeneous regions, which is problematic for dose mapping. For this reason, these algorithms must be tested for their overall spatial accuracy. The developed phantom is an objective tool for this purpose.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

SU-E-T-545: MLC Distance Travelled as a Predictor for Motor Failure

Purpose: To study the frequency of Multi-Leaf Collimator (MLC) leaf failures, investigate methods to predict them and reduce linac downtime. Methods: A Varian HD120 MLC was used in our study. The hyperterminal MLC errors logged from 06/2012 to 12/2014 were collected. Along with the hyperterminal errors, the MLC motor changes and all other MLC interventions by the linear accelerator engineer were recorded. The MLC dynalog files were also recorded on a daily basis for each treatment and during linac QA. The dynalog files were analyzed to calculate root mean square errors (RMS) and cumulative MLC travel distance per motor. An in-house MatLab code was used to analyze all dynalog files, record RMS errors and calculate the distance each MLC traveled per day. Results: A total of 269 interventions were recorded over a period of 18 months. Of these, 146 included MLC motor leaf change, 39 T-nut replacements, and 84 MLC cleaning sessions. Leaves close to the middle of each side required the most maintenance. In the A bank, leaves A27 to A40 recorded 73% of all interventions, while the same leaves in the B bank counted for 52% of the interventions. On average, leaves in the middle of the bank had their motors changed approximately every 1500m of travel. Finally, it was found that the number of RMS errors increased prior to an MLC motor change. Conclusion: An MLC dynalog file analysis software was developed that can be used to log daily MLC usage. Our eighteen-month data analysis showed that there is a correlation between the distance an MLC travels, the RMS and the life of the MLC motor. We plan to use this tool to predict MLC motor failures and with proper and timely intervention, reduce the downtime of the linac during clinical hours.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

SU-E-T-679: Retrospective Analysis of the Sensitivity of Planar Dose Measurements To Gamma Analysis Criteria

Purpose: IMRT QA using planar dose measurements is still a widely used method for checking the accuracy of treatment plans. A pass/fail judgment is made using gamma analysis based on a single endpoint. Using more stringent criteria is a way to increase the sensitivity to planning and delivery errors. Before such implementation, it is necessary to understand how the sensitivity to different gamma criteria settings affects gamma passing rates (GPR). Methods: 752 IMRT QA measurements were re-analyzed with varying distance to agreement (DTA) and dose difference (DD) percentages using a Matlab program. Other quantifying information such as the mean dose difference in the treatment target (defined as points that are greater than 80% of maximal dose) were stored in a relational database for retrospective analysis. Results: The average and standard deviation of GPR (%) fell from 99.84 ± (0.43) to 89.61 ± (6.08) when restricting DD from 5 − 1% respectively, as compared to a drop from 99.15 ± (1.19) to 95.00 ± (4.43), when restricting the DTA from 5 − 1 mm respectively. The mean dose difference (%) in the treatment target between measured and calculated dose was −1.96 ± (0. 83), −0.09 ± (0.98), and 1.44 ± (0. 86) for each of our institution’s three matched linear accelerators (LINAC 1, 2, and 3 respectively). For plans that are approximately 2.7 sigma below the mean GPR, an average of 78.4% of those plans were measured on LINAC 1 or 3, while only 48% of the total plans were run on those machines. Conclusion: The data demonstrates that when restricting gamma criterion, such as the DD, the greatest indicator of reduced GPR in our institution is which matched LINAC the plan was measured on. While small, these differences manifest themselves to levels comparable to other treatment related differences and possibly confound the gamma analysis.

61 RADIATION PROTECTION AND DOSIMETRY↗