GPUE  v1.0
GPU Gross-Pitaevskii Equation numerical solver for Bose-Einstein condensates
velField.m File Reference

Go to the source code of this file.

Functions

end imagesc (y(valsx), x(valsy), sqrt(v_x(valsy, valsx).^ 2+v_y(valsy, valsx).^ 2))
 
 set (gca, 'TickLabelInterpreter', 'latex')
 
 set (gca, 'DefaultTextInterpreter', 'Latex') set(gca
 
Latin Modern xlabel (' $x$(m)')
 
 ylabel (' $y$(m)')
 
 set (h, 'TickLabelInterpreter', 'latex') ylabel(h
 
 V (m/s)') if length(lims~=0) axis(lims)
 

Variables

 function []
 
and assumes dx
 
and assumes too high and very sparse Start at
 
and assumes too high and very sparse Start and % increase until happy % x
 
and assumes too high and very sparse Start and % increase until happy % y are the grid spacings along the x and y axis % normed specifies whether to normalise the vector directions if yes
 
and assumes too high and very sparse Start and % increase until happy % y are the grid spacings along the x and y axis % normed specifies whether to normalise the vector directions if % otherwise % lims is [xMin xMax yMin yMax] Hides the edge garbage Otherwise
 
and assumes too high and very sparse Start and % increase until happy % y are the grid spacings along the x and y axis % normed specifies whether to normalise the vector directions if % otherwise % lims is [xMin xMax yMin yMax] Hides the edge garbage phase = angle(wfc0)
 
 np1 =unwrap(phase,[],1)
 
 np2 =unwrap(phase,[],2)
 
 hbar =1.05457e-34
 
 v_y = (hbar/m)*(velnp1y)
 
 v_x = (hbar/m)*(velnp2x)
 
 valsx =[1:incr:size(wfc0,1)]
 
 % [200:1:300]
 
 valsy =valsx
 
% if normed
 
else L =1
 
hold on q = quiver(y(valsx),x(valsy),v_x(valsy,valsx)./L,v_y(valsy,valsx)./L,'AutoScaleFactor',0.5,'Color','w')
 
axis square
 
hold off
 
 FontName
 
Latin Modern Roman
 
Latin Modern FontSize
 
 h =colorbar
 

Function Documentation

◆ imagesc()

end imagesc ( y(valsx ,
x(valsy ,
sqrt(v_x(valsy, valsx).^ 2+v_y(valsy, valsx).^ 2)   
)

◆ set() [1/3]

set ( gca  ,
'TickLabelInterpreter ,
'latex  
)

◆ set() [2/3]

set ( gca  ,
'DefaultTextInterpreter'  ,
'Latex'   
)

◆ set() [3/3]

set ( h  ,
'TickLabelInterpreter ,
'latex  
)

◆ V()

V ( m s)
pure virtual

◆ xlabel()

Latin Modern xlabel ( ' $x$(m)'  )

◆ ylabel()

ylabel ( ' $y$(m)'  )

Variable Documentation

◆ %

%[200:1:300]

Definition at line 26 of file velField.m.

◆ at

and assumes too high and very sparse Start at

Definition at line 5 of file velField.m.

◆ dx

and assumes dx
Initial value:
==dy
% m is the mass of the atomic species
% incr is the increment over which to calculate field direction. Too low
% and it may be very dense
tuple dy
Definition: en.py:60
% Used for determining interaction strength of system % % Constants required for simulation m
Definition: GPE_2d.m:13
% If at least one of the indices is
incr
Definition: en.py:57
tuple mass
Definition: observables.py:72

Definition at line 5 of file velField.m.

◆ FontName

FontName

Definition at line 43 of file velField.m.

◆ FontSize

Latin Modern FontSize

Definition at line 43 of file velField.m.

◆ function

function[]
Initial value:
= velField(wfc0,dx,m,incr,x,y,normed, lims)
%velField calculates velocity field of wavefunction
% velField plots the magnitude and direction of the velocity field of the
% given wavefunction wfc0 in 2D.
% dx is the increment along x
and assumes dx
Definition: velField.m:5
and print the results to file Essentially wraps % voronoi2dCellColour and saves the resulting plots
Definition: VtxCorr.m:1
% Used for determining interaction strength of system % % Constants required for simulation m
Definition: GPE_2d.m:13
% if normed
Definition: velField.m:30
and assumes too high and very sparse Start and % increase until happy % x
Definition: velField.m:5
% If at least one of the indices is
incr
Definition: en.py:57

Definition at line 1 of file velField.m.

◆ h

Definition at line 45 of file velField.m.

◆ hbar

hbar =1.05457e-34

Definition at line 22 of file velField.m.

◆ L

else L =1

Definition at line 33 of file velField.m.

◆ normed

% if normed
Initial value:
==1
valsx
Definition: velField.m:26
v_y
Definition: velField.m:23
end if sqrt(sum([x(ii), y(ii)].^ 2))< radius %% ignore edges if(length(DT.vertexAttachments
else L
Definition: velField.m:33
valsy
Definition: velField.m:27
v_x
Definition: velField.m:24

Definition at line 30 of file velField.m.

◆ np1

np1 =unwrap(phase,[],1)

Definition at line 16 of file velField.m.

◆ np2

np2 =unwrap(phase,[],2)

Definition at line 17 of file velField.m.

◆ off

hold off

Definition at line 39 of file velField.m.

◆ Otherwise

and assumes too high and very sparse Start and % increase until happy % y are the grid spacings along the x and y axis % normed specifies whether to normalise the vector directions if % otherwise % lims is [xMin xMax yMin yMax] Hides the edge garbage Otherwise

Definition at line 5 of file velField.m.

◆ phase

and assumes too high and very sparse Start and % increase until happy % y are the grid spacings along the x and y axis % normed specifies whether to normalise the vector directions if % otherwise % lims is [xMin xMax yMin yMax] Hides the edge garbage phase = angle(wfc0)

Definition at line 15 of file velField.m.

◆ q

hold on q = quiver(y(valsx),x(valsy),v_x(valsy,valsx)./L,v_y(valsy,valsx)./L,'AutoScaleFactor',0.5,'Color','w')

Definition at line 38 of file velField.m.

◆ Roman

Latin Modern Roman

Definition at line 43 of file velField.m.

◆ square

axis square

Definition at line 38 of file velField.m.

◆ v_x

v_x = (hbar/m)*(velnp2x)

Definition at line 24 of file velField.m.

◆ v_y

v_y = (hbar/m)*(velnp1y)

Definition at line 23 of file velField.m.

◆ valsx

valsx =[1:incr:size(wfc0,1)]

Definition at line 26 of file velField.m.

◆ valsy

valsy =valsx

Definition at line 27 of file velField.m.

◆ x

and assumes too high and very sparse Start and % increase until happy % x

Definition at line 5 of file velField.m.

◆ yes

and assumes too high and very sparse Start and % increase until happy % y are the grid spacings along the x and y axis % normed specifies whether to normalise the vector directions if yes

Definition at line 5 of file velField.m.