HOME | DOWNLOAD | DOCUMENTATION | FAQ |
> Home > Documentation
Version History
Version 3.4 (August 30, 2016)
- Add pyshtools to PyPI repository. Can now be installed using pip install pyshtools.
- Add new function SHBiasKMask which is the arbitrary window counterpart to the spherical cap window SHBiasK.
- Add get_biasedpowerspectrum() method to SHWindows for arbitrary windows.
- Add copy() method to all classes, which returns a deep copy of the instance.
- Add __sub__, __add__, __rsub__, __radd__, __mul__, __div__, __truediv__, and __pow__ operators for two sets of coefficient or grid classes, or one coefficient or grid class and a scalar.
- Add nwinrot option when rotating spherical cap windows in SHWindow that will rotate only the first nwinrot windows.
- Add degrees option to get_lats() and get_lons() methods.
- Add the constructor from_file() to initialize an SHGrid with a numpy formatted data file. Add option to read coeffs and grids from a numpy formatted binary file. Add tofile() methods to output raw grid and coefficient data as either text or binary formatted files.
- Update Intro 1 notebook and add example to Intro notebook 2 showing how to use arbitrary localization windows.
- Convert notebooks to html and add links to web documentation.
- Add option fixed_power to SHCoeffs.from_random() method to generate random coefficients that fit exactly the expected power spectrum.
- Add Earth topography coefficients referenced to mean sea level to the example files, expanded to degree 300: srtmp300.msl.
Version 3.3 (August 11, 2016)
- Added full support for Python 3.
- Fixed a critical, but rare, bug in MakeGridDH, MakeGravGridDH, MakeMagGridDH, and MakeGravGradGridDH. In these routines, the rows of the output grid are calculated by inverse Fourier transmforming a vector that depends upon the spherical harmonic coefficients. This vector, when using the complex-to-real FFTW routines, includes one element that corresponds to the coefficients with m=lmax+1. This element of the array was not properly initialized to zero in the Fortran code, and if the Fortran compiler did not explicitly zero all new arrays, this could have resulted in incorrect output. Given that this term is 1 index above the Nyquist frequency (m=lmax), if this element were not zero, each column of the grid would contain a component that oscillates from -1 to +1 (sclaed by the magnitude of the element). Even when this element was not initialized, a subsequent spherical harmonic expansion of the grid would usually give correct results.
- Added a setup.py file for easy installation.
- Makefiles were extensively modified to simplify the fortran and python builds: make all2 and make all3 were removed and replaced by flags that make use of the precompiler to resolve underscore problems; make install now places compiled module files in /usr/local/includes.
- Makefiles were improved to minimize problems when installing both fortran and fortran-mp components. When making the latter, all object and module files in the directory src are first deleted.
- The namespace of pyshtools was reorganized to list the routines by submodule name. A list of all routines is given in the submodule shtools.
- Add full support for the python classes SHCoeffs, SHGrid, and SHWindow.
- Added an info() method to shtools constants that prints an info string.
- Change python wrapper so that the output arrays in the SHExpand routines correspond to the optional input variable lmax_calc.
- Changed the primary input parameter of SHMTCouplingMatrix to be a matrix of power spectra of the localization windows instead of a matrix of spherical-harmonic coefficients of spherical-cap localization windows. Furthermore, the output dimensions of the matrix have been switched.
- Added two fortran routines for performing multitaper spectral analyses when using windows generaged from a mask, SHMultiTaperMaskSE and SHMultiTaperCSE.
- Minor modifications to the python example scripts.
- Minor bug fix to the scripts that create the unix man documentation.
- Add Python notebook tutorials.
- Created an SHTOOLS development fund, funded by bitcoin donations.
Version 3.2 (June 20, 2016)
- Added the optional argument centralmeridian to Curve2Mask that accounts for cases where the curve makes a complete circle in longitude about the planet.
- Fixed in bug in the python implementation where the outputs error and corr of SHAdmitCorr were switched.
- Added OpenMP support. When compiling with make fortran-mp, make fortran2-mp, or make fortran3-mp, saved variables in the subroutines are defined as being threadprivate.
- Optimized performance of the routines SHMultiTaperSE, SHMultiTaperCSE, and SHLocalizedAdmitCorr.
- Minor documentation fixes.
Version 3.1 (July 21, 2015)
- Added OSX installation support via brew.
- Added make install that copies files to /usr/local
- Reformatted all Fortran documentation files and rewrote the Python documentation and man pages.
- Added the routines CilmMinus and CilmMinusRhoH, which are the counterparts to CilmPlus and CilmPlusRhoH.
- Removed the following routines from the fortran documentation and Python wrappers: DhAj, NGLQ, NGLQSH, NGLQSHN.
- Removed the following redundant routines from SHTOOLS: ComputeD0, SHMTVarOpt0, SHSjkPG0.
- Renamed the routine PreCompute to SHGLQ.
- Renamed the routine YilmIndex to YilmIndexVector.
- Renamed the routine Hilm to BAtoHilm, and HilmRhoH to BAtoHilmRhoH.
- Renamed the routine Wl to DownContFilterMA, and WlCurv to DownContFilterMC.
- Added minimal support for three python classes: SHGrid, SHCoeffs, and SHWindow.
- Added the routine SHMTCouplingMatrix.
Version 3.0 (March 10, 2015)
- Added full python compatibility. This includes python wrappers for all SHTOOLS functions, python documentation for all functions, and a simple test suite to ensure that the SHTOOLS functions are working correctly.
- Moved the project from Sourceforge to GitHub. The Sourceforge project will no longer be maintained.
- Fixed bugs in SHrtoc and SHctor when the optional parameter CONVENTION was set to its default values, causing zeros to be returned.
- Modified ComputeDMap and SJHReturnTapersMap such that the parameter SAMPLING is optional.
- Fixed a bug in Curve2Mask that could give rise to vertical lines when the longitudes of the profile points were decreasing.
- Removed the routines Import_Wisdom_From_File and Export_Wisdom_From_File.
Version 2.9 (August 11, 2014)
- Modified CilmPlus and CilmRhoH so that they are accurate for powers beyond 13.
- Fixed a bug in MakeGeoidGrid that could give rise to runtime errors when the maximum dimension of CILM was less than LMAX.
- Transferred the discussion forum and download page to the Sourceforge project page.
- Minor fixes to documentation and web pages.
Version 2.8 (November 20, 2012)
- Following the suggestions of Leonid Petrov, the routines ShExpandDH, ShExpandDHC, ShExpandGLQ, ShExpandGLQC, MakeGridDH, MakeGridDHC, MakeGridGLQ, and MakeGridGLQC were optimized, giving rise to performance gains of about two.
- Following the suggestions of Leonid Petrov, the routine MakeGridPoint was optimized, giving rise to large performance gains.
- Replaced the routine MakeGravGrid2D with MakeGravGridDH. This routine is orders of magnitude faster than its predecessor, and can calculate all three components of the gravity field and gravity disturbance of the 2190-degree EGM2008 field in under 300 seconds.
- Added the routine MakeGravGradGridDH which calculates the components of the gravity "gradient" tensor on a flattened ellipsoid. This routine can calculate the degree-2190 EGM2008 gravity tensor in under 720 seconds.
- Replaced the routine MakeMagGrid2D with MakeMagGridDH, which is orders of magnitude faster than its predecessor.
- Added the routine NormalGravity to calculate the normal gravity on an ellipsoid.
- Modified SHLocalizedAdmitCorr to increase performance.
- Modified MakeCircleCoord such that the output coordinates are organized in a clockwise manner with the first index corresponding to due north.
- Introduced a different scaling procedure in CilmPlus and CilmRho that should be more accurate at high powers.
Version 2.7 (July 4, 2012)
- Improved memory management of CilmPlus and Hilm, avoiding potential runtime crashes at high degrees.
- Added routines CilmPlusRhoH and HilmRhoH (modified from CilmPlus and Hilm), which take into account lateral variations in density.
- Improved makefiles, including options to make and run example programs.
- Updated many of the example programs.
Version 2.6 (November 21, 2011)
- Fixed a bug in ShExpandDH, ShExpandDHC, ShExpandGLQ, ShExpandGLQC, MakeGridDH, MakeGridDHC, MakeGridGLQ, MakeGridGLQC that could have given inaccurate results if called more than twice in a program with the same N but different values of NORM.
- Fixed a bug in ShExpandDHC and SHExpandGLQC that could have given rise to a runtime crash or errors in the degree-0 coefficients when called with LMAX_CALC = 0.
- Removed routine SHReadGRACE and renamed SHReadCHAMP to SHRead2. SHRead2 now reads in time derivatives of terms, their epochs, and the starting and ending date of the gravity solution.
- Added the routine MakeEllipseCoord for calculating the (lat, long) coordinates of an ellipse placed on a sphere.
- Fixed a major bug in SHBiasK that would give wrong answers when SAVE_CG was set equal to 1.
- Fixed a minor bug in SHExpandGLQ and SHExpandGLQC that could have caused a runtime crash.
- Modified Hilm to accept Driscoll and Healy sampled grids, and added the optional parameter LMAX_CALC for calculating the spherical harmonic transforms up to a smaller value than the bandwidth corresponding to the input grids.
- Fixed a bug in CilmPlus that could cause an execution stop for Driscoll and Healy sampled grids when N was not equal to 2*LMAX+2.
- Modified the example code MarsCrustalThickness to use Driscoll and Healy sampled grids and to better remove high frequency aliases.
- Improved portability of routines RandomN and RandomGaussian.
Version 2.5 (August 20, 2009)
- Added the routines ComputeDMap and SHReturnTapersMap to compute the space-concentration kernel and its eigenvalues and eigenfunctions for an arbitrary concentration region. The routine Curve2Mask can be used to compute a mask of the concentration region given a set of coordinates that specifies a single closed curve.
- Added the optional argument NORMAL_GRAVITY to MakeGravGrid2d in order to remove the normal gravity from the total gravitational acceleration on the ellipsoid.
- Added routines SHCilmToVector and SHVectorToCilm for converting between 3D arrays of spherical harmonic coefficients, Cilm, and a 1D vector of ordered coefficients. The index corresponding to a given I, L and M can be calculated by YilmIndex.
- Added optional parameter K to EigValVecSym in order to only calculate and return the K largest eigenvalues and corresponding eigenvectors.
- Fixed a bug in ShExpandDH and SHExpandGLQ that could have given rise to a crash or errors in the degree-0 coefficients when called with LMAX_CALC = 0.
- Added the routine SHBiasAdmitCorr that calculates the expected localized admittance and correlation functions given the known cross-power spectra of two functions and the power spectra of the localizing windows.
- Added the optional parameter SAVE_CG to the routines SHBias and SHBiasK that allows the Clebsch-Gordon coefficients to be precomputed and saved for future calls.
- Added the optional parameter MTDEF to the routine SHLocalizedAdmitCorr, allowing the multitaper admittance and correlation to be computed in one of two ways. In the first, the multitaper cross-power spectra are first calculated, and then from these the admittance and correlation. In the second, the admittance and correlation are calculated for each individual taper, and then these are averaged.
- Added the option to use Driscoll-Healy or Gauss-Legendre quadrature grids with the routine CilmPlus. It is now necessary to specify the GRIDTYPE, and the argument W is now optional.
- Added the option DREF to CilmPlus in order to set the reference radius to calculate the relief and spherical harmonic coefficients. If this parameter is absent, the mean radius of the input grid will be used.
- Reordered the GRIDTYPES in MakeGeoidGrid in order to be consistent with CilmPlus.
- Fixed two bugs in SHrtoc and SHctor that only arose when the default value of 0 was explicitly specified for the optional parameter switchcs.
- Modified Wigner3j, based on the suggestions of Roelof Rietbroek, so that certain unused variables did not get divided by zero.
- Fixed a memory problem in the example code MarsCrustalThickness that might have prevented the file from compiling, or to have crashed at runtime.
- Fixed a memory allocation problem that could have caused a runtime crash in SphericalCapCoef.
- Removed routines SHDeltaX and SHDeltaL.
Version 2.4 (May 4, 2008)
- Added routines SHExpandDHC, MakeGridDHC, SHExpandGLQC, and MakeGridGLQC for performing complex spherical harmonic transforms and reconstructions.
- Added the optional parameter CNORM to the routines PlmBar, PlmBar_d1, PlmON, PlmON_d1, PlmSchmidt, PlmSchmidt_d1 and PreCompute in order to use the complex normalization.
- Added complex power spectrum routines SHPowerLC, SHPowerDensityLC, SHCrossPowerLC, SHCrossPowerDensityLC, SHPowerSpectrumC, SHPowerSpectrumDensityC, SHCrossPowerSpectrumC, and SHCrossPowerSpectrumDensityC.
- Fixed bug in MakeGrid2d, MakeMagGrid2d, and MakeGravGrid2d where the optional input of a west boundary was not properly converted to radians, and where 360 degrees was mistakenly added to EAST when EAST was greater than WEST (instead of the converse).
- Added the option to calculate all three components and magnitude of the gravitational field in MakeGravGrid2D. The output is now in SI units (as opposed to mGals before), and the radial gravitational acceleration is now positive when directed upwards (as opposed to downwards before).
- Added routines SHMagPowerL and SHMagPowerSpectrum for computing the power spectrum of the magnetic field given the spherical harmonic coefficients of the magnetic potential.
- Changed routine names of SHPowerSpectra to SHPowerSpectrum, SHPowerSpectraDensity to SHPowerSpectrumDensity, SHCrossPowerSpectra to SHCrossPowerSpectrum, and SHCrossPowerSpectrumDensity to SHCrossPowerSpectrumDensity.
- Fixed a bug where the spherical harmonic transforms and reconstructions would crash for LMAX=0.
Version 2.3 (October 4, 2007)
- Switched the naming convention of source files with a "2" appended at their end. If you previously compiled using Make all, you must now use Make all2, or vice versa.
- MakeGrid2D now requires the optional parameter semi-major axis A whenever the optional flattening parameter F is specified. Furthermore, the reference ellipsoid is now calculated exactly, instead of just using the degree-2 Legendre expansion coefficient as was done before.
- MakeGravGrid2D and MakeMagGrid2D now calculate fields on the ellipsoid using exact expressions for the ellipsoid, instead of just the degree-2 Legendre expansion coefficient as was done before. The semi-major axis must now be input, instead of the average planetary radius.
- Added optional parameter LMAX_CALC to MakeGridGLQ and MakeGridDH that specifies the maximum spherical harmonic degree to evaluate the function up to. The required parameter LMAX determines the grid spacing of the output grid, as well as the maximum expansion degree if LMAX_CALC is not present.
- Added the optional parameter LMAX_CALC to SHExpandGLQ and SHExpandDH that specifies the maximum degree of the spherical harmonic expansion. LMAX determines the grid spacing of the input grids, and corresponds to the maximum degree of the expansion if LMAX_CALC is not specified.
- MakeGeoid2D has been rewritten and renamed to MakeGeoidGrid. This new routine allows one to calculate the geoid using either first, second, or third order approximations, and the map values can be specified on either a 2D Cartesian grid, a Gauss-Legendre quadrature grid, a N by N grid, or a N by 2N grid. The geoid can either be referenced to a flattened ellipsoid, or a spherical interface.
- Fixed an error in the documentation to SHRotateCoef and SHRotateRealCoef. The rotations use the "y convention" where the second rotation is with respect to the y axis, and not the x axis as was originally stated.
- Added the optional parameters NORTH, SOUTH, EAST and WEST to the routines MakeGrid2D, MakeMagGrid2D and MakeGravGrid2D that specifies the maximum and minimum latitudes and longitudes of the output grid.
- Large temporary arrays are now dynamically allocated in order to resolve problems with Absoft ProFortran 10 and stack overflows.
Version 2.2 (May 20, 2007)
- Added routine MakeGravGrid that will create a map of the radial gravity or potential referenced to a flattened ellipsoid.
- Added routine SHConfidence that determines the probability two sets of spherical harmonic coefficients are correlated at a given degree with a given correlation coefficient.
- Added optional parameter f to the routine MakeGrid2D in order to reference the function to a flattened ellipsoid.
- Added optional parameter K1LINSIG to SHLocalizedAdmitCorr that allows one to calculate the expected uncertainty of the admittance function using a single localization window. The gravity and topography coefficients are assumed to be related by a linear (degree-dependent) transfer function and that the lack of correlation is a result of uncorrelated noise.
- Added optional argument CNUM to MakeCircleCoord that gives the number of elements in the output arrays.
- Fixed minor bugs in PlmSchmidt, PlmSchmidt_d1, PlmBar, PlmBar_d1 and PlmON_d1 that could have caused runtime crashes under certain circumstances.
- Fixed a minor bug in MakeGrid2d that could have caused a runtime crash.
- Fixed a minor bug in PlmON_d1 that could cause compilation errors.
- Fixed a minor bug in SHSjkPG that would cause g95 compiled code to crash at runtime.
- Saved allocated memory in the associated Legendre function routines (PlmSchmidt, PlmSchmidt_d1, PlmON, PlmON_d1, PlmBar, and PlmBar_d1) can be deallocated by calling with a spherical harmonic degree of -1.
- All routines that call associated Legendre functions now deallocate the saved allocated memory.
Version 2.1 (October 20, 2006)
- Modified PlmBar, PlmBar_d1, PlmON, PlmON_d1, PlmSchmidt and PlmSchmidt_d1 to save certain variables that are repeatedly calculated with multiple calls using the same maximum degree.
- MakeGridDH, MakeGridGLQ, SHExpandDH, and SHExpandGLQ have been rewritten so that the Legendre functions are computed within these routines, as opposed to being computed externally. This results in a reduction in computation time by a factor of about three to four.
- Added optional argument sampling that allows for the calculation of either equally sampled (NxN) and equally spaced (Nx2N) grids.
- Fixed a bug in MakeGridDH that would cause this routine to crash if the optional argument norm was not specified.
- Fixed SHMTVarOpt so that it would calculate all elements of the covariance matrix, not just those where the two tapers had the same value of m.
- Added optional argument nocross to SHMTVarOpt and SHMTVarOpt0 in order to calculate only the diagonal terms of the covariance matrix Fij.
- Fixed a bug in SHReturnTapersM where the eigenvalues and eigenvectors were wrong when m was not equal to zero.
- Added subroutine SHMTDebais that inverts for the global power spectrum given a localized multitaper spectrum estimate.
- Minor tweaks and improved documentation.
Version 2.0.1 (August 7, 2006)
- Updated the FAQ to describe common linking and compilation problems.
- make all and make all2 will compile the source code with and without underscores appended to the LAPACK subroutine names, respectively.
- Created the user group SHTOOLS at Google groups. This is probably the best place to look for answers to common problems.
- The optional argument csphase was mistakenly declared for the function PLegendre_d1 instead of PLegendreA_d1 in the SHTOOLS.f95 module file. This might have caused SHReturnTapersM to crash for m=0.
- Fixed a minor bug in the example file TestExpandDH.f95 that might have caused this program to crash.
- Fixed a minor array allocation problem in SHBiasK and SHBias that might have caused this routine to crash.
- Added example programs for determining the timing and accuracy of the spherical harmonic transform routines.
- Modified the FFTW wisdom example to only create a wisdom file for a single spherical harmonic degree.
- Minor tweaks.
Version 2.0 (July 21, 2006)
- Complete man and web based documentation has been written. In addition, all routines now check that the input arrays are dimensioned properly, and output informative errors if this is not the case.
- Subroutines Cilm_plus and DH_aj have been rename to CilmPlus and DHaj, respectively.
- MakeGridCoordGLQ renamed and completely written as GLQGridCoord.
- Subroutine ShannonNumber has been removed.
- Arguments of MakeGridDH have been reordered.
- Renamed SHReturnTapers to SHReturnTapers0.
- Renamed SHReturnAllTapers to SHReturnTapers.
- Added optional argument NORM to SHMultiply.
Version 1.7 (June 4, 2006)
- Modified SHExpandDH, MakeGridDH, SHExpandGLQ and MakeGridGLQ to make use of the symmetry relationship of Legendre functions about the equator. This leads to an increase in speed by almost a factor of 2 (note that the documentation under the section "Speed" was generated using the old routines).
- Renamed SHAdmitCoher to SHAdmitCorr.
- Modified SHDeltaX to work with non-zonal windows.
- Fixed bug in SHMultiply that caused this routine to crash.
- Added routine MakeCircleCoord for projecting circles on a map.
- Removed subroutines SHSjkPP, SHMTVarPP, SHBiasMT, SHBiasMTM, LmaxSimons, dpnlgndr and dpnlgndr_d1.
- Added functions SHSjkPG0 and SHjkPG.
- Added subroutines SHMTVarOpt0 and SHMTVarOpt for calculating optimal weights and minimum variances of multitaper spectral estimates using zonal and non-zonal tapers, respectively.
- Added subroutine SHReturnAllTapers to return all tapers of the spherical cap concentration problem (not just the zonal ones as in SHReturnTapers).
- SHExpandDH has been modified (as well as the example file TestExpandDH) to allow for grids that are N by 2N. This ensures that higher frequencies in longitude are not aliased into lower frequencies when using data expressed on equally spaced grids.
- SHBias was modified such that the routine now expects as input the window's power spectrum, and NOT the window's coefficients! Furthermore, the routine now assumes (as stated in the program) that the variable incspectra is zero beyond ldata.
- SHBiasK, SHMultiTaperCSE, SHMultiTaperSE, and SHLocalizedAdmitCorr have been rewritten to work with non-zonal tapers.
- SHrtoc and SHctor now allow to switch between Condon-Shortley phase conventions.
- SHRotateRealCoef now checks to see what the default value of the Condon-Shortley phase is before converting the real coefficients to complex form.
- PLegendreA and PLegendreA_d1 have been modified such that the default is to exclude the Condon-Shortley phase of (-1)m.
- All assocaiated Legendre function routines can use the Condon-Shortely phase factor by specifying the optional parameter csphase=-1.
- The default Condon-Shortley phase convention can be set at compile time by changeing the variable CSPHASE_DEFAULT in SHTOOLS.f95.
- Added the optional argument csphase to the routines SHExpandLSQ, MakeGridDH, SHExpandDH, MakeGrid2D, MakeGridPoint, MakeGridGLQ, SHExpandGLQ, PreCompute, SHMutiTaperCSE, SHMultiTaperSE, and SHMultiply. (Nevertheless, MakeMagGrid2D, MakeGeoid2D, Hilm, and Cilm_Plus never use the Condon-Shortely phase, regardless of what the default value is.)
Version 1.6 (November 12, 2005)
- SHTOOLS is now copyrighted according to the terms of the revised BSD license.
- Added "ortho-normalized" Legendre functions PlmON, PlmON_d1, PlON, and PlON_d1.
- Added subroutine SHReadCHAMP and SHReadGRACE that reads CHAMP and GRACE formatted spherical harmonic files, respectively.
- Fixed bug in routines SHrtoc and SHctor when the optional parameter convention was set equal to 2. This case was never used in previous versions of SHTOOLS.
- More constants have been added to the module Planets_Const.
- Added subroutine MakeGeoid2D that will create a 2D equally spaced grid of the geoid, calculated to second order, that includes the pseudo-rotational potential.
- Fixed a bug in SHExpandGLQ and SHExpandDH that only concerned the degree-0 term when using Schmidt-normalized harmonics.
- Fixed a bug in SHExpandDH that gave meaningless results for very last degree of the spherical harmonic coefficients.
Version 1.5 (October 7, 2005)
- PlmBar, PlmBar_d1, PlmSchmidt, and PlmSchmidt_d1 have been rewritten in order to give more accurate results for large m near the poles. The square roots of the integers are also precomputed leading to an increase in speed of about a factor of two.
- Added documentation concerning the accuracy and speed of the routines.
- Added make files to create a static archive of the compiled routines.
- Added optional argument chi2 to SHExpandLSQ, which is the residual sum of squares misfit for an overdetermined inversion.
- PreCompute, MakeGridGLQ, SHExpandGLQ, MakeGrid2D, MakeGridPoint, SHExpandDH and SHExpandLSQ now allow the option of using orthonormalized, unnormalized, or Schmidt normalized spherical harmonics.
- Added subroutine SHReadJPL that reads JPL formatted spherical harmonic files.
- Modified SHExpandDH so that it takes up much less memory when working with large spherical harmonic degrees.
- Modified PreCompute, SHExpandGLQ, MakeGridGLQ, Cilm_Plus, Hilm, and SHMultiply, as well as the example files MarsCrustalThickness and TestExpandGLQ, such that the array plx is optional. When working with large spherical harmonic degrees, this array takes up a lot of memory and should not be precomputed. If this array is not specified when calling the subroutine, then the optional array zero (obtained from PreCompute) must be specified instead.
Version 1.4 (August 7, 2005)
- Added optional argument taper_number to the function SHFindLWin.
- Added subroutines SHSjkPP, SHMTVarPP, SHMultiTaperSE, and SHMultiTaperCSE.
- SHLocalizedAdmitCoher now requires the mandatory calling parameter K which is the number of tapers to be used in the spectral estimations.
- Changed name of SHLocalizedAdmitCoher to SHLocalizedAdmitCorr.
- Added planets_const.
- Added optional parameter to SHrtoc and SHctor in order to choose between Varshalovich et al. and 4-pi geodesy complex coefficients.
- MakeGrid2D and MakeMagGrid2d have been speeded up by a factor of about 10.
Version 1.3 (August 17, 2004)
- Added routines SHExpandLSQ, RandomN, RandomGaussian, SHMultiply, ReadWisdom, ComputeD0, ComputeDm, ComputeDG82, SphericalCapCoef, LmaxSimons, SHDeltaX, SHDeltaL, SHReturnTapers, EigValVecSym, EigValVecSymTri, EigValSym, SHFindLWin, SHLocalizedAdmitCoher, SHAdmitCoher, Wigner3j, SHBias, SHBiasK, SHBiasMT, SHBiasMTM.
- Added optional argument skip to SHRead.
- PreCompute now tries to read a default system-wide FFTW wisdom file at /etc/fftw/wisdom, as well as an optional user defined FFTW wisdom file.
- Removed suffixes _FFTW from filenames and routines that use the FFTW routines.
Version 1.2
- Added routines MakeGridPoint, PlmSchmidt_d1, PlBar, PlBar_d1, PlSchmidt, PlSchmidt_d1, PlmBar_d1, PLegendreL, PLegendreL_d1, PLegendreA, and PLegendreA_d1.
- Added function declarations for wl and wl_curv to the interface block in SHTOOLS.
- Renamed the function index to PlmIndex.
- Renamed module FFTW to FFTW3.
- Added the optional argument lmax to the spherical harmonic conversion routines SHrtoc, SHctor, SHcilmtocindex, SHcindextocilm.
- Added the optional argument norm (which defines the spherical harmonic normalization) to the routines MakeGrid2D and MakeGridPoint.
Version 1.1 (May 26, 2004)
- Added routines SHExpandDH_FFTW and DH_aj.
> Home > Documentation
Institut de Physique du Globe de Paris | University of Sorbonne Paris Cité | © 2016 SHTOOLS |