Commit cc5a4523 by Javier

### doc html updated

parent dfdd630d
This diff is collapsed.
This diff is collapsed.
 Observables · juobs Documentation

Observables

juobs.meffFunction
meff(corr::Vector{uwreal}, plat::Vector{Int64}; pl::Bool=true, data::Bool=false )

meff(corr::Corr, plat::Vector{Int64}; pl::Bool=true, data::Bool=false)

Computes effective mass for a given correlator corr at a given plateau plat. Correlator can be passed as an Corr struct or Vector{uwreal}.

The flags pl and data allow to show the plots and return data as an extra result.

corr_pp = corr_obs.(data)
m = meff(corr_pp[1], [50, 60], pl=false)
source
juobs.dec_const_pcvcFunction
dec_const_pcvc(corr::Vector{uwreal}, plat::Vector{Int64}, m::uwreal, mu::Vector{Float64}, y0::Int64 ; pl::Bool=true, data::Bool=false)meff(corr::Corr, plat::Vector{Int64}; pl::Bool=true, data::Bool=false)

dec_const_pcvc(corr::Corr, plat::Vector{Int64}, m::uwreal; pl::Bool=true, data::Bool=false)

Computes decay constant using the PCVC relation for twisted mass fermions. The decay constant is computed in the plateau plat. Correlator can be passed as an Corr struct or Vector{uwreal}. If it is passed as a uwreal vector, vector of twisted masses mu and source position y0 must be specified.

The flags pl and data allow to show the plots and return data as an extra result.

corr_pp = corr_obs.(data)
m = meff(corr_pp[1], [50, 60], pl=false)
f = dec_const_pcvc(corr_pp[1], [50, 60], m, pl=false)
source
juobs.comp_t0Function
comp_t0(Y::YData, plat::Vector{Int64}; L::Int64, pl::Bool=false, rw::Union{Matrix{Float64}, Nothing}=nothing, npol::Int64=2)

comp_t0(Y::Vector{YData}, plat::Vector{Int64}; L::Int64, pl::Bool=false, rw::Union{Vector{Matrix{Float64}}, Nothing}=nothing, npol::Int64=2)

Computes t0 using the energy density of the action Ysl(Yang-Mills action). t0 is computed in the plateau plat. A polynomial interpolation in t is performed to find t0, where npol is the degree of the polynomial (linear fit by default)

The flag pl allows to show the plot.

#Single replica

t0 = comp_t0(Y, [38, 58], L=32)
t0_r = comp_t0(Y, [38, 58], L=32, rw=rw)

#Two replicas

t0 = comp_t0([Y1, Y2], [38, 58], L=32, pl=true)
t0_r = comp_t0(Y, [38, 58], L=32, rw=[rw1, rw2], pl=true)
source

read_mesons(path::String, g1::Union{String, Nothing}=nothing, g2::Union{String, Nothing}=nothing; id::Union{Int64, Nothing}=nothing)

read_mesons(path::Vector{String}, g1::Union{String, Nothing}=nothing, g2::Union{String, Nothing}=nothing; id::Union{Int64, Nothing}=nothing)

This function read a mesons dat file at a given path and returns a vector of CData structures for different masses and Dirac structures. Dirac structures g1 and/or g2 can be passed as string arguments in order to filter correaltors. ADerrors id can be specified as argument. If is not specified, the id is fixed according to the ensemble name (example: "H400"-> id = 400)

Examples:

source

Reads openQCD ms dat files at a given path. This method return YData:

• t(t): flow time values

• Ysl(icfg, x0, t): the time-slice sums of the densities of the Yang-Mills action

• vtr: vector that contains trajectory number

• id: ensmble id

Examples:

source

Reads openQCD ms1 dat files at a given path. This method returns a matrix W[irw, icfg] that contains the reweighting factors, where irw is the rwf index and icfg the configuration number. The function is compatible with the output files of openQCD v=1.2, 1.4 and 1.6. Version can be specified as argument.

Examples:

source

Reads openQCD pbp.dat files at a given path. This method returns a matrix md[irw, icfg] that contains the derivatives $dS/dm$, where $md[irw=1] = dS/dm_l$ and $md[irw=2] = dS/dm_s$

$Seff = -tr(log(D+m))$

$dSeff/ dm = -tr((D+m)^-1)$

Examples:

source
juobs.truncate_data!Function
truncate_data!(data::YData, nc::Int64)

truncate_data!(data::Vector{YData}, nc::Vector{Int64})

truncate_data!(data::Vector{CData}, nc::Int64)

truncate_data!(data::Vector{Vector{CData}}, nc::Vector{Int64})

Truncates the output of read_mesons and read_ms taking the first nc configurations.

Examples:

#Single replica
truncate_data!(dat, nc)
truncate_data!(Y, nc)

#Two replicas
dat = read_mesons([path1, path2], "G5", "G5")
truncate_data!(dat, [nc1, nc2])
truncate_data!(Y, [nc1, nc2])
source
 Search · juobs Documentation

 Search · juobs Documentation

This diff is collapsed.
This diff is collapsed.
 using Documenter, juobs push!(LOAD_PATH, "../src") makedocs(sitename = "juobs Documentation") makedocs(sitename = "juobs Documentation", doctest=true, repo = "https://gitlab.ift.uam-csic.es/jugarrio/juobs", pages = [ "Home" => "index.md", "Reader" => "reader.md", "Tools" => "tools.md", "Observables" => "obs.md", "Linear Algebra" => "linalg.md"], format = Documenter.HTML(prettyurls=false))
 # DOCUMENTATION @contents  ## READER @docs read_mesons read_ms read_ms1 read_md truncate_data!  ## Contents ## TOOLS @docs corr_obs md_sea lin_fit fit_routine @contents Pages = ["reader.md", "tools.md", "obs.md", "linalg.md"] Depth = 3  ## LINEAR ALGEBRA @docs uweigvals uweigvecs uweigen get_matrix energies getall_eigvals getall_eigvecs  ## OBSERVABLES @docs meff dec_const_pcvc comp_t0 
 # Linear Algebra @docs uweigvals uweigvecs uweigen get_matrix energies getall_eigvals getall_eigvecs 
docs/src/obs.md 0 → 100644
 # Observables @docs meff dec_const_pcvc comp_t0 
 # Reader @docs read_mesons read_ms read_ms1 read_md truncate_data! 
docs/src/tools.md 0 → 100644
 # Tools @docs corr_obs md_sea lin_fit fit_routine 
 ... ... @@ -254,7 +254,7 @@ This function computes the eigenvalues and the eigenvectors of a matrix of uwre If a second matrix b is given as input, it returns the generalised eigenvalues and eigenvectors instead. It takes as input: - a::Matrix{uwreal} : a matrix of uwreal - a::Matrix{uwreal} : a matrix of uwreal - b::Matrix{uwreal} : a matrix of uwreal, optional ... ...
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!