Skip to content
Open
Show file tree
Hide file tree
Changes from 5 commits
Commits
Show all changes
39 commits
Select commit Hold shift + click to select a range
9324194
add function for computing GL FSE correction factors
kostrzewa Feb 20, 2020
0444b7b
remove outdated documentation and unnecessary parameters from gl_fse_…
kostrzewa Feb 21, 2020
d965f4c
Merge remote-tracking branch 'origin/master' into gl_fse
kostrzewa Jul 18, 2020
18b1b36
have explicit source file for generic functions
kostrzewa Jul 28, 2020
2540be1
use .Rdata as extension in analysis_online
kostrzewa Jul 28, 2020
9475263
begin work on 'tseries' class
kostrzewa Jul 28, 2020
3d95596
vignette in which the new 'tseries' class is applied to study the gra…
kostrzewa Jul 28, 2020
eeaa93a
add 'tseries' vignette
kostrzewa Jul 28, 2020
4c88db7
In the analysis of online measurements, indicate M_ps extracted from …
kostrzewa Apr 7, 2021
02afcc0
store WIP of timeseries toolset and new implementation of gradient fl…
kostrzewa Jul 21, 2021
5bdf9d9
two digits for error on scale in gradient flow evolution plots
kostrzewa Jul 21, 2021
c87a928
print statement which shouldn't exist in readcmifiles
kostrzewa Dec 15, 2021
aa4cdab
add to parameters of as a way to have different starting points bet…
kostrzewa Jan 13, 2022
85f11a8
speed up readcmifiles by a factor of more than 100
kostrzewa Mar 5, 2022
ce93f04
Merge remote-tracking branch 'origin/master' into tseries
kostrzewa Mar 7, 2022
5aa3c81
resolve documentation warnings due to missing names in R/hadron-packa…
kostrzewa Mar 7, 2022
8550080
speed up gradient flow I/O by about a factor of 200
kostrzewa Mar 7, 2022
c26cbc2
remove 'skip' from analysis_online, replace this functionality for 'o…
kostrzewa Apr 20, 2022
73b407c
be more helpful when omeas.offset is used and no trajectories are found
kostrzewa Apr 20, 2022
597efe4
note dependency on reshape2
kostrzewa Jul 23, 2022
7d7d38f
account for dimesions of gradient flow reference data
kostrzewa Jul 23, 2022
9e8d133
add `trajectory_length` parameter for `analysis_online` which can be …
kostrzewa Jul 23, 2022
accc42c
in analysis_tmlqcd_gradient_flow, plot and summarize the topological …
kostrzewa Jul 26, 2022
b8690b9
draw error rectangles for mpcac plateau only if the error is not too …
kostrzewa Aug 20, 2022
e642956
introduce some restrictions on the plotting of error polygons in the …
kostrzewa Nov 11, 2022
7288fbc
add more documentation to 'analysis_tmlqcd_gradient_flow'
kostrzewa Dec 16, 2022
ade423b
typo in evals.stepsize
kostrzewa Jan 19, 2023
34272fc
Allow the definition point for gluonic scales (scale_definition) to b…
kostrzewa Jan 20, 2023
ed17427
parenthesis problem in analysis_online
kostrzewa Jan 21, 2023
cafec48
do not use the bootstrap samples to estimate the mean of mpi_ov_fpi
kostrzewa Nov 18, 2023
306fa22
in boostrap.nlsfit, make sure that x divides y
kostrzewa Feb 4, 2024
4a4809e
generalise the generator functions for CVC / meson_2pt keys to suppor…
kostrzewa Feb 4, 2024
3018fd4
nyom HDF5 files differ from CVC HDF5 files in that they use a compoun…
kostrzewa Feb 4, 2024
065c2b2
do not use 'sample', use 'istoch' instead
kostrzewa Feb 5, 2024
94922b3
remove unused dbg parameter
kostrzewa Feb 9, 2024
8824f0f
ensure that timeseries dat contains 'y' and 't'
kostrzewa Feb 9, 2024
829f13a
Merge remote-tracking branch 'origin/gl_fse' into tseries
kostrzewa Jun 19, 2024
5876354
various small updates
kostrzewa Oct 2, 2025
cac77b9
fix some of the roxygen warnings
kostrzewa Oct 2, 2025
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 3 additions & 3 deletions R/analysis_online.R
Original file line number Diff line number Diff line change
Expand Up @@ -120,7 +120,7 @@ analysis_online <- function(L, Time, t1, t2, beta, kappa, mul,
}

# store analysis results in practical R format, replacing entries as new data is added
resultsfile <- "omeas.summary.RData"
resultsfile <- "omeas.summary.Rdata"

resultsum <- list()
if(file.exists(resultsfile)){
Expand Down Expand Up @@ -227,8 +227,8 @@ analysis_online <- function(L, Time, t1, t2, beta, kappa, mul,
if(trajlabel){
filelabel <- sprintf("%s_traj%06d-%06d",filelabel,min(omeas.cnums),max(omeas.cnums))
}
message("Writing online measurements RData to ", sprintf("onlineout.%s.RData",filelabel), "\n")
save(onlineout,file=sprintf("onlineout.%s.RData",filelabel))
message("Writing online measurements Rdata to ", sprintf("onlineout.%s.Rdata",filelabel), "\n")
save(onlineout,file=sprintf("onlineout.%s.Rdata",filelabel))

plotcounter <- plotcounter+1
dpaopp_filename <- sprintf("%02d_dpaopp_%s",plotcounter,filelabel)
Expand Down
14 changes: 0 additions & 14 deletions R/bootstrap.nlsfit.R
Original file line number Diff line number Diff line change
Expand Up @@ -1278,20 +1278,6 @@ plot.bootstrapfit <- function(x, ..., col.line="black", col.band="gray", opacity
lines(x=X, y=Y, col=col.line, lty=lty, lwd=lwd)
}

#' residual_plot
#'
#' generic residual_plot method
#'
#' @param x the object to plot
#' @param ... additional parameters to be passed on to specialised functions
#'
#' @return
#' No return value.
#'
#' @export
residual_plot <- function (x, ...) {
UseMethod("residual_plot", x)
}

#' @export
residual_plot.bootstrapfit <- function (x, ..., error_fn = x$error.function, operation = `/`) {
Expand Down
4 changes: 2 additions & 2 deletions R/cf.R
Original file line number Diff line number Diff line change
Expand Up @@ -419,7 +419,7 @@ is_empty.cf <- function (.cf) {
#' with elements `boot`, `boot.R`, `boot.l`, `sim`, `endcorr`,
#' `resampling_method`, `boot_dim`, `icf` and, optionally
#' `iboot_dim` (if both `cf1` and `cf2` contain imaginary parts).
resampling_is_compatible <- function(cf1, cf2){
resampling_is_compatible.cf <- function(cf1, cf2){

res <- list()
res$boot <- ( inherits(cf1, 'cf_boot') & inherits(cf2, 'cf_boot') )
Expand Down Expand Up @@ -453,7 +453,7 @@ resampling_is_compatible <- function(cf1, cf2){
#' with elements `boot`, `boot.R`, `boot.l`, `sim`, `endcorr`,
#' `resampling_method`, `boot_nrow`, `icf` and, optionally
#' `iboot_nrow` (if both `cf1` and `cf2` contain imaginary parts).
resampling_is_concatenable <- function(cf1, cf2){
resampling_is_concatenable.cf <- function(cf1, cf2){
res <- list()
res$boot <- ( inherits(cf1, 'cf_boot') & inherits(cf2, 'cf_boot') )
res$seed <- (cf1$seed == cf2$seed)
Expand Down
114 changes: 114 additions & 0 deletions R/generics.R
Original file line number Diff line number Diff line change
@@ -0,0 +1,114 @@
#' residual_plot
#'
#' generic residual_plot method
#'
#' @param x the object to plot
#' @param ... additional parameters to be passed on to specialised functions
#'
#' @return
#' No return value.
#'
#' @export
residual_plot <- function (x, ...) {
UseMethod("residual_plot", x)
}

#' generic function to extract a fitted mass
#'
#' @description
#' One of the main analysis tasks in \link{hadron} is the estimation
#' of energy levels or masses from correlation functions. The
#' corresponding analysis functions return objects, typically lists,
#' containing the masses or energy levels. `extract_mass` is a
#' generic function to extrac such fitted mass values.
#'
#' @param object Object to extract the mass from.
#'
#' @return Numeric. The mass value.
#'
#' @export
extract_mass <- function (object) {
UseMethod('extract_mass')
}

#' generic function to check if resampling samples are compatible
#'
#' @description
#' When binary operations are performed on resampled data, it is
#' necessary to check if the samples are compatible, at the very
#' least they must have the same dimensions.
#'
#' @details
#' Note that since R's double dispatch doesn't really work with S3
#' classes, the class of \code{x} decides which method is called.
#'
#' @param x lhs object
#' @param y rhs object
#'
#' @return
#' list of booleans which correspond to the results of checks
#' for equality of different properties of the resampling samples
#'
#' @export
resampling_is_compatible <- function(x,y){
UseMethod('resampling_is_compatible', x)
}

#' generic function to check if resampling samples are concatenable
#'
#' @description
#' When we want to combine resampled data along an axis orthogonal
#' to the axis of samples (for example, if we want to turn two
#' vectors with 'boot.R' samples into a matrix with 2 columns),
#' then we need to check if the number of samples for both
#' data are the same.
#'
#' @details
#' Note that since R's double dispatch doesn't really work with S3
#' classes, the class of \code{x} decides which method is called.
#'
#' @param x lhs object
#' @param y rhs object
#'
#' @return
#' list of booleans which correspond to the results of checks
#' for equality of different properties of the resampling samples
#'
#' @export
resampling_is_concatenable <- function(x,y){
UseMethod('resampling_is_concatenable', x)
}

#' generic function to multiply two objects with each other
#'
#' @details
#' Note that since R's double dispatch doesn't really work with S3
#' classes, the class of \code{x} decides which method is called.
#'
#' @description
#'
#'
#' @export
mul <- function(x, y, ...){
UseMethod('mul', x)
}

#' generic function to add two objects to each other
#'
#' This function provides
#'
#' @export
add <- function(x, y, ...){
UseMethod('add', x)
}

#' @export
subtract <- function(x, y, ...){
UseMethod('subtract', x)
}

#' generic function to divide objects by each other
#' @export
div <- function(x, y, ...){
UseMethod('div', x)
}
18 changes: 0 additions & 18 deletions R/removeTemporal.cf.R
Original file line number Diff line number Diff line change
Expand Up @@ -291,24 +291,6 @@ dispersion_relation <- function (energy, momentum_d, extent_space, plus = TRUE,
return (energy_out)
}

#' generic function to extract a fitted mass
#'
#' @description
#' One of the main analysis tasks in \link{hadron} is the estimation
#' of energy levels or masses from correlation functions. The
#' corresponding analysis functions return objects, typically lists,
#' containing the masses or energy levels. `extract_mass` is a
#' generic function to extrac such fitted mass values.
#'
#' @param object Object to extract the mass from.
#'
#' @return Numeric. The mass value.
#'
#' @export
extract_mass <- function (object) {
UseMethod('extract_mass')
}

#' specialisation of \link{extract_mass} to objects of type
#' `effectivemassfit`
#'
Expand Down
Loading