Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Mauve diprotodon #15

Open
wants to merge 11 commits into
base: main
Choose a base branch
from
Open

Mauve diprotodon #15

wants to merge 11 commits into from

Conversation

GroundhogState
Copy link
Collaborator

Additions to core_BEC_analysis

Jacob 2020-09-22

coldatom_tools

doppler_broadening

Returns the multiplicative factor to obtain FWHM of a delta-function distribution
at frequency f0 for atoms of mass m at temperature T. Particle mass input optional, default to mHe. eg:
doppler_width = doppler_broadening(T,m)*f0

doppler_limit

Accepts a transition linewidth (Hz) and returns the Doppler cooling limit (in K)

deBroglie_wavelength

Accepts temperature and particle mass (optional, defaults to mHe). Removed redundant lamnda_dB().

txy_tools

show_txy_raw

A function for flexible visualization and investigation of TXY data. Functions similarly to dld_front_panel but can be fed TXY data structures from memory rather than importing. Useful for quick checks of BEC structure or for quantitative analysis of count flux in cartesian histograms.
Mandatory inputs:

  • data_in: either an N x 3 array or a struct with fields {shot_num,num_counts,counts_txy} - compatible with the standard output of import_tdc_data.
    Outputs:
  • 3D histogram counts, associated bin edges, and normalized flux density averaged over all shots (assumes num_shots=1 if passed an Nx3 array).
  • 1D and 2D histograms obtained by summing over
  • Visual presentation of 1D and 2D histograms. Can be suppressed by setting draw_plots argument to false.
    Optional inputs include log_plot, blur_size, num_bins and others.

getlims

Returns the bounding coordinate values of a set of N data points in D dimensions, returning an Nx2 array.

CLI_tools

sshow

Displays a tree diagram of a given structure s including the size and type of the leaves (assuming leaves terminate in an array of some kind). Removed redundancy (in /command_line_format_tools and /CLI_tools)

simple_utilities

intCDF

Returns the cumulative distribution function of data x over
a domain d. Y = intCDF(x,d) returns a row vector of the same size as x,
whose entries Y(i) counts the number of elements of d that are less than x(i).
Inputs:
x - data - 1D real vector
d - domain - 1D real vector
Example

domain = 0:5;
data = [1,3,3.1]
isequal(intCDF(domain,data),[0,0,1,1,3,3])

add_unc

Accepts tuples of the form [x,x_unc] and returns the pair [val,unc]. Would be good to write a more general function that can accept [val,unc] pairs arrays and propagate uncertainties automatically through more complicated functions.

Updates to existing functions

bec_properties

  • Accepts keyword argument inputs for temperature, mass, and a_scat_len (defaults to 0, mHe, and aHe_scat from hebec_constants)
  • Accepts arrays of BEC_number and omega has been input
  • Corrected ambiguity in the choice of number (condensed vs. total) in some formulae
  • Returns additional quantities:
    • First-order Tan two-body contact
    • Lee-Huang-Yang correction to energy
    • Speed of sound (at peak density)
    • Oscillator length
    • Healing length
    • Condensed fraction
    • De Broglie wavelength (1 if temperature not input)

coldatom_tools

g_bose

  • Accepts additional optional argument tol_err which specifies the convergence cutoff. Defaults to 1e-16 but considerable speedup can be had by using a less stringent tolerance.

Removed:

nucellf
A shorthand for adding the ('UniformOutput',false) kwarg pair to cellfun. Can use inbuilt function cellfun(@(x) f(x), x, 'uni', 0) instead.

Notes:

cli_header contains some of the functionality of cli_format_text.
cli_format_text is strictly more general (and recommended). Removing cli_header from core_BEC and replacing calls throughout all projects is an unnecessary expense as they are (intended to be) user-interface tools rather than functional components.

sltext is recommended for removal as has no known usages

…ysis

Merging to include OuterProduct, sshow, r2k and k2r
damnit vim
… colour generation, fake shot (thermal + QD) generator, some tools in dev'
… inputs. integrated into show_txy_raw; could potentially amend input structure of txy2vzxy to enable seamless replacement of txy_to_vel
…n edges to centres in quick command), could probably be tensorified. Also broke out dev/find_sat_peak including a test script, could be extended to higher dims
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant