PK8H~z.fityk-v1.3.0/objects.inv# Sphinx inventory version 2 # Project: Fityk # Version: 1.3.0 # The remainder of this file is compressed using zlib. xڭMs6Ɂ5Ī37#" C?VbjHxwc\i l{}z<+߃rd[u=3UT7eUy(yj^#|CEF(H=?ɁaPFd;SԖL^1 fۓtOfx}3Юpuu%>ZEuuowa N973` UVBݢQi!lSc|(T5yz34!$ SզgN` fAD|; ڸq5[.F;lRd߀5%F"Gij1d 9'*My f;B.zZnzԔIWԜHfNlJa]tս*C}:zUZ1bœ-fP,zl'3s~o*mrVW6^6N?# 1: Sz~`wO㝢`mr8| W &G3"a6juO߫eʵdS3ћo88yry<Y|3hbnʭPo:($n(uՖ?=32CSRZXO^p%KĿլiGNEпZ&±&8qM]na Sm Cj?DQQ텔|*V!Wâ:8ыƌ#zt?IIZCa`ֵ _OnfC0tH0p$UX f#a7v65Lk90v\m]OMgq~V"{^߽/Lzr釥^MRG>+1k8z۸v]~tpQ晏Xi5/f#:#HYLR(R0osw-aa78^ z`˱%b⨨j\FHXnP66"/eXDK$ HCrtR =&Ǚ)̎_('2ďaɈ6:c u+tG_?%EL&Mh͹fl578PK8Hm$!fityk-v1.3.0/.buildinfo# Sphinx build info version 1 # This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done. config: tags: PK8Hw?E22fityk-v1.3.0/index.html Contents — Fityk 1.3.0 manual
PK8H|/w>v>vfityk-v1.3.0/fityk-manual.html Fityk 1.3.0 manual

Fityk 1.3.0 - User’s Manual

Introduction

What’s This Program for?

Fityk is a program for nonlinear fitting of analytical functions (especially peak-shaped) to data (usually experimental data). The most concise description: peak fitting software. There are also people using it to remove the baseline from data, or to display data only.

It is reportedly used in crystallography, chromatography, photoluminescence and photoelectron spectroscopy, infrared and Raman spectroscopy, to name but a few. Although the author has a general understanding only of experimental methods other than powder diffraction, he would like to make it useful to as many people as possible.

Fityk offers various nonlinear fitting methods, simple background subtraction and other manipulations to the dataset, easy placement of peaks and changing of peak parameters, support for analysis of series of datasets, automation of common tasks with scripts, and much more. The main advantage of the program is flexibility - parameters of peaks can be arbitrarily bound to each other, e.g. the width of a peak can be an independent variable, the same as the width of another peak, or can be given by complex (and general for all peaks) formula.

GUI vs CLI

The program comes in two versions: the GUI (Graphical User Interface) version - more comfortable for most users, and the CLI (Command Line Interface) version (named cfityk to differentiate).

technical note

The GUI version is written using the wxWidgets library. It runs on Unix species with GTK+, on MS Windows and (with some problems) on MacOS X.

If the CLI version was compiled with the GNU Readline Library, command line editing, TAB-expanding and command history will be available. On Unix, the gnuplot program can be used for data visualization.

About...

Fityk is free software; you can redistribute and modify it under the terms of the GPL, version 2 or (at your option) any later version.

To download the latest version of the program or to contact the author visit http://fityk.nieto.pl/.

Reference for academic papers: M. Wojdyr, J. Appl. Cryst. 43, 1126-1128 (2010) [reprint]

This manual is written in ReStructuredText. All corrections and improvements are welcome. Use the Show Source link to get the source of the page, edit it and send me either the modified version or a patch.

Alternatively, go to GitHub, open corresponding rst file, press Fork and edit this file button, do edits in your web browser and click Propose file change.

The following people have contributed to this manual (in chronological order): Marcin Wojdyr (maintainer), Stan Gierlotka, Jaap Folmer, Michael Richardson.

One of the fitting methods uses MPFIT library (MINPACK-1 Least Squares Fitting Library in C), which includes software developed by the University of Chicago, as Operator of Argonne National Laboratory.

Getting Started

Graphical Interface

That’s how the GUI looks like:

_images/fityk-with-tooltip.png

The input field with the output window provide a console-like interface to the program. All important operation performed using the GUI are translated to textual commands and displayed in the output windows.

The main plot can display data points, model that is to be fitted to the data and component functions of the model. Use the pop-up menu (click right button on the plot) to configure it. Some properties of the plot (e.g. colors of data points) can be changed from the sidebar.

The helper plot below can show various statistics. The default one one is the difference between the data and the model. The plot can be configured from a pop-up menu. If you do not find it useful you may hide it (GUI ‣ Show).

Configuration of the GUI (visible windows, colors, etc.) can be saved using GUI ‣ Save current config.

On the main plot, the meaning of the left and right mouse button depends on current mode. There are hints on the status bar. In normal mode, the left button is used for zooming and the right invokes the pop-up menu.

On the helper plot, selecting a horizontal range with the left button will show this range, automatically adjusting vertical range. The middle button shows the whole dataset (like Zoom All in the toolbar). The right button displays a menu.

Minimal Example

Let us analyze a diffraction pattern of NaCl. Our goal is to determine the position of the center of the highest peak. It is needed for calculating the pressure under which the sample was measured, but this later detail in the processing is irrelevent for the time being.

The data file used in this example is distributed with the program and can be found in the samples directory.

First load data from file nacl01.dat. You can do this by typing:

@0 < nacl01.dat

The CLI version of the program is all about typing commands. From time to time it is also handy to type a command in the GUI, but usually the GUI provides more intuitive, mouse-driven way to perform the same operation – it is mentioned in the the grey boxes below.

In the GUI

select Data ‣ Load File from the menu (or Load Data from the toolbar) and choose the file.

If you use the GUI, you can zoom-in to the biggest peak using left mouse button on the residual plot (the plot below the main plot). To zoom out, press the View whole toolbar button.

Now all data points are active. Because only the biggest peak is of interest for the sake of this example, the remaining points can be deactivated:

A = (x > 23.0 and x < 26.0)

In the GUI

change to the range mode (toolbar: Data-Range Mode) and deactivate not needed points with the right mouse button.

As our example data has no background to worry about, our next step is to define a peak with reasonable initial values and fit it to the data. We will use Gaussian. To see its formula, type: info Gaussian (or i Gaussian) or look for it in the section Built-In Functions.

To add a peak, either set the initial parameters manually:

F += Gaussian(~60000, ~24.6, ~0.2)

In the GUI

it is also possible to set the initial parameters with the mouse: change the GUI mode to Add-Peak Mode, click on the plot and drag the mouse to select the position, height and width of a new peak.

or let the program guess it:

guess Gaussian

In the GUI

select Gaussian from the list of functions on the toolbar and press Auto Add.

If the functions are not named explicitely (like in this example), they get automatic names %_1, %_2, etc.

Now let us fit the function. Type: fit.

In the GUI

select Fit ‣ Run from the menu or press Fit.

Important

Fitting minimizes the weighted sum of squared residuals (see Nonlinear Optimization). The default weights of points are not equal.

To see the peak parameters, type: info prop %_1.

In the GUI

move the cursor to the top of the peak and try out the context menu (the right mouse button), or check the parameters on the sidebar.

That’s it!

You can save all the issued commands to a file:

info history > myscript.fit

and later use it as a macro:

exec myscript.fit

In the GUI

use Session ‣ Save History and Session ‣ Execute script, correspondingly.

Mini-Language

Fityk comes with own domain-specific language (DSL), which is humbly called mini-language. All operations are driven by commands of the mini-language.

Do not worry

you do not need to learn the syntax of the mini-language. It is possible to use menus and dialogs in the GUI and avoid typing commands.

When you use the GUI and perform an action using the menu, you can see the corresponding command in the output window. It is one of less than 30 fityk commands. The commands have relatively simple syntax and perform single actions, such as loading data from file, adding function, assigning variable, fitting, or writing some values to a file.

It is possible to write a script (macro) as a sequence of such commands. This can automate common tasks, although some complex tasks still need to be programmed in a general-purpose language. That is why Fityk comes with embedded Lua (lightweight programming language) and also with bindings to Python and several other languages.

Now a quick glimpse at the syntax. Here, the =-> prompt marks an input:

=-> print pi
3.14159
=-> # this is a comment -- from `#' to the end of line
=-> p '2+3=', 2+3 # p stands for print
2+3 = 5
=-> set numeric_format='%.9f' # show 9 digits after dot
=-> pr pi, pi^2, pi^3 # pr, pri and prin also stand for print
3.141592654 9.869604401 31.006276680

Usually, one line is one command, but if it is really needed, two or more commands can be put in one line like this:

=-> $a = 3; $b = 5 # two commands separated with `;'

If the user works simultaneously with multiple datasets, she can refer to a dataset using its number: the first dataset is @0, the second – @1, etc:

=-> fit # perform fitting of the default dataset (the first one)
=-> @2: fit # fit the third dataset (@2)
=-> @*: fit # fit all datasets, one by one

All settings in the program are changed using the command set:

set key = value

For example:

=-> set logfile = 'C:\log.fit' # log all commands to this file
=-> set verbosity = 1 # make output from the program more verbose
=-> set epsilon = 1e-14

The last example changes the ε value, which is used to test floating-point numbers a and b for equality (it is well known that due to rounding errors the equality test for two numbers should have some tolerance, and the tolerance should be tailored to the application): |a−b| < ε.

To change a setting only for one command, add with key=value before the command:

=-> with epsilon = 0.1 print pi == 3.14 # abusing epsilon
1
=-> print pi == 3.14 # default epsilon = 10^-12
0

Writing informally, each line has a syntax:

[[@...:] [with ...] command [";" command]...] [#comment]

In scripts and in the CLI backslash () at the end of the line means that the next line is continuation.

All the commands are described in the next chapters.

Important

The rest of this section can be useful as reference, but it is recommended to skip it when reading the manual for the first time.

To keep the description above concise, some details were skipped.

The datasets listed before the colon (:) make a foreach loop:

=-> $a=0
=-> @0 @0 @0: $a={$a+1}; print $a
1
2
3

@* stands for all datasets, from @0 to the last one.

There is a small difference between two commands in two lines and two commands separated by ;. The whole line is parsed before the execution begins and some checks for the second command are performed before the first command is run:

=-> $a=4; print $a # print gives unexpected error
Error: undefined variable: $a

=-> $b=2
=-> $b=4; print $b # $b is already defined at the check time
4

Grammar

The grammar is expressed in EBNF-like notation:

  • (*this is a comment*)
  • A* means 0 or more occurrences of A.
  • A+ means 1 or more occurrences of A.
  • A % B means A (B A)* and the % operator has the highest precedence. For example: term % "+" comment is the same as term ("+" term)* comment.
  • The colon : in quoted string means that the string can be shortened, e.g. "del:ete" means that any of del, dele, delet and delete can be used.

The functions that can be used in p_expr and v_expr are available here and here, respectively. v_expr contains only a subset of functions from p_expr (partly, because we need to calculate symbolical derivatives of v_expr)

Line structure

line      ::=  [statement] [comment]
statement ::=  [Dataset+ ":"] [with_opts] command % ";"
with_opts ::=  "w:ith" (Lname "=" value) % ","
comment   ::=  "#" AllChars*

Commands

The kCmd* names in the comments correspond to constants in the code.

command ::=  (
             "deb:ug" RestOfLine              | (*kCmdDebug*)
             "def:ine" define                 | (*kCmdDefine*)
             "del:ete" delete                 | (*kCmdDelete*)
             "del:ete" delete_points          | (*kCmdDeleteP*)
             "e:xecute" exec                  | (*kCmdExec*)
             "f:it" fit                       | (*kCmdFit*)
             "g:uess" guess                   | (*kCmdGuess*)
             "i:nfo" info_arg % "," [redir]   | (*kCmdInfo*)
             "l:ua" RestOfLine                | (*kCmdLua*)
             "=" RestOfLine                   | (*kCmdLua*)
             "pl:ot" [range] [range] Dataset* [redir] | (*kCmdPlot*)
             "p:rint" print_args [redir]      | (*kCmdPrint*)
             "quit"                           | (*kCmdQuit*)
             "reset"                          | (*kCmdReset*)
             "s:et" (Lname "=" value) % ","   | (*kCmdSet*)
             "sleep" expr                     | (*kCmdSleep*)
             "title" "=" filename             | (*kCmdTitle*)
             "undef:ine" Uname % ","          | (*kCmdUndef*)
             "use" Dataset                    | (*kCmdUse*)
             "!" RestOfLine                   | (*kCmdShell*)
             Dataset "<" load_arg             | (*kCmdLoad*)
             Dataset "=" dataset_expr         | (*kCmdDatasetTr*)
             Funcname "=" func_rhs            | (*kCmdNameFunc*)
             param_lhs "=" v_expr             | (*kCmdAssignParam*)
             Varname "=" v_expr               | (*kCmdNameVar*)
             Varname "=" "copy" "(" var_id ")" | (*kCmdNameVar*)
             model_id ("="|"+=") model_rhs    | (*kCmdChangeModel*)
             (p_attr "[" expr "]" "=" p_expr) % "," | (*kCmdPointTr*)
             (p_attr "=" p_expr) % ","        | (*kCmdAllPointsTr*)
             "M" "=" expr                     ) (*kCmdResizeP*)

Other rules

define         ::=  Uname "(" (Lname [ "=" v_expr]) % "," ")" "="
                       ( v_expr |
                         component_func % "+" |
                         "x" "<" v_expr "?" component_func ":" component_func
                       )
component_func ::=  Uname "(" v_expr % "," ")"
delete         ::=  (Varname | func_id | Dataset | "file" filename) % ","
delete_points  ::=  "(" p_expr ")"
exec           ::=  filename |
                    "!" RestOfLine |
                    "=" RestOfLine
fit            ::=  [Number] [Dataset*] |
                    "undo" |
                    "redo" |
                    "history" Number |
                    "clear_history"
guess          ::=  [Funcname "="] Uname ["(" (Lname "=" v_expr) % "," ")"] [range]
info_arg       ::=  ...TODO
print_args     ::=  [("all" | ("if" p_expr ":")]
                    (p_expr | QuotedString | "title" | "filename") % ","
redir          ::=  (">"|">>") filename
value          ::=  (Lname | QuotedString | expr) (*value type depends on the option*)
model_rhs      ::=  "0" |
                    func_id |
                    func_rhs |
                    model_id |
                    "copy" "(" model_id ")"
func_rhs       ::=  Uname "(" ([Lname "="] v_expr) % "," ")" |
                    "copy" "(" func_id ")"
load_arg       ::=  filename Lname* |
                    "."
p_attr         ::=  ("X" | "Y" | "S" | "A")
model_id       ::=  [Dataset "."] ("F"|"Z")
func_id        ::=  Funcname |
                    model_id "[" Number "]"
param_lhs      ::=  Funcname "." Lname |
                    model_id "[" (Number | "*") "]" "." Lname
var_id         ::=  Varname |
                    func_id "." Lname
range          ::=  "[" [expr] ":" [expr] "]"
filename       ::=  QuotedString | NonblankString

Mathematical expressions

expr        ::=  expr_or ? expr_or : expr_or
expr_or     ::=  expr_and % "or"
expr_and    ::=  expr_not % "and"
expr_not    ::=  "not" expr_not | comparison
comparison  ::=  arith % ("<"|">"|"=="|">="|"<="|"!=")
arith       ::=  term % ("+"|"-")
term        ::=  factor % ("*"|"/")
factor      ::=  ('+'|'-') factor | power
power       ::=  atom ['**' factor]
atom        ::=  Number | "true" | "false" | "pi" |
                 math_func | braced_expr | ?others?
math_func   ::=  "sqrt" "(" expr ")" |
                 "gamma" "(" expr ")" |
                  ...
braced_expr ::=  "{" [Dataset+ ":"] p_expr "}"

The atom rule also accepts some fityk expressions, such as $variable, %function.parameter, %function(expr), etc.

p_expr and v_expr are similar to expr, but they use additional variables in the atom rule.

p_expr recognizes n, M, x, y, s, a, X, Y, S and A. All of them but n and M can be indexed (e.g. x[4]). Example: (x+x[n-1])/2.

v_expr uses all unknown names (Lname) as variables (example: a+b*x^2). Only a subset of functions (math_func) from expr is supported. The tilde (~) can be used to create simple-variables (~5), optionally with a domain in square brackets (~5[1:6]).

Since v_expr is used to define variables and user-defined functions, the program calculates symbolically derivatives of v_expr. That is why not all the function from expr are supported (they may be added in the future).

dataset_expr supports very limited set of operators and a few functions that take Dataset token as argument (example: @0 - shirley_bg(@0)).

Lexer

Below, some of the tokens produced by the fityk lexer are defined.

The lexer is context-dependend: NonblankString and RestOfLine are produced only when they are expected in the grammar.

Uname is used only for function types (Gaussian) and pseudo-parameters (%f.Area).

Dataset        ::=  "@"(Digit+|"+"|"*")
Varname        ::=  "$" Lname
Funcname       ::=  "%" Lname
QuotedString   ::=  "'" (AllChars - "'")* "'"
Lname          ::=  (LowerCase | "_") (LowerCase | Digit | "_")*
Uname          ::=  UpperCase AlphaNum+
Number         ::=  ?number read by strtod()?
NonblankString ::=  (AllChars - (WhiteSpace | ";" | "#" ))*
RestOfLine     ::=  AllChars*

Data

Loading Data

Data files are read using the xylib library.

In the GUI

click Load Data. If it just works for your files, you may go straight to Active and Inactive Points.

Points are loaded from files using the command:

dataslot < filename[:xcol:ycol:scol:block] [filetype options...]

where

  • dataslot should be replaced with @0, unless many datasets are to be used simultaneously (for details see: Working with Multiple Datasets),
  • xcol, ycol, scol (supported only in text file) are columns corresponding to x, y and std. dev. of y. Column 0 means index of the point: 0 for the first point, 1 for the second, etc.
  • block - selects one or more blocks of data from a multi-block file such as VAMAS
  • filetype usually can be omitted, because in most of the cases the filetype can be detected; the list of supported filetypes is at the end of this section
  • options depend on a filetype and usually are omitted

If the filename contains blank characters, a semicolon or comma, it should be put inside single quotation marks (together with colon-separated indices, if any).

Multiple y columns and/or blocks can be specified, see the examples below:

@0 < foo.vms
@0 < foo.fii text first_line_header
@0 < foo.dat:1:4:: # x,y - 1st and 4th columns
@0 < foo.dat:1:3,4:: # load two dataset (with y in columns 3,4)
@0 < foo.dat:1:3..5:: # load three dataset (with y in columns 3,4,5)
@0 < foo.dat:1:4..6,2:: # load four dataset (y: 4,5,6,2)
@0 < foo.dat:1:2..:: # load 2nd and all the next columns as y
@0 < foo.dat:1:2:3: # read std. dev. of y from 3rd column
@0 < foo.dat:0:1:: # x - 0,1,2,..., y - first column
@0 < 'foo.dat:0:1::' # the same
@0 < foo.raw::::0,1 # load two first blocks of data (as one dataset)

Information about loaded data can be obtained with:

info data

Supported Filetypes

text
ASCII text, multicolumn numeric data. The details are given in the next section.
dbws
format used by DBWS (program for Rietveld analysis) and DMPLOT.
cpi
Sietronics Sieray CPI format
uxd
Siemens/Bruker UXD format (powder diffraction data)
bruker_raw
Simens-Bruker RAW format (version 1,2,3)
rigaku_dat
Rigaku dat format (powder diffraction data)
vamas
VAMAS ISO-14976 (only experiment modes: “SEM” or “MAPSV” or “MAPSVDP” and only “REGULAR” scan mode are supported)
philips_udf
Philips UDF (powder diffraction data)
philips_rd
Philips RD raw scan format V3 (powder diffraction data)
spe
Princeton Instruments WinSpec SPE format (only 1-D data is supported)
pdcif
CIF for powder diffraction

The full list is available at: http://xylib.sourceforge.net/.

Reading Text Files

The xylib library can read TSV or CSV formats (tab or comma separated values). In fact, the values can be separated by any whitespace character or by one of ,;: punctations, or by any combination of these.

Empty lines and comments that start with hash (#) are skipped.

Since there is a lot of files in the world that contain numeric data mixed with text, unless the strict option is given any text that can not be interpreted as a number is regarded a start of comment (the rest of the line is ignored).

Note that the file is parsed regardless of blocks and columns specified by the user. The data read from the file are first stored in a table with m columns and n rows. If some of the lines have 3 numbers in it, and some have 5 numbers, we can either discard the lines that have 3 numbers or we can discard the numbers in 4th and 5th column. Usually the latter is done, but there are exceptions. The shorter lines are ignored

  • if it is the last line in the file (probably the program was terminated while writing the file),
  • if it contains only one number, but the prior lines had more numbers (this may be a comment that starts with a number)
  • if all the (not ignored) prior lines and the next line are longer

These rule were introduced to read free-format log files with textual comments inserted between lines with numeric data.

For now, xylib does not handle well nan’s and inf’s in the data.

Data blocks and columns may have names. These names are used to set a title of the dataset (see Working with Multiple Datasets for details). If the option first_line_header is given and the number of words in the first line is equal to the number of data columns, each word is used as a name of corresponding column. If the number of words is different, the first line is used as a name of the block. If the last_line_header option is given, the line preceding the first data line is used to set either column names or the block name.

If the file starts with the “LAMMPS (” string, the last_line_header option is set automatically. This is very helpful when plotting data from LAMMPS log files.

Active and Inactive Points

We often have the situation that only a part of the data from a file is of interest. In Fityk, each point is either active or inactive. Inactive points are excluded from fitting and all calculations. (Since active points do not need to be in one region, we do not use the region of interest term here, but such region can be easy selected). A data transformation:

A = boolean-condition

can be used to change the state of points.

In the GUI

data points can be activated and disactivated with mouse in the data-range mode (toolbar: Data-Range Mode).

Standard Deviation (or Weight)

When fitting data, we assume that only the y coordinate is subject to statistical errors in measurement. This is a common assumption. To see how the y‘s standard deviation, σ, influences fitting (optimization), look at the weighted sum of squared residuals formula in Nonlinear Optimization. We can also think about weights of points – every point has a weight assigned, that is equal w_i=1/\sigma_i^2.

Standard deviation of points can be read from file together with the x and y coordinates. Otherwise, it is set either to max(y1/2, 1) or to 1, depending on the default_sigma option. Setting std. dev. as a square root of the value is common and has theoretical ground when y is the number of independent events. You can always change the standard deviation, e.g. make it equal for every point with the command: S=1. See Data Point Transformations for details.

Note

It is often the case that user is not sure what standard deviation should be assumed, but it is her responsibility to pick something.

Data Point Transformations

Every data point has four properties: x coordinate, y coordinate, standard deviation of y and active/inactive flag. These properties can be changed using symbols X, Y, S and A, respectively. It is possible to either change a single point or apply a transformation to all points. For example:

  • Y[3]=1.2 assigns the y coordinate of the 4th point (0 is first),
  • Y = -y changes the sign of the y coordinate for all points.

On the left side of the equality sign you can have one of symbols X, Y, S, A, possibly with the index in brackets. The symbols on the left side are case insensitive.

The right hand side is a mathematical expression that can have special variables:

  • lower case letters x, y, s, a represent properties of data points before transformation,
  • upper case X, Y, S, A stand for the same properties after transformation,
  • M stands for the number of points.
  • n stands for the index of currently transformed point, e.g., Y=y[M-n-1] means that n-th point (n=0, 1, ... M-1) is assigned y value of the n-th point from the end.

Before the transformation a new array of points is created as a copy of the old array. Operations are applied sequentially from the first point to the last one, so while Y[n+1] and y[n+1] have always the same value, Y[n-1] and y[n-1] may differ. For example, the two commands:

Y = y[n] + y[n-1]
Y = y[n] + Y[n-1]

differ. The first one adds to each point the value of the previous point. The second one adds the value of the previous point after transformation, so effectively it adds the sum of all previous points. The index [n] could be omitted (Y = y + y[n-1]). The value of undefined points, like y[-1] and Y[-1], is explained later in this section.

Expressions can contain:

  • real numbers in normal or scientific format (e.g. 1.23e5),
  • constants pi, true (1), false (0)
  • binary operators: +, -, *, /, ^,
  • boolean operators: and, or, not,
  • comparisions: >, >=, <, <=, ==, !=.
  • one argument functions:
    • sqrt
    • exp
    • log10
    • ln
    • sin
    • cos
    • tan
    • sinh
    • cosh
    • tanh
    • atan
    • asin
    • acos
    • erf
    • erfc
    • gamma
    • lgamma (=ln(|gamma()|))
    • abs
    • round (rounds to the nearest integer)
  • two argument functions:
    • mod (modulo)
    • min2
    • max2 (max2(3,5) gives 5),
    • randuniform(a, b) (random number from interval (a, b)),
    • randnormal(mu, sigma) (random number from normal distribution),
    • voigt(a, b) = \frac{b}{\pi} \int_{-\infty}^{+\infty} \frac{\exp(-t^2)}{b^2+(a-t)^2} dt
  • ternary ?: operator: condition ?  expression1 : expression2, which returns expression1 if condition is true and expression2 otherwise.

A few examples.

  • The x scale of diffraction pattern can be changed from 2θ to Q:

    X = 4*pi * sin(x/2*pi/180) / 1.54051 # Cu 2θ -> Q
    
  • Negative y values can be zeroed:

    Y = max2(y, 0)
    
  • All standard deviations can be set to 1:

    S = 1
    
  • It is possible to select active range of data:

    A = x > 40 and x < 60 # select range (40, 60)
    

All operations are performed on real numbers. Two numbers that differ less than ε (the value of ε is set by the option epsilon) are considered equal.

Points can be created or deleted by changing the value of M. For example, the following commands:

M=500; x=n/100; y=sin(x)

create 500 points and generate a sinusoid.

Points are kept sorted according to their x coordinate. The sorting is performed after each transformation.

Note

Changing the x coordinate may change the order and indices of points.

Indices, like all other values, are computed in the real number domain. If the index is not integer (it is compared using ε to the rounded value):

  • x, y, s, a are interpolated linearly. For example, y[2.5] is equal to (y[2]+[3])/2. If the index is less than 0 or larger than M-1, the value for the first or the last point, respectively, is returned.
  • For X, Y, S, A the index is rounded to integer. If the index is less than 0 or larger than M-1, 0 is returned.

Transformations separated by commas (,) form a sequance of transformations. During the sequance, the vectors x, y, s and a that contain old values are not changed. This makes possible to swap the axes:

X=y, Y=x

The special index(arg) function returns the index of point that has x equal arg, or, if there is no such point, the linear interpolation of two neighbouring indices. This enables equilibrating the step of data (with interpolation of y and σ):

X = x[0] + n * (x[M-1]-x[0]) / (M-1), Y = y[index(X)], S = s[index(X)]

It is possible to delete points for which given condition is true, using expression delete(condition):

delete(not a) # delete inactive points

# reduce twice the number of points, averaging x and adding y
x = (x[n]+x[n+1])/2
y = y[n]+y[n+1]
delete(mod(n,2) == 1)

If you have more than one dataset, you may need to specify to which dataset the transformation applies. See Working with Multiple Datasets for details.

The value of a data expression can be shown using the print command. The precision of printed numbers is governed by the numeric_format option.

print M # the number of points
print y[index(20)] # value of y for x=20

Aggregate Functions

Aggregate functions have syntax:

aggregate(expression [if condition])

and return a single value, calculated from values of all points for which the given condition is true. If the condition is omitted, all points in the dataset are taken into account.

The following aggregate functions are recognized:

  • min() — the smallest value,

  • max() — the largest value,

  • argmin() — (stands for the argument of the minimum)

    the x value of the point for which the expression in brackets has the smallest value,

  • argmax() — the x value of the point for which the expression

    in brackets has the largest value,

  • sum() — the sum,

  • count() — the number of points for which the expression is true,

  • avg() — the arithmetic mean,

  • stddev() — the standard deviation,

  • darea() — a function used to normalize the area (see the example below). It returns the sum of expression*(x[n+1]-x[n-1])/2. In particular, darea(y) returns the interpolated area under data points.

Examples:

p avg(y) # print the average y value
p max(y) # the largest y value
p argmax(y) # the position of data maximum
p max(y if x > 40 and x < 60)   # the largest y value for x in (40, 60)
p max(y if a) # the largest y value in the active range
p min(x if y > 0.1)] # x of the first point with y > 0.1
p count(y>100) # the number of points that have y above 100
p count(y>avg(y)) # aggregate functions can be nested
p y[min(n if y > 100)] # the first (from the left) value of y above 100

# take the first 2000 points, average them and subtract as background
Y = y - avg(y if n<2000)

Y = y / darea(y) # normalize data area

# make active only the points on the left from the first
# point with y > 0.1
a = x < min(x if y > 0.1)]

Functions and Variables in Data Transformation

You may postpone reading this section and read about the Models first.

Variables ($foo) and functions (%bar) can be used in data expressions:

Y = y / $foo  # divides all y's by $foo
Y = y - %f(x) # subtracts function %f from data
Y = y - @0.F(x) # subtracts all functions in F

# print the abscissa value of the maximum of the model
# (only the values in points are considered,
#  so it's not exactly the model's maximum)
print argmax(F(x))

# print the maximum of the sum of two functions
print max(%_1(x) + %_2(x))

# Fit constant x-correction (i.e. fit the instrumental zero error), ...
Z = Constant(~0)
fit
X = x + Z(x)        # ... correct the data
Z = 0               # ... and remove the correction from the model.

In the GUI

in the Baseline Mode (Baseline Mode), functions Spline and Polyline are used to subtract manually selected background. Clicking Strip Background results in a command like this:

%bg0 = Spline(14.2979,62.1253, 39.5695,35.0676, 148.553,49.9493)
Y = y - %bg0(x)

Clicking the same button again undoes the subtraction:

Y = y + %bg0(x)

The function edited in the Baseline Mode is always named %bgX, where X is the index of the dataset.

Values of the function parameters (e.g. %fun.a0) and pseudo-parameters Center, Height, FWHM and Area (e.g. %fun.Area) can also be used. Pseudo-parameters are supported only by functions, which know how to calculate these properties.

It is also possible to calculate some properties of %functions:

  • %f.numarea(x1, x2, n) gives area integrated numerically from x1 to x2 using trapezoidal rule with n equal steps.
  • %f.findx(x1, x2, y) finds x in interval (x1, x2) such that %f(x)=y using bisection method combined with Newton-Raphson method. It is a requirement that %f(x1) < y < %f(x2).
  • %f.extremum(x1, x2) finds x in interval (x1, x2) such that %f’(x)=0 using bisection method. It is a requirement that %f’(x1) and %f’(x2) have different signs.

A few examples:

print %fun.numarea(, 0, 100, 10000) # shows area of function %fun
print %_1(%_1.extremum(40, 50)) # shows extremum value

# calculate FWHM numerically, value 50 can be tuned
$c = {%f.Center}
p %f.findx($c, $c+50, %f.Height/2) - %f.findx($c, $c-50, %f.Height/2)
p %f.FWHM # should give almost the same.

Working with Multiple Datasets

Let us call a set of data that usually comes from one file – a dataset. It is possible to work simultaneously with multiple datasets. Datasets have numbers and are referenced by @ with the number, (e.g. @3). The user can specify which dataset the command should be applied to:

@0: M=500    # change the number of points in the first dataset
@1 @2: M=500 # the same command applied to two datasets
@*: M=500    # and the same applied to all datasets

If the dataset is not specified, the command applies to the default dataset, which is initially @0. The use command changes the default dataset:

use @2 # set @2 as default

To load dataset from file, use one of the commands:

@n < filename:xcol:ycol:scol:block filetype options...

@+ < filename:xcol:ycol:scol:block filetype options...

The first one uses existing data slot and the second one creates a new slot. Using @+ increases the number of datasets, and the command delete @n decreases it.

The dataset can be duplicated (@+ = @n) or transformed, more on this in the next section.

Each dataset has a separate model, that can be fitted to the data. This is explained in the next chapter.

Each dataset also has a title (it does not have to be unique, however). When loading file, a title is automatically created:

  • if there is a name associated with the column ycol, the title is based on it;
  • otherwise, if there is a name associated with the data block read from file, the title is set to this name;
  • otherwise, the title is based on the filename

Titles can be changed using the command:

@n: title = 'new-title'

To print the title of the dataset, type @n: info title.

Dataset Transformations

There are a few transformations defined for a whole dataset or for two datasets. The syntax is @n = ... or @+ = .... The the right hand side expression supports the following operations:

-@n
negation of all y values,
d * @n
(e.g. 0.4*@0) y values are multiplied by d,
@n + @m
returns @n with added y values from interpolated @m,
@n - @m
returns @n with subtracted y values from interpolated @m,
@n and @m
returns points from both datasets (re-sorted),

and functions:

sum_same_x(@n)
Merges points which have distance in x is smaller than epsilon. x of the merged point is the average, and y and σ are sums of components.
avg_same_x(@n)
The same as sum_same_x, but y and σ are set as the average of components.
shirley_bg(@n)
Calculates Shirley background (useful in X-ray photoelectron spectroscopy).

Examples:

@+ = @0 # duplicate the dataset
@+ = @0 and @1 # create a new dataset from @0 and @1
@0 = @0 - shirley_bg(@0) # remove Shirley background
@0 = @0 - @1 # subtract @1 from @0
@0 = @0 - 0.28*@1 # subtract scaled dataset @1 from @0

Exporting Data

Command:

print all: expression, ... > file.tsv

can export data to an ASCII TSV (tab separated values) file.

In the GUI

Data ‣ Export

To export data in a 3-column (x, y and standard deviation) format, use:

print all: x, y, s > file.tsv

Any expressions can be printed out:

p all: n+1, x, y, F(x), y-F(x), %foo(x), sin(pi*x)+y^2 > file.tsv

It is possible to select which points are to be printed by replacing all with if followed by a condition:

print if a: x, y # only active points are printed
print if x > 30 and x < 40: x, y # only points in (30,40)

The option numeric_format controls the format and precision of all numbers.

Models

From Numerical Recipes, chapter 15.0:

Given a set of observations, one often wants to condense and summarize the data by fitting it to a “model” that depends on adjustable parameters. Sometimes the model is simply a convenient class of functions, such as polynomials or Gaussians, and the fit supplies the appropriate coefficients. Other times, the model’s parameters come from some underlying theory that the data are supposed to satisfy; examples are coefficients of rate equations in a complex network of chemical reactions, or orbital elements of a binary star. Modeling can also be used as a kind of constrained interpolation, where you want to extend a few data points into a continuous function, but with some underlying idea of what that function should look like.

This chapter shows how to construct the model.

Complex models are often a sum of many functions. That is why in Fityk the model F is constructed as a list of component functions and is computed as F = \sum_i f_i.

Each component function f_i is one of predefined functions, such as Gaussian or polynomial. This is not a limitation, because the user can add any function to the predefined functions.

To avoid confusion, the name function will be used only when referring to a component function, not when when referring to the sum (model), which mathematically is also a function. The predefined functions will be sometimes called function types.

Function f_i=f_i(x; \boldsymbol{a}) is a function of x, and depends on a vector of parameters \boldsymbol{a}. The parameters \boldsymbol{a} will be fitted to achieve agreement of the model and data.

In experiments we often have the situation that the measured x values are subject to systematic errors caused, for example, by instrumental zero shift or, in powder diffraction measurements, by displacement of sample in the instrument. If this is the case, such errors should be a part of the model. In Fityk, this part of the model is called x-correction. The final formula for the model is:

F(x; \boldsymbol{a}) = \sum_i f_i(x+Z(x; \boldsymbol{a}); \boldsymbol{a})

where Z(x; \boldsymbol{a}) = \sum_i z_i(x; \boldsymbol{a}) is the x-correction. Z is constructed as a list of components, analogously to F, although in practice it has rarely more than one component.

Each component function is created by specifying a function type and binding variables to type’s parameters. The next section explains what are variables in Fityk, and then we get back to functions.

Variables

Variables have names prefixed with the dollar symbol ($) and are created by assigning a value:

$foo=~5.3           # simple-variable
$bar=5*sin($foo)    # compound-variable
$c=3.1              # constant (the simplest compound-variable)

The numbers prefixed with the tilde (~) are adjustable when the model is fitted to the data. Variable created by assigning ~number (like $foo in the example above) will be called a simple-variable.

All other variables are called compound-variables. Compound variables either depend on other variables ($bar above) or are constant ($c).

Important

Unlike in popular programming languages, variable can store either a single numeric (floating-point) value or a mathematical expression. Nothing else. In case of expression, if we define $b=2*$a the value of $b will be recalculated every time $a changes.

To assign a value (constant) of another variable, use: $b={$a}. Braces return the current value of the enclosed expression. The left brace can be preceded by the tilde (~). The assignment $b=~{$a} creates a simple variable.

Compound-variables can be build using operators +, -, *, /, ^ and the functions sqrt, exp, log10, ln, sin, cos, tan, sinh, cosh, tanh, atan, asin, acos, erf, erfc, lgamma, abs, voigt. This is a subset of the functions used in data transformations.

The braces may contain any data expression:

$x0 = {x[0]}
$min_y = {min(y if a)}
$c = {max2($a, $b)}
$t = {max(x) < 78 ? $a : $b}

Sometimes it is useful to freeze a variable, i.e. to prevent it from changing while fitting:

$a = ~12.3 # $a is fittable (simple-variable)
$a = {$a}  # $a is not fittable (constant)
$a = ~{$a} # $a is fittable (simple-variable) again

In the GUI

a variable can be switched between constant and simple-variable by clicking the padlock button on the sidebar. The icons open lock and lock show that the variable is fittable and frozen, respectively.

If the assigned expression contains tildes:

$bleh=~9.1*exp(~2)

it automatically creates simple-variables corresponding to the tilde-prefixed numbers. In the example above two simple-variables (with values 9.1 and 2) are created. Automatically created variables are named $_1, $_2, $_3, and so on.

Variables can be deleted using the command:

delete $variable

Domains

Simple-variables may have a domain, which is used for two things when fitting.

Most importantly, fitting methods that support bound constraints use the domain as lower and/or upper bounds. See the section Bound Constraints for details.

The other use is for randomizing parameters (simple-variables) of the model. Methods that stochastically initialize or modify parameters (usually generating a set of initial points) need well-defined domains (minimum and maximum values for parameters) to work effectively. Such methods include Nelder-Mead simplex and Genetic Algorithms, but not the default Lev-Mar method, so in most cases you do not need to worry about it.

The syntax is as follows:

$a = ~12.3 [0:20] # initial values are drawn from the (0, 20) range
$a = ~12.3 [0:]   # only lower bound
$a = ~12.3 [:20]  # only upper bound
$a = ~15.0        # domain stays the same
$a = ~15.0 []     # no domain
$a = ~{$a} [0:20] # domain is set again

If the domain is not specified but it is required (for the latter use) by the fitting method, we assume it to be ±p% of the current value, where p can be set using the domain_percent option.

Function Types and Functions

Function types have names that start with upper case letter, e.g. Linear or Voigt.

Functions have names prefixed with the percent symbol, e.g. %func. Every function has a type and variables bound to its parameters.

Functions can be created by giving the type and the correct number of variables in brackets, e.g.:

%f1 = Gaussian(~66254., ~24.7, ~0.264)
%f2 = Gaussian(~6e4, $ctr, $b+$c)
%f3 = Gaussian(height=~66254., hwhm=~0.264, center=~24.7)

Every expression which is valid on the right-hand side of a variable assignment can be used as a variable. If it is not just a name of a variable, an automatic variable is created. In the above examples, two variables were implicitely created for %f2: first for value 6e4 and the second for $b+$c).

If the names of function’s parameters are given (like for %f3), the variables can be given in any order.

Function types can can have specified default values for some parameters. The variables for such parameters can be omitted, e.g.:

=-> i Pearson7
Pearson7(height, center, hwhm, shape=2) = height/(1+((x-center)/hwhm)^2*(2^(1/shape)-1))^shape
=-> %f4 = Pearson7(height=~66254., center=~24.7, fwhm=~0.264) # no shape is given
New function %f4 was created.

Functions can be copied. The following command creates a deep copy (i.e. all variables are also duplicated) of %foo:

%bar = copy(%foo)

Functions can be also created with the command guess, as described in Guessing Initial Parameters.

Variables bound to the function parameters can be changed at any time:

=-> %f = Pearson7(height=~66254., center=~24.7, fwhm=~0.264)
New function %f was created.
=-> %f.center=~24.8
=-> $h = ~66254
=-> %f.height=$h
=-> info %f
%f = Pearson7($h, $_5, $_3, $_4)
=-> $h = ~60000 # variables are kept by name, so this also changes %f
=-> %p1.center = %p2.center + 3 # keep fixed distance between %p1 and %p2

Functions can be deleted using the command:

delete %function

Built-In Functions

The list of all functions can be obtained using i types. Some formulae here have long parameter names (like “height”, “center” and “hwhm”) replaced with a_i

Gaussian:

y = a_0
    \exp\left[-\ln(2)\left(\frac{x-a_1}{a_2}\right)^{2}\right]

a_2 here is half width at half maximum (HWHM=FWHM/2, where FWHM stands for full width...), which is proportional to standard deviation: a_2=\sqrt{2\ln2}\sigma.

SplitGaussian:

y(x;a_0,a_1,a_2,a_3) = \begin{cases}
\textrm{Gaussian}(x;a_0,a_1,a_2) & x\leq a_1\\
\textrm{Gaussian}(x;a_0,a_1,a_3) & x>a_1\end{cases}

GaussianA:

y = \sqrt{\frac{\ln(2)}{\pi}}\frac{a_0}{a_2}
    \exp\left[-\ln(2)\left(\frac{x-a_1}{a_2}\right)^{2}\right]

Lorentzian:

y = \frac{a_0}{1+\left(\frac{x-a_1}{a_2}\right)^2}

SplitLorentzian:

y(x;a_0,a_1,a_2,a_3) = \begin{cases}
\textrm{Lorentzian}(x;a_0,a_1,a_2) & x\leq a_1\\
\textrm{Lorentzian}(x;a_0,a_1,a_3) & x>a_1\end{cases}

LorentzianA:

y = \frac{a_0}{\pi a_2\left[1+\left(\frac{x-a_1}{a_2}\right)^2\right]}

Pearson VII (Pearson7):

y = \frac{a_0} {\left[1+\left(\frac{x-a_1}{a_2}\right)^2
                        \left(2^{\frac{1}{a_3}}-1\right)\right]^{a_3}}

split Pearson VII (SplitPearson7):

y(x;a_{0},a_{1},a_{2},a_{3},a_{4},a_{5}) = \begin{cases}
 \textrm{Pearson7}(x;a_0,a_1,a_2,a_4) & x\leq a_1\\
 \textrm{Pearson7}(x;a_0,a_1,a_3,a_5) & x>a_1\end{cases}

Pearson VII Area (Pearson7A):

y = \frac{a_0\Gamma(a_3)\sqrt{2^{\frac{1}{a_3}}-1}}
         {a_2\Gamma(a_3-\frac{1}{2})\sqrt{\pi} \left[
            1 + \left(\frac{x-a_1}{a_2}\right)^2
                \left(2^{\frac{1}{a_3}}-1\right)
         \right]^{a_3}}

Pseudo-Voigt (PseudoVoigt):

y = a_0 \left[(1-a_3)\exp\left(-\ln(2)\left(\frac{x-a_1}{a_2}\right)^2\right)
              + \frac{a_3}{1+\left(\frac{x-a_1}{a_2}\right)^2}
        \right]

Pseudo-Voigt is a name given to the sum of Gaussian and Lorentzian. a_3 parameters in Pearson VII and Pseudo-Voigt are not related.

split Pseudo-Voigt (SplitPseudoVoigt):

y(x;a_{0},a_{1},a_{2},a_{3},a_{4},a_{5}) = \begin{cases}
 \textrm{PseudoVoigt}(x;a_0,a_1,a_2,a_4) & x\leq a_1\\
 \textrm{PseudoVoigt}(x;a_0,a_1,a_3,a_5) & x>a_1\end{cases}

Pseudo-Voigt Area (PseudoVoigtA):

y = a_0 \left[\frac{(1-a_3)\sqrt{\ln(2)}}{a_2\sqrt{\pi}}
              \exp\left(-\ln2\left(\frac{x-a_1}{a_2}\right)^2\right)
              + \frac{a_3}{\pi a_2
                           \left[1+\left(\frac{x-a_1}{a_2}\right)^2\right]}
        \right]

Voigt:

y = \frac
    {a_0 \int_{-\infty}^{+\infty}
             \frac{\exp(-t^2)}{a_3^2+(\frac{x-a_1}{a_2}-t)^2} dt}
    {\int_{-\infty}^{+\infty}
             \frac{\exp(-t^2)}{a_3^2+t^2} dt}

The Voigt function is a convolution of Gaussian and Lorentzian functions. a_0 = heigth, a_1 = center, a_2 is proportional to the Gaussian width, and a_3 is proportional to the ratio of Lorentzian and Gaussian widths.

Voigt is computed according to R.J.Wells, Rapid approximation to the Voigt/Faddeeva function and its derivatives, Journal of Quantitative Spectroscopy & Radiative Transfer 62 (1999) 29-48. The approximation is very fast, but not very exact.

FWHM is estimated using an approximation called modified Whiting (Olivero and Longbothum, 1977, JQSRT 17, 233): 0.5346 w_L + \sqrt{0.2169 w_L^2 + w_G^2}, where w_G=2\sqrt{\ln(2)} |a_2|, w_L=2 |a_2| a_3.

VoigtA:

y = \frac{a_0}{\sqrt{\pi}a_2}
    \int_{-\infty}^{+\infty}
        \frac{\exp(-t^2)}{a_3^2+(\frac{x-a_1}{a_2}-t)^2} dt

split Voigt (SplitVoigt):

y(x;a_{0},a_{1},a_{2},a_{3},a_{4},a_{5}) = \begin{cases}
 \textrm{Voigt}(x;a_0,a_1,a_2,a_4) & x\leq a_1\\
 \textrm{Voigt}(x;a_0,a_1,a_3,a_5) & x>a_1\end{cases}

Exponentially Modified Gaussian (EMG):

y = \frac{ac\sqrt{2\pi}}{2d}
    \exp\left(\frac{c^2}{2d^2}-\frac{x-b}{d}\right)
    \left[\frac{d}{\left|d\right|}
          +\textrm{erf}\left(\frac{x-b}{\sqrt{2}c}
                             - \frac{c}{\sqrt{2}d}\right)
    \right]

The exponentially modified Gaussian is a convolution of Gaussian and exponential probability density. a = Gaussian heigth, b = location parameter (Gaussian center), c = Gaussian width, d = distortion parameter (a.k.a. modification factor or time constant).

LogNormal:

y = h \exp\left\{ -\ln(2) \left[
                                \frac{\ln\left(1+2b\frac{x-c}{w}\right)}{b}
                         \right]^{2} \right\}

Doniach-Sunjic (DoniachSunjic):

y = \frac{h\left[\frac{\pi a}{2}
                 + (1-a)\arctan\left(\frac{x-E}{F}\right)\right]}
         {F+(x-E)^2}

Polynomial5:

y = a_0 + a_1 x +a_2 x^2 + a_3 x^3 + a_4 x^4 + a_5 x^5

Sigmoid:

y = L + \frac{U-L}{1+\exp\left(-\frac{x-x_{mid}}{w}\right)}

FCJAsymm:

Axial asymmetry peak shape in the Finger, Cox and Jephcoat model, see J. Appl. Cryst. (1994) 27, 892 and J. Appl. Cryst. (2013) 46, 1219.

Variadic Functions

Variadic function types have variable number of parameters. Two variadic function types are defined:

Spline(x1, y1, x2, y2, ...)
Polyline(x1, y1, x2, y2, ...)

This example:

%f = Spline(22.1, 37.9, 48.1, 17.2, 93.0, 20.7)

creates a function that is a natural cubic spline interpolation through points (22.1, 37.9), (48.1, 17.2), ....

The Polyline function is a polyline interpolation (spline of order 1).

Both Spline and Polyline functions are primarily used for the manual baseline subtraction via the GUI.

The derivatives of Spline function are not calculated, so this function is not refined by the default, derivative-based fitting algorithm.

Since the Polyline derivatives are calculated, it is possible to perform weighted least squares approximation by broken lines, although non-linear fitting algorithms are not optimal for this task.

User-Defined Functions (UDF)

User-defined function types can be added using command define, and then used in the same way as built-in functions.

Example:

define MyGaussian(height, center, hwhm) = height*exp(-ln(2)*((x-center)/hwhm)^2)
  • The name of new type must start with an upper-case letter, contain only letters and digits and have at least two characters.

  • The name of the type is followed by parameters in brackets.

  • Parameter name must start with lowercase letter and, contain only lowercase letters, digits and the underscore (‘_’).

  • The name “x” is reserved, do not put it into parameter list, just use it on the right-hand side of the definition.

  • There are special names of parameters that Fityk understands:

    • if the functions is peak-like (bell-shaped): height, center, hwhm, area,
    • if the functions is S-shaped (sigmoidal) or step-like: lower, upper, xmid, wsig,
    • if the function is more like linear: slope, intercept, avgy.

    The initial values of these parameters can be guessed (command guess) from the data. hwhm means half width at half maximum, the other names are self-explaining.

  • Each parameter may have a default value (see the examples below). The default value can be either a number or an expression that depends on the parameters listed above (e.g. 0.8*hwhm). The default value always binds a simple-variable to the parameter.

UDFs can be defined in a few ways:

  • by giving a full formula, like in the example above,
  • as a re-parametrization of existing function (see the GaussianArea example below),
  • as a sum of already defined functions (see the GLSum example below),
  • as a splitted (bifurcated) function: x < expression ? Function1(...) : Function2(...) (see the SplitL example below).

When giving a full formula, the right-hand side of the equality sign is similar to the definiton of variable, but the formula can also depend on x. Hopefully the examples can make the syntax clear:

# this is how some built-in functions could be defined
define MyGaussian(height, center, hwhm) = height*exp(-ln(2)*((x-center)/hwhm)^2)
define MyLorentzian(height, center, hwhm) = height/(1+((x-center)/hwhm)^2)
define MyCubic(a0=height,a1=0, a2=0, a3=0) = a0 + a1*x + a2*x^2 + a3*x^3

# supersonic beam arrival time distribution
define SuBeArTiDi(c, s, v0, dv) = c*(s/x)^3*exp(-(((s/x)-v0)/dv)^2)/x

# area-based Gaussian can be defined as modification of built-in Gaussian
# (it is the same as built-in GaussianA function)
define GaussianArea(area, center, hwhm) = Gaussian(area/hwhm/sqrt(pi/ln(2)), center, hwhm)

# sum of Gaussian and Lorentzian, a.k.a. PseudoVoigt (should be in one line)
define GLSum(height, center, hwhm, shape) = Gaussian(height*(1-shape), center, hwhm)
+ Lorentzian(height*shape, center, hwhm)

# split-Gaussian, the same as built-in SplitGaussian (should be in one line)
define SplitG(height, center, hwhm1=fwhm*0.5, hwhm2=fwhm*0.5) =
  x < center ? Lorentzian(height, center, hwhm1)
             : Lorentzian(height, center, hwhm2)

There is a simple substitution mechanism that makes writing complicated functions easier. Substitutions must be assigned in the same line, after the keyword where.

Example:

define ReadShockley(sigma0=1, a=1) = sigma0 * t * (a - ln(t)) where t=x*pi/180

# more complicated example, with nested substitutions
define FullGBE(k, alpha) = k * alpha * eta * (eta / tanh(eta) - ln (2*sinh(eta))) where eta = 2*pi/alpha * sin(theta/2), theta=x*pi/180

How it works internally

The formula is parsed, derivatives of the formula are calculated symbolically, expressions are simplified and bytecode for virtual machine (VM) is created.

When fitting, the VM calculates the value of the function and derivatives for every point.

Defined functions can be undefined using command undefine:

undefine GaussianArea

It is common to add own definitions to the init file. See the section Starting fityk and cfityk for details.

Cutoff

With default settings, the value of every function is calculated at every point. Peak functions, such as Gaussian, often have non-negligible values only in a small fraction of all points, so if you have many narrow peaks (like here), the basic optimization is to calculate values of each peak function only near the function’s center. If the option function_cutoff is set to a non-zero value, each function is evaluated only in the range where its values are greater than the function_cutoff.

This optimization is supported only by some built-in functions.

Model, F and Z

As already discussed, each dataset has a separate model that can be fitted to the data. As can be seen from the formula at the beginning of this chapter, the model is defined as a set functions f_i and a set of functions z_i. These sets are named F and Z respectively. The model is constructed by specifying names of functions in these two sets.

In many cases x-correction Z is not used. The fitted curve is thus the sum of all functions in F.

Command:

F += %function

adds %function to F, and

Z += %function

adds %function to Z.

A few examples:

# create and add function to F
%g = Gaussian(height=~66254., hwhm=~0.264, center=~24.7)
F += %g

# create unnamed function and add it to F
F += Gaussian(height=~66254., hwhm=~0.264, center=~24.7)

# clear F
F = 0

# clear F and put three functions in it
F = %a + %b + %c

# show info about the first and the last function in F
info F[0], F[-1]

The next sections shows an easier way to add a function (command guess).

If there is more than one dataset, F and Z can be prefixed with the dataset number (e.g. @1.F).

The model can be copied. To copy the model from @0 to @1 we type one of the two commands:

@1.F = @0.F        # shallow copy
@1.F = copy(@0.F)  # deep copy

The former command uses the same functions in both models: if you shift a peak in @1, it will be also shifted in @0. The latter command (deep copy) duplicates all functions and variables and makes an independent model.

In the GUI

click the button Copy-Model on the sidebar to make a deep copy.

It is often required to keep the width or shape of peaks constant for all peaks in the dataset. To change the variables bound to parameters with a given name for all functions in F, use the command:

F[*].param = variable

Examples:

# Set hwhm of all functions in F that have a parameter hwhm to $foo
# (hwhm here means half-width-at-half-maximum)
F[*].hwhm = $foo

# Bound the variable used for the shape of peak %_1 to shapes of all
# functions in F
F[*].shape = %_1.shape

# Create a new simple-variable for each function in F and bound the
# variable to parameter hwhm. All hwhm parameters will be independent.
F[*].hwhm = ~0.2

In the GUI

buttons =W and =S on the sidebar make, respectively, the HWHM and shape of all functions the same. Pressing the buttons again will make all the parameters independent.

Guessing Initial Parameters

The program can automatically set initial parameters of peaks (using peak-detection algorithm) and lines (using linear regression). Choosing initial parameters of a function by the program will be called guessing.

It is possible to guess peak location and add it to F with the command:

guess [%name =] PeakType [(initial values...)] [[x1:x2]]

Examples:

# add Gaussian in the given range
@0: guess Gaussian [22.1:30.5]

# the same, but name the new function %f1
@0: guess %f1 = Gaussian [22.1:30.5]

# search for the peak in the whole dataset
@0: guess Gaussian

# add one Gaussian to each dataset
@*: guess Gaussian

# set the center and shape explicitely (determine height and width)
guess PseudoVoigt(center=$ctr, shape=~0.3) [22.1:30.5]
  • Name of the function is optional.
  • Some of the parameters can be specified in brackets.
  • If the range is omitted, the whole dataset will be searched.

Fityk offers a simple algorithm for peak-detection. It finds the highest point in the given range (center and height), and than tries to find the width of the peak (hwhm, and area = height × hwhm).

If the highest point is at boundary of the given range, the points from the boundary to the nearest local minimum are ignored.

The values of height and width found by the algorithm are multiplied by the values of options height_correction and width_correction, respectively. The default value for both options is 1.

Another simple algorithm can roughly estimate initial parameters of sigmoidal functions.

The linear traits slope and intercept are calculated using linear regression (without weights of points). avgy is calculated as average value of y.

In the GUI

select a function from the list of functions on the toolbar and press Auto Add to add (guess) the selected function.

To choose a data range change the GUI mode to Add-Peak Mode and select the range with the right mouse button.

Displaying Information

The info command can be show useful information when constructing the model.

info types
shows the list of available function types.
info FunctionType
(e.g. info Pearson7) shows the formula (definition).
info guess [range]
shows where the guess command would locate a peak.
info functions
lists all defined functions.
info variables
lists all defined variables.
info F
lists components of F.
info Z
lists components of Z.
info formula
shows the full mathematical formula of the fitted model.
info simplified_formula
shows the same, but the formula is simplified.
info gnuplot_formula
shows same as formula, but the output is readable by gnuplot, e.g. x^2 is replaced by x**2.
info simplified_gnuplot_formula
shows the simplified formula in the gnuplot format.
info peaks
show a formatted list of parameters of functions in F.
info peaks_err
shows the same data, additionally including uncertainties of the parameters.
info models
a script that reconstructs all variables, functions and models.

The last two commands are often redirected to a file (info peaks > filename).

The complete list of info arguments can be found in Information Display.

In the GUI

most of the above commands has clickable equivalents.

Curve Fitting

Nonlinear Optimization

This is the core. We have a set of observations (data points) to which we want to fit a model that depends on adjustable parameters. Let me quote Numerical Recipes, chapter 15.0, page 656):

The basic approach in all cases is usually the same: You choose or design a figure-of-merit function (merit function, for short) that measures the agreement between the data and the model with a particular choice of parameters. The merit function is conventionally arranged so that small values represent close agreement. The parameters of the model are then adjusted to achieve a minimum in the merit function, yielding best-fit parameters. The adjustment process is thus a problem in minimization in many dimensions. [...] however, there exist special, more efficient, methods that are specific to modeling, and we will discuss these in this chapter. There are important issues that go beyond the mere finding of best-fit parameters. Data are generally not exact. They are subject to measurement errors (called noise in the context of signal-processing). Thus, typical data never exactly fit the model that is being used, even when that model is correct. We need the means to assess whether or not the model is appropriate, that is, we need to test the goodness-of-fit against some useful statistical standard. We usually also need to know the accuracy with which parameters are determined by the data set. In other words, we need to know the likely errors of the best-fit parameters. Finally, it is not uncommon in fitting data to discover that the merit function is not unimodal, with a single minimum. In some cases, we may be interested in global rather than local questions. Not, “how good is this fit?” but rather, “how sure am I that there is not a very much better fit in some corner of parameter space?”

Our function of merit is the weighted sum of squared residuals (WSSR), also called chi-square:

\chi^{2}(\mathbf{a})
  =\sum_{i=1}^{N} \left[\frac{y_i-y(x_i;\mathbf{a})}{\sigma_i}\right]^{2}
  =\sum_{i=1}^{N} w_{i}\left[y_{i}-y(x_{i};\mathbf{a})\right]^{2}

Weights are based on standard deviations, w_i=1/\sigma_i^2. You can learn why squares of residuals are minimized e.g. from chapter 15.1 of Numerical Recipes.

The most popular method for curve-fitting is Levenberg-Marquardt. Fityk can also use a few general-purpose optimization methods. These methods are slower, some of them are orders of magnitude slower. Sometimes alternative methods find global minimum when the L-M algorithm is stuck in a local minimum, but in majority of cases the default L-M method is superior.

Uncertainty of Parameters

(It is easier for me to find a quote than to express it myself).

From the book J. Wolberg, Data Analysis Using the Method of Least Squares: Extracting the Most Information from Experiments, Springer, 2006, p.50:

(...) we turn to the task of determining the uncertainties associated with the a_k‘s. The usual measures of uncertainty are standard deviation (i.e., σ) or variance (i.e., σ2) so we seek an expression that allows us to estimate the \sigma_{a_k}‘s. It can be shown (...) that the following expression gives us an unbiased estimate of \sigma_{a_k}:

\sigma_{a_k}^{2}=\frac{S}{n-p}C_{kk}^{-1}

Note that \sigma_{a_k} is a square root of the value above. In this formula n-p, the number of (active) data points minus the number of independent parameters, is equal to the number of degrees of freedom. S is another symbol for \chi^2.

Terms of the C matrix are given as (p. 47 in the same book):

C_{jk}=\sum_{i=1}^n w_i \frac{\partial f}{\partial a_j} \frac{\partial f}{\partial a_k}

\sigma_{a_k} above is often called a standard error. Having standard errors, it is easy to calculate confidence intervals. But all these values should be used with care. Now another book will be cited: H. Motulsky and A. Christopoulos, Fitting Models to Biological Data Using Linear and Nonlinear Regression: A Practical Guide to Curve Fitting, Oxford University Press, 2004. This book can be downloaded for free as a manual to GraphPad Prism 4.

The standard errors reported by most nonlinear regression programs (...) are “approximate” or “asymptotic”. Accordingly, the confidence intervals computed using these errors should also be considered approximate.

It would be a mistake to assume that the “95% confidence intervals” reported by nonlinear regression have exactly a 95% chance of enclosing the true parameter values. The chance that the true value of the parameter is within the reported confidence interval may not be exactly 95%. Even so, the asymptotic confidence intervals will give you a good sense of how precisely you have determined the value of the parameter.

The calculations only work if nonlinear regression has converged on a sensible fit. If the regression converged on a false minimum, then the sum-of-squares as well as the parameter values will be wrong, so the reported standard error and confidence intervals won’t be helpful.

Bound Constraints

Simple-variables can have a domain. Fitting method mpfit (Lev-Mar implementation) and the methods from the NLOpt library use domains to constrain the parameters – they never let the parameters go outside of the domain during fitting.

In the literature, bound constraints are also called box constraints or, more generally, inequality constraints. Now a quotation discouraging the use of constraints. Peter Gans, Data Fitting in the Chemical Sciences by the Method of Least Squares, John Wiley & Sons, 1992, chapter 5.2.2:

Before looking at ways of dealing with inequality constraints we must ask a fundamental question: are they necessary? In the physical sciences and in least-squares minimizations in particular, inequality constraints are not always justified. The most common inequality constraint is that some number that relates to a physical quantity should be positive, pj > 0. If an unconstrained minimalization leads to a negative value, what are we to conclude? There are three possibilities; (a) the refinement has converged to a false minimum; (b) the model is wrong; (c) the parameter is not well defined by the data and is not significantly different from zero. In each of these three cases a remedy is at hand that does not involve constrained minimization: (a) start the refinement from good first estimates of the parameters; (b) change the model; (c) improve the quality of the data by further experimental work. If none of these remedies cure the problem of non-negativity constraints, then something is seriously wrong with the patient, and constrained minimization will probably not help.

Setting the domain is described in the section Domains.

For a convenience, the box_constraints option can globally disable (and re-enable) the constraints.

Levenberg-Marquardt

This is a standard nonlinear least-squares routine, and involves computing the first derivatives of functions. For a description of the algorithm see Numerical Recipes, chapter 15.5 or Siegmund Brandt, Data Analysis, chapter 10.15. Essentially, it combines an inverse-Hessian method with a steepest descent method by introducing a λ factor. When λ is equal to 0, the method is equivalent to the inverse-Hessian method. When λ increases, the shift vector is rotated toward the direction of steepest descent and the length of the shift vector decreases. (The shift vector is a vector that is added to the parameter vector.) If a better fit is found on iteration, λ is decreased.

Two implementation of this method are available: one from the MPFIT library, based on the old good MINPACK code (default method since ver. 1.3.0), and a custom implementation (default method in earlier fityk versions).

To switch between the two implementation use command:

set fitting_method = mpfit               # switch to MPFIT
set fitting_method = levenberg_marquardt # switch to fityk implem. of L-M

The following stopping criteria are available for mpfit:

  • the relative change of WSSR is smaller than the value of the ftol_rel option (default: 10^-10),
  • the relative change of parameters is smaller than the value of the xtol_rel option (default: 10^-10),

and for levenberg_marquardt:

  • the relative change of WSSR is smaller than the value of the lm_stop_rel_change option twice in row,
  • λ is greater than the value of the lm_max_lambda option (default: 10^15), which normally means WSSR is not changing due to limited numerical precision.

Nelder-Mead Downhill Simplex

To quote chapter 4.8.3, p. 86 of Peter Gans, Data Fitting in the Chemical Sciences by the Method of Least Squares:

A simplex is a geometrical entity that has n+1 vertices corresponding to variations in n parameters. For two parameters the simplex is a triangle, for three parameters the simplex is a tetrahedron and so forth. The value of the objective function is calculated at each of the vertices. An iteration consists of the following process. Locate the vertex with the highest value of the objective function and replace this vertex by one lying on the line between it and the centroid of the other vertices. Four possible replacements can be considered, which I call contraction, short reflection, reflection and expansion.[...] It starts with an arbitrary simplex. Neither the shape nor position of this are critically important, except insofar as it may determine which one of a set of multiple minima will be reached. The simplex than expands and contracts as required in order to locate a valley if one exists. Then the size and shape of the simplex is adjusted so that progress may be made towards the minimum. Note particularly that if a pair of parameters are highly correlated, both will be simultaneously adjusted in about the correct proportion, as the shape of the simplex is adapted to the local contours.[...] Unfortunately it does not provide estimates of the parameter errors, etc. It is therefore to be recommended as a method for obtaining initial parameter estimates that can be used in the standard least squares method.

This method is also described in previously mentioned Numerical Recipes (chapter 10.4) and Data Analysis (chapter 10.8).

There are a few options for tuning this method. One of these is a stopping criterium nm_convergence. If the value of the expression 2(Mm)/(M+m), where M and m are the values of the worst and best vertices respectively (values of objective functions of vertices, to be precise!), is smaller then the value of nm_convergence option, fitting is stopped. In other words, fitting is stopped if all vertices are almost at the same level.

The remaining options are related to initialization of the simplex. Before starting iterations, we have to choose a set of points in space of the parameters, called vertices. Unless the option nm_move_all is set, one of these points will be the current point – values that parameters have at this moment. All but this one are drawn as follows: each parameter of each vertex is drawn separately. It is drawn from a distribution that has its center in the center of the domain of the parameter, and a width proportional to both width of the domain and value of the nm_move_factor parameter. Distribution shape can be set using the option nm_distribution as one of: uniform, gaussian, lorentzian and bound. The last one causes the value of the parameter to be either the greatest or smallest value in the domain of the parameter – one of the two bounds of the domain (assuming that nm_move_factor is equal 1).

NLopt

A few methods from the NLopt library are available:

  • nlopt_nm – Nelder-Mead method, similar to the one described above,
  • nlopt_lbfgs – low-storage BFGS,
  • nlopt_var2 – shifted limited-memory variable-metric,
  • nlopt_praxis – PRAXIS (PRincipal AXIS),
  • nlopt_bobyqa – BOBYQA,
  • nlopt_sbplx – Sbplx (based on Subplex),

All NLopt methods have the same stopping criteria (in addition to the common criteria):

  • an optimization step changes the WSSR value by less than the value of the ftol_rel option (default: 10^-10) multiplied by the WSSR,
  • an optimization step changes every parameter by less than the value of the xtol_rel option (default: 10^-10) multiplied by the absolute value of the parameter.

Scripts

Working with Scripts

Scripts can be executed using the command:

exec filename

The file can be either a fityk script (usually with extension fit), or a Lua script (extension lua).

Note

Fityk can save its state to a script (info state > file.fit). It can also save all commands executed (directly or via GUI) in the session to a script (info history > file.fit).

Since a Fityk script with all the data inside can be a large file, the files may be stored compressed and it is possible to directly read gzip-compressed fityk script (.fit.gz).

Embedded Lua interpreter can execute any program in Lua 5.1. One-liners can be run with command lua:

=-> lua print(_VERSION)
Lua 5.1
=-> lua print(os.date("Today is %A."))
Today is Thursday.
=-> lua for n,f in F:all_functions() do print(n, f, f:get_template_name()) end
0       %_1     Constant
1       %_2     Cycle

(The Lua print function in fityk is redefined to show the output in the GUI instead of writing to stdout).

Like in the Lua interpreter, = at the beginning of line can be used to save some typing:

=-> = os.date("Today is %A.")
Today is Thursday.

Similarly, exec= also interprets the rest of line as Lua expressions, but this time the resulting string is executed as a fityk command:

=-> = string.format("fit @%d", math.random(0,5))
fit 17
=-> exec= string.format("fit @%d", math.random(0,5))
# fitting random dataset (useless example)

The Lua interpreter in Fityk has defined global object F which enables interaction with the program:

=-> = F:get_info("version")
Fityk 1.2.1

Now the first example that can be useful. For each dataset write output of the info peaks command to a file named after the data file, with appended ”.out”:

=-> @*: lua F:execute("info peaks >'%s.out'" % F:get_info("filename"))

This and other methods of F are documented in the next section.

Here is an example of Lua-Fityk interactions:

-- load data from files file01.dat, file02.dat, ... file13.dat
for i=1,13 do
    F:execute("@+ < file%02d.dat:0:1::" % i)
end

-- print some statistics about the loaded data
n = F:get_dataset_count()
print(n .. " datasets loaded.")

total_max_y = 0
for i=0, n-1 do
    max_y = F:calculate_expr("max(y)", i)
    if max_y > total_max_y then
        total_max_y = max_y
    end
end
print("The largest y: " .. total_max_y)

If a fityk command executed from Lua script fails, the whole script is stopped, unless you catch the error:

-- wrap F:execute() in pcall to handle possible errors
status, err = pcall(function() F:execute("fit") end)
if status == false then
    print("Error: " .. err)
end

The Lua interpreter was added in ver. 1.0.3. If you have any questions about it, feel free to ask.

Older, but still supported way to automate Fityk is to prepare a stand-alone program that writes a valid fityk script to the standard output. To run such program and execute the output use command:

exec ! program [args...]

Fityk library API

Fityk comes with embedded Lua interpreter and this language is used in this section. The API for other supported languages is similar. Lua communicates with Fityk using object F of type Fityk. To call the methods listed below use F:method(), for example F:get_dof() (not Fityk.get_dof()).

Note

Other supported languages include C++, C, Python, Perl, Ruby and Java. Except for C, all APIs are similar.

Unlike in built-in Lua, in other cases it is necessary to create an instance of the Fityk class first. Then you use this object in the same way as F is used below.

The fityk.h header file is the best reference. Additionally, C++ and Python have access to functions from the ui_api.h header. These functions are used in command line versions of fityk (cfityk or its equivalent – samples/cfityk.py).

Examples of scripts in all the listed languages and in the samples directory.

Here is the most general function:

Fityk.execute(cmd)

Executes a fityk command. Example: F:execute("fit").

The % operator for the string type is pre-set to support Python-like formatting:

= "%d pigs" % 3
= "%d %s" % {3, "pigs"}

Input / output

Fityk.input(prompt)

Query user. In the CLI user is asked for input in the command line, and in the GUI in a pop-up box. As a special case, if the prompt contains string “[y/n]” the GUI shows Yes/No buttons instead of text entry.

Example: TODO

Fityk.out(s)

Print string in the output area. The print() function in built-in Lua is redefined to do the same.

Settings

Fityk.set_option_as_string(opt, val)

Set option opt to value val. Equivalent of fityk command set opt=val.

Fityk.set_option_as_number(opt, val)

Set option opt to numeric value val.

Fityk.get_option_as_string(opt)

Returns value of opt (string).

Fityk.get_option_as_number(opt)

Returns value of opt (real number).

Data

Fityk.load(spec[, d])

Load data to @*d* slot. The first argument is either a string with path or LoadSpec struct that apart from the path has also the following optional members: x_col, y_col, sig_col, blocks, format, options. The meaning of these parameters is the same as described in Loading Data.

Fityk.load_data(d, xx, yy, sigmas[, title])

Load data to @*d* slot. xx and yy must be numeric arrays of the same size, sigmas must either be empty or have the same size. title is an optional data title (string).

Fityk.add_point(x, y, sigma[, d])

Add one data point ((x, y) with std. dev. set to sigma) to an existing dataset d. If d is not specified, the default dataset is used.

Example: F:add_point(30, 7.5, 1).

Fityk.get_dataset_count()

Returns number of datasets (n >= 1).

Fityk.get_default_dataset()

Returns default dataset. Default dataset is set by the “use @n” command.

Fityk.get_data([d])

Returns points for dataset d.

  • in C++ – returns vector<Point>
  • in Lua – userdata with array-like methods, indexed from 0.

Each point has 4 properties: x, y, sigma (real numbers) and is_active (bool).

Example:

points = F:get_data()
for i = 0, #points-1 do
    p = points[i]
    if p.is_active then
        print(i, p.x, p.y, p.sigma)
    end
end

1       4.24    1.06    1
2       6.73    1.39    1
3       8.8     1.61    1
...

General Info

Fityk.get_info(s[, d])

Returns output of the fityk info command as a string. If d is not specified, the default dataset is used (the dataset is relevant for few arguments of the info command).

Example: F:get_info("history") – returns a multiline string containing all fityk commands issued in this session.

Fityk.calculate_expr(s[, d])

Returns output of the fityk print command as a number. If d is not specified, the default dataset is used.

Example: F:calculate_expr("argmax(y)", 0).

Fityk.get_view_boundary(side)

Get coordinates of the plotted rectangle, which is set by the plot command. Return numeric value corresponding to given side, which should be a letter L(eft), R(ight), T(op) or B(ottom).

Model info

Fityk.get_parameter_count()

Returns number of simple-variables (parameters that can be fitted)

Fityk.all_parameters()

Returns array of simple-variables.

  • in C++ – vector<double>
  • in Lua – userdata with array-like methods, indexed from 0.
Fityk.all_variables()

Returns array of all defined variables.

  • in C++ – vector<Var*>
  • in Lua – userdata with array-like methods, indexed from 0.

Example:

variables = F:all_variables()
for i = 0, #variables-1 do
    v = variables[i]
    print(i, v.name, v:value(), v.domain.lo, v.domain.hi,
          v:gpos(), v:is_simple())
end

Var.is_simple() returns true for simple-variables.

Var.gpos() returns position of the variable in the global array of parameters (Fityk.all_parameters()), or -1 for compound variables.

Fityk.get_variable(name)

Returns variable $name.

Fityk.all_functions()

Returns array of functions.

  • in C++ – vector<Func*>
  • in Lua – userdata with array-like methods, indexed from 0.

Example:

f = F:all_functions()[0] -- first functions
print(f.name, f:get_template_name())          -- _1        Gaussian
print(f:get_param(0), f:get_param(1))         -- height  center
print("$" .. f:var_name("height"))            -- $_4
print("center:", f:get_param_value("center")) -- center: 24.72235945525
print("f(25)=", f:value_at(25))               -- f(25)=  4386.95533969
Fityk.get_function(name)

Return the function with given name, or NULL if there is no such function.

Example:

f = F:get_function("_1")
print("f(25)=", f:value_at(25))               -- f(25)=  4386.95533969
Fityk.get_components(d[, fz])

Returns %functions used in dataset d. If fz is Z, returns zero-shift functions.

Example:

func = F:get_components(1)[3] -- get 4th (index 3) function in @1
print(func)                   -- <Func %_6>
vname = func:var_name("hwhm")
print(vname)                  -- _21
v = get_variable(vname)
print(v, v:value())           -- <Var $_21>       0.1406587
Fityk.get_model_value(x[, d])

Returns the value of the model for dataset @d at x.

Fit statistics

Fityk.get_wssr([d])

Returns WSSR (weighted sum of squared residuals).

Fityk.get_ssr([d])

Returns SSR (sum of squared residuals).

Fityk.get_rsquared([d])

Returns R-squared.

Fityk.get_dof([d])

Returns the number of degrees of freedom (#points - #parameters).

Fityk.get_covariance_matrix([d])

Returns covariance matrix.

Examples

List peak center of series of data:

-- file list-max.lua
prev_x = nil
for n = 0, F:get_dataset_count()-1 do
    local path = F:get_info("filename", n)
    local filename = string.match(path, "[^/\\]+$") or ""
    -- local x = F:calculate_expr("argmax(y)", n)
    local x = F:calculate_expr("F[0].center", n)
    s = string.format("%s: max at x=%.4f", filename, x)
    if prev_x ~= nil then
        s = s .. string.format("  (%+.4f)", x-prev_x)
    end
    prev_x = x
    print(s)
end
=-> exec list-max.lua
frame-000.dat: max at x=-0.0197
frame-001.dat: max at x=-0.0209  (-0.0012)
frame-002.dat: max at x=-0.0216  (-0.0007)
frame-003.dat: max at x=-0.0224  (-0.0008)

Write to file values of the model F(x) at chosen x’s (in this example x = 0, 1.5, 3, ... 150):

-- file tabular-f.lua
file = io.open("output.dat", "w")
for x = 0, 150, 1.5 do
    file:write(string.format("%g\t%g\n", x, F:get_model_value(x)))
end
file:close()
=-> exec tabular-f.lua
=-> !head -5 output.dat
0       12.1761
1.5     12.3004
3       10.9096
4.5     9.12635
6       8.27044

All the Rest

Settings

The syntax is simple:

  • set option = value changes the option,
  • info set option shows the current value,
  • info set lists all available options.

In the GUI

the options can be set in a dialog (Session ‣ Settings).

The GUI configuration (colors, fonts, etc.) is changed in a different way (GUI ‣ ...) and is not covered here.

It is possible to change the value of the option temporarily:

with option1=value1 [,option2=value2]  command args...

For example:

info set fitting_method  # show the current fitting method
set fitting_method = nelder_mead_simplex # change the method
# change the method only for this one fit command
with fitting_method = levenberg_marquardt fit
# and now the default method is back Nelder-Mead

# multiple comma-separated options can be given
with fitting_method=levenberg_marquardt, verbosity=quiet fit

The list of available options:

autoplot
See autoplot.
cwd
Current working directory or empty string if it was not set explicitely. Affects relative paths.
default_sigma
Default y standard deviation. See Standard Deviation (or Weight). Possible values: sqrt max(y1/2, 1) and one (1).
domain_percent
See the section about variables.
epsilon
The ε value used to test floating-point numbers a and b for equality (it is well known that due to rounding errors the equality test for two numbers should have some tolerance, and the tolerance should be tailored to the application): |a−b| < ε. Default value: 10-12. You may need to decrease it when working with very small numbers.
fit_replot
Refresh the plot when fitting (0/1).
fitting_method
See Fitting Related Commands.
function_cutoff
See description in the chapter about model.
height_correction
See Guessing Initial Parameters.
lm_*
Setting to tune the Levenberg-Marquardt fitting method.
logfile
String. File where the commands are logged. Empty – no logging.
log_output
When logfile is set, log output together with input (0/1).
max_fitting_time
Stop fitting when this number of seconds of processor time is exceeded. See Fitting Related Commands.
max_wssr_evaluations
See Fitting Related Commands.
nm_*
Setting to tune the Nelder-Mead downhill simplex fitting method.
numeric_format
Format of numbers printed by the info command. It takes as a value a format string, the same as sprintf() in the C language. For example set numeric_format='%.3f' changes the precision of numbers to 3 digits after the decimal point. Default value: %g.
on_error
Action performed on error. If the option is set to stop (default) and the error happens in script, the script is stopped. Other possible values are nothing (do nothing) and exit (finish program – ensures that no error can be overlooked).
pseudo_random_seed
Some fitting methods and functions, such as randnormal in data expressions use a pseudo-random number generator. In some situations one may want to have repeatable and predictable results of the fitting, e.g. to make a presentation. Seed for a new sequence of pseudo-random numbers can be set using the option pseudo_random_seed. If it is set to 0, the seed is based on the current time and a sequence of pseudo-random numbers is different each time.
refresh_period
During time-consuming computations (like fitting) user interface can remain not changed for this time (in seconds). This option was introduced, because on one hand frequent refreshing of the program’s window notably slows down fitting, and on the other hand irresponsive program is a frustrating experience.
verbosity
Possible values: -1 (silent), 0 (normal), 1 (verbose), 2 (very verbose).
width_correction
See Guessing Initial Parameters.

Data View

The command plot controls the region of the graph that is displayed:

plot [[xrange] yrange] [@n, ...]

xrange and yrange has syntax [min:max]. If the boundaries are skipped, they are automatically determined using the given datasets.

In the GUI

there is hardly ever a need to use this command directly.

The CLI version on Unix systems visualizes the data using the gnuplot program, which has similar syntax for the plot range.

Examples:

plot [20.4:50] [10:20] # show x from 20.4 to 50 and y from 10 to 20
plot [20.4:] # x from 20.4 to the end,
             # y range will be adjusted to encompass all data
plot         # all data will be shown

The values of the options autoplot and fit_replot change the automatic plotting behaviour. By default, the plot is refreshed automatically after changing the data or the model (autoplot=1). It is also possible to replot the model when fitting, to show the progress (see the options fit_replot and refresh_period).

Redirecting the plot command to a file saves a plot as an image:

plot [20.4:50] [10:20] > myplot.png

For now, it works only in fityk (not cfityk) and is less flexible than Session ‣ Save as Image.

Information Display

First, there is an option verbosity which sets the amount of messages displayed when executing commands.

There are three commands that print explicitely requested information:

  • info – used to show preformatted information
  • print – mainly used to output numbers (expression values)
  • debug – used for testing the program itself

The output of info and print can be redirected to a file:

info args > filename    # truncate the file
info args >> filename   # append to the file
info args > 'filename'  # the filename can (and sometimes must) be in quotes

The redirection can create a file, so there is also a command to delete it:

delete file filename

info

The following info arguments are recognized:

  • TypeName – definition
  • $variable_name – formula and value
  • %function_name – formula
  • F – the list of functions in F
  • Z – the list of functions in Z
  • compiler – options used when compiling the program
  • confidence level @n – confidence limits for given confidence level
  • cov @n – covariance matrix
  • data – number of points, data filename and title
  • dataset_count – number of datasets
  • errors @n – estimated uncertainties of parameters
  • filename – dataset filename
  • fit – goodness of fit
  • fit_history – info about recorded parameter sets
  • formula – full formula of the model
  • functions – the list of %functions
  • gnuplot_formula – full formula of the model, gnuplot style
  • guess – peak-detection and linear regression info
  • guess [from:to] – the same, but in the given range
  • history – the list of all the command issued in this session
  • history [m:n] – selected commands from the history
  • history_summary – the summary of command history
  • models – script that re-constructs all variables, functions and models
  • peaks – formatted list of parameters of functions in F.
  • peaks_err – the same as peaks + uncertainties
  • prop %function_name – parameters of the function
  • refs $variable_name – references to the variable
  • set – the list of settings
  • set option – the current value of the option
  • simplified_formula – simplified formula
  • simplified_gnuplot_formula – simplified formula, gnuplot style
  • state – generates a script that can reproduce the current state of the program. The scripts embeds all datasets.
  • title – dataset title
  • types – the list of function types
  • variables – the list of variables
  • version – version number
  • view – boundaries of the visualized rectangle

Both info state and info history can be used to restore the current session.

In the GUI

Session ‣ Save State and Session ‣ Save History.

print

The print command is followed by a comma-separated list of expressions and/or strings:

=-> p pi, pi^2, pi^3
3.14159 9.8696 31.0063
=-> with numeric_format='%.15f' print pi
3.141592653589793
=-> p '2+3 =', 2+3
2+3 = 5

The other valid arguments are filename and title. They are useful for listing the same values for multiple datasets, e.g.:

=-> @*: print filename, F[0].area, F[0].area.error

print can also print a list where each line corresponds to one data point, as described in the section Exporting Data.

As an exception, print expression > filename does not work if the filename is not enclosed in single quotes. That is because the parser interprets > as a part of the expression. Just use quotes (print 2+3 > 'tmp.dat').

debug

Only a few debug sub-commands are documented here:

  • der mathematic-function – shows derivatives:

    =-> debug der sin(a) + 3*exp(b/a)
    f(a, b) = sin(a)+3*exp(b/a)
    df / d a = cos(a)-3*exp(b/a)*b/a^2
    df / d b = 3*exp(b/a)/a
    
  • df x – compares the symbolic and numerical derivatives of F in x.

  • lex command – the list of tokens from the Fityk lexer

  • parse command – show the command as stored after parsing

  • expr expression – VM code from the expression

  • rd – derivatives for all variables

  • %function – bytecode, if available

  • $variable – derivatives

Other Commands

  • reset – reset the session
  • sleep sec – makes the program wait sec seconds.
  • quit – works as expected; if it is found in a script it quits the program, not only the script.
  • ! – commands that start with ! are passed (without the !) to the system() call (i.e. to the operating system).

Starting fityk and cfityk

On startup, the program runs a script from the $HOME/.fityk/init file (on MS Windows XP: C:\Documents and Settings\USERNAME\.fityk\init). Following this, the program executes command passed with the --cmd option, if given, and processes command line arguments:

  • if the argument starts with =->, the string following =-> is regarded as a command and executed (otherwise, it is regarded as a filename),
  • if the filename has extension ”.fit” or the file begins with a “# Fityk” string, it is assumed to be a script and is executed,
  • otherwise, it is assumed to be a data file; columns and data blocks can be specified in the normal way, see Loading Data.

There are also other parameters to the CLI and GUI versions of the program. Option “-h” (“/h” on MS Windows) gives the full listing:

wojdyr@ubu:~/fityk/src$ ./fityk -h
Usage: fityk \[-h] \[-V] \[-c <str>] \[-I] \[-r] \[script or data file...]
-h, --help            show this help message
-V, --version         output version information and exit
-c, --cmd=<str>       script passed in as string
-g, --config=<str>    choose GUI configuration
-I, --no-init         don't process $HOME/.fityk/init file
-r, --reorder         reorder data (50.xy before 100.xy)


wojdyr@ubu:~/foo$ cfityk -h
Usage: cfityk \[-h] \[-V] \[-c <str>] \[script or data file...]
-h, --help            show this help message
-V, --version         output version information and exit
-c, --cmd=<str>       script passed in as string
-I, --no-init         don't process $HOME/.fityk/init file
-n, --no-plot         disable plotting (gnuplot)
-q, --quit            don't enter interactive shell

The example of non-interactive using CLI version on Linux:

wojdyr@ubu:~/foo$ ls *.rdf
dat_a.rdf  dat_r.rdf  out.rdf
wojdyr@ubu:~/foo$ cfityk -q -I "=-> set verbosity=-1, autoplot=0" \
> *.rdf "=-> @*: print min(x if y > 0)"
in @0 dat_a: 1.8875
in @1 dat_r: 1.5105
in @2 out: 1.8305
PKHGR^+fityk-v1.3.0/_images/fityk-with-tooltip.pngPNG  IHDRGNbsRGBbKGD pHYs+ IDATxwx?WriAh қJ""( "( "QQl4KS4;^]./ܓܖwywv{3(ʨQ@ ,1gAӰaC @ @  9" (1ԮJԺOW}aÄ@ ̛7O8A pEGnهG5o˽^z O @ iRJB#)ͷTJ᫛*;O[n:XY@ ۷o[|yNG () ,FʡS0K*Uz9∪BFp=p=lwF&= fdUxU7̮VI\g'i<7i.=/^6U@ U?Ht2)R'pU)2{P?y:<5=1]Xw.p#ROzE_!:! W } O?m321\t?1'bKpQOI!e|D$a3[el_c ij$c nO]Vஏ˯u50i[Y!vV̟陸PPHwkD ObN`70UuE]@ (a$IW#Iqb୲vK([w&HX;6,5qh*vCScr")f ^|X灌$˸F׎*GYx4G6i *$1Ԥyfϳl$@d,eOUj$ c>%?4㥘|21_n&ίNeHC-J1ܬ@myz/Mc-j?cQCIi n|[!b%@ Q(8bXv*xxz?d ;RT~ tCZ-^r8(/&3gҳ< k|;5~#_|,> I6KFV/kY0ax7IDuF|yg 8~`PRJ2Smʻd/aL WU/71gwb"R2=+OQK;{'>b-SAUW%N4"OutsXWaOش8];RCoJ mPK}Q:5Ƹ3y?j[.8cdzN>ZȚ ۩~OF|cef$$YEsUYFvLHO2ny$&@P(0]۞|cKv^AGQ?ydd˹W0%ƯR]ty '7Rc Hr\'nn6TBJ%i9!mҩ1޿#ANn9l:G1~=3vt|;Iߝ"oBLmMs Y C#ϻs%h4CZlgI>M>uyAl M<<!]ΌfIP65NUz::?pM[Fݭ=b]RMxѴ[hB(N#ݒv* 9Ǟl7ܼ 5Qt*ܼ=Jv* 9GdrUYF-igu *@Pf(O=.?"F q{)8'eYFuX: ]^y#)] [wN8D @T C]svĵkh7gcp[`m]mh}E5rzg{׿?lˆ{g0nqrL{]sңht4 `onI'bBwR̛I+/ wzsN9q>ZUsJ|@gh ӧ<į}Qs/SFGh-=iEk//|8w97ܩ /mFFޚp.;Bw*n&Ԍ0(W ZΧi]iĭgO8n-a vv#;Z磬3 뢑4dI-Gs&RjSm=jt[zc +&٨7&ǶH*~1zR}1eO`zz3>PqI.:/89s!tT~ BU0`=FޕT)BBEg5TjUS R|SoTZ5e}׺}fw|Z"c`tQl8}5%o:=Сj̪4cDz XUMvzGr*4 }MS+IiӼ u@>$j^)I3?CDLdNUdd| GMTiM)Kֲ7YyCƺq}PӎatlȈ#IǓ,wٔO?@L13YALFڼhIv EV&*!I ) (!U'}Qjwv;4m[ F Feh& iFwkώ5 LL>3_$?Lˇ;! 5Aw+Bjl" %V/3d YoYY}U[ou)Kv[}K<)\&bPӆh! 27p\g'$ %amnT=$ r;8IOXuWʳuEFVeu5o I7.dd$T9[]~Q$&*nT ITxf#! e.*hBB=#VUyAK l}ƒ%K7hР}ҥKs<,]4O{2!8rQc!qFvZm!rAIr`Ȳb::@gF3|,\?7M(w~ UZv$q| vX2vlӜƪ.gU4+þ7AơSdڱ&I'6/&{1 ID:~?A)[YYvI'67&-wEfP5Mi!i ' | <) MKJ6wl娡6 ǫ#$פ[uff)eGٹK+^fŻP&<3׸bED٦u(Hmޕ"u,c\G#ICt@ 7PgBJ`mOźN-xv'֢wGޚ]ݧWƼz_sL_to֜3uk07h~'&a]KhAǯu gxB GV^B#1 \ #+9Ɨ@ AY"0 3Qk#Vn r[xtaXRJ:̢_ӯ뷢5r kOjJ#ﶎt3$>tlٝ%-vJ2mt BɓbP+žKih|W>U ܋wv $Po@v糬swiK*h%)hUWwkNLq\ _Fn ?6Nߦyvt]nYCT6y}cAE{GIsOՏRґ\$9$@q{zr̷xX:7G `;H6gK_eW} zn@ ;XlR_햴ls@Z@lo,abH!ݢ-1X;ƗJ @Pv>RgqM'k䔎m{AC(FaÆ{˵jv ưaXH!j/|u214,ξ.bL!!1@B^#!y7YQ+}7ӟ#/zmH5~&RI6*HO|?#r߸ư٧pB-9JG\oIIJ!hWiyzU!)$`a&56d3hт))4 ?/4SH33^߫B ;Tv;mIK%5d3u9ׂ=O .s۬]c #ktIR d72|PR҈Ő;xŜJbb*F%2j;^^1H"(i~xƂyh$u,1 *L-*9cqFqx#,L+y[%I7]Cev#Iɤ wO/5`L!>nSvS"-O3&3C%@PZ̛7y\9!yp2Α#G3g(lAtȍ@VAu<湎s}_/u)z~n]xPU[oG- d8OF_)ñ54F7wk0dhuWo o?w[UB_4d<dS/Ov -M/fvfؚmސʳr-ڭ;a1*//rs SV9E]5_IoXR-Aj7ޞNc){<:;Nc*tQm;nxxϺT 6O!Xo!@P(YK3đ"L̂3i*[ t}w\ R_vK@H/eI\r @ r^[&HdPbEz7T*Hʸ{ib,;㷿Njr'm~^ra>/\ǙsG[8rOź@ W#Nyr;!uY@ +oE9"(I9Ta@ @ '#G@ @.5M @ @  pW@ @F#@ @ ا("!@ ?9"@ ?G@ @ )jL|f6Fٽɽ^>@ @#G 1m@{jDHt2c,h5=J:ޥ;EeUokmM:`"eX_^X 6Pw,{SC$e l trM1jD!fЛ5BvWKGF3gER (SSY-LIA06XtN2<_?p˂@ @ 9~\M/!Z埿NbFBɏ?HJm| =~%dbcOsylxo=K6+lnK *`!gy,+XXp߅8$ZR94 _{es H:kh 7+.%..'Op5sjF~}o:u[~ML!%-3;+9dC)>@ K7Ά>oTԘ:Ѽ1+fh,iY㴨L`P0:ݱN L`: y3]㛮Uh/{+ F~Gd^9>nS3eSmLpnKFj׬Ii@^_C}2uFYtN/C-kZR^Z$HAxSSRP0_fWJnf07lHDƧa<0aYy`Vר)vUSpkyWxqe<&s~XX׷]>~$P6nkP9^߹]b Ƽc,^ݎ& Ѩ<1}~ḙNF\94vXeBb?L @ 8҇Ę`cľ'Wƾ2>~~at;F&`wD|U#ݵovKfs>?}i%ǯ\d ~{uGn&˦}x{xڷ?>LO/:;ͣk˱@1rG򡣎zon懃*ھܟZnjU^9z+NX^8F "@ ܣZ_}'Ou@ͦh~?ZaY+u _;~K*zVϢͷU*1*؈A* R-)&6EQ0Fe_j~:Zc3GCyZG T{h(Ά3))Y8O)vp4 Uڄ{8w1RGR߶Ly{kƭ[h6]( $<(sslqlSD2,ۖ{V]ìZc"Tj̓yT\[ݽr{Y94QsyIV\'A zاH^@yZ9LCROs>sEYilLQ[>8-w&b+?dg|DK/:ʸ@ @p/#y#|n8dc,煴]KNUjbQ=8.bL,v?;QS­Fc~Rzc۫ jN:1hu\ IDATZ<)KJJ1h_WCD3 SxtP{ěrx uwfkhxWAͱ|ӯR!>fAQѭqc5*3/*yI r]iQkx񉅔޺ ӟdBl+[&TLSK~qD @ )wWe=/N6ᥐ#|6{'IjœC׸yzqѺel]:nc~<[,۾:=GP\)Or>ٍUc#6;O_n)9*3kؐ Ұa}+h -'cvcʱ1SDb?=R,i-^HG $N VH*4jC!kX $@QeFyi|TFs.@ =+(V⒐*6:'Y+sNǫЏ~&z5M{S(hBgD{{Cpo$'g*J!YY;o'umx4õ=74乡!c$ILs%d %rˇz}&v_#3Siw Ԯ澣1e$fezAhS.)Lh+C0qy1b 7ػ?zU] #>#>#>#>@P*H^°~/3 >vv|hTѬw11o63Vµ?%'茁xwJv[ؿnM/;ϯ|b=>'lCbw9Μ\-3>@o8mzTQԥ|?ĈBٟщ̟zBCD37KF@J|8WN(()Y13+*A^'r)2ǶM70:ɣX)GaL]>H33w<_'Nrt~g,TH'S*v<˩l5eC 2err¸kbY6:Y`}D@ Q}өZd]M4?2u`42to0^|Ѱ\_!A=54sw4U'>Íжc}ik<1'^!mu?o4;xn}ǯRk~AvͶM]fZ }Mƃc^48燏GUC/9 ?:}^cs7ݟүtj@7z9jE/[KSMʲ6!@ ȗ8⬓m2" 0F$HRLoJj*8A1 4_,(6D @ %I8mZf\OT ))hQ2,#ɲ 9jD߾me 1-H`8qb@3Q(( H A F,pYQk4 gb#s͑lH uJ$: 9rgċN7%CQH#A()Dq!@ HMN&22G#/Gez=uգJ{z {(dyp*zߝ)hH"өSg1LY]] Ǖ$8QWcHOZ5K.]R*J8_dl5GG8>B]= S{7t\DˆȑhgbA8m\Kr,v(T|z6Z%xq℻ELnݼA*9c_xb"I]w&-έGSڴm,t҅ R\9, ofϞ=4h]zÇЪu;RSKѷ@p-\"%SϗlJ]be8' n>9o.M+|o!Bfhѣtҕa&/tVjk I;+j UO1qj)Y9S@Pqȹ(O<5 ={c2ذnmb8֑ʷ%(TVK)y7hڼU|g&bhgftMt18Iՙldℋ;9![кuky_K@ر#˖-j\x6dF:} wK(bp$G"aĜأ w;h\+WNӨԶ1&<驅Mqeɓ'|p>`ٳg-[,r !ؘnøǍ.fx53ZkڨY]Ù8V&Z1*Wk~h4;zo))7׹YRZ˗X 6ٳțSf6{֯[ˬ!1!"6r~%8RޥNh5h%s l}|7uĶ- Q(slq߆Y&Ehp̘s\@!qJ0aP2OF4o2߱VZu.^8_[a4j܄Gy{ wwwt:^^^t:j5$=z4T\+/UW͛Ԕ4=FxR&LXd[6Yl$Zn윏)51ONAY,_hU~6Ӓ3&ITp?KyQ$'g2JbbbļǗO&q=%}?eÝU*KJ *nLM!G1jqq?Ųt1t>CǢУW_0HrbX*~.ܿ[d\S@mJD܈X*ϵ^O _~Jtt4Ϝv/go(=osԬμK%v.%i@qā^&f9#|~+H8v[~vgۖب4)Suzk*u%[´iӘ6mZ+:GQ#G$?Ɍ&N(\gJUw&ҹ9{{ohԘ֭[ӯ_?<==EQX,j<==aСDFFZ W\.w<|HĂڗ8g>a"7z_6w9./ >LqŋyYp! `ŊՋ6mJf-ؿwwB~4٥Kآ^\hH*%|m_ VǿDݦ;Ic`#}|wQE>qWҟ GIh4h4dYDQgF#&Lw&39p =Z${:wl߷dwپbO޽ۗիW=Aၯ/Hd_tq0xzzRzuγ>z im޼9cF#Νc_`ۃDt SW@Z*VkNܹHեXWƱc/,&ix;ؾ.^^QD (;:.JU8һ-DqLfqcmO|'zAz~S!CgСȲ̐!CP!C`64h#GsNݛ%|F\j3*Xvo-:8/sqmͨ@Jb@} ӮزX-q2y2D3 6䥘et&Z}./IEQPTY:·c08sm;;wرc pɔcg8c˲ly_>۷/Q{*T[ouqfb̘1ꟺӷo_^6@ƣM0QFу`n޸R)) /@pp0>>ֺ:ۥVÔm[rr2~~~YGW^t:~W#GJhZ>Ύkj\zde_ͩy,W7"X47gaYlԙ-g]y,C;EwK&Em[EuoȲLV/k?Gpa6n܈?y}n:[[(nh04KM#~gFk=c]玁<=sbK|}.k)R9E鱄 k#@vܙQha(E^"yPgFR;ͦ |Ҭ<}}cŠML-'￧q.mg5o7Dt;/ բ?@Qr_._LB|,BrJj'L[VTq$f[:} eT҅ T[t6$===ET!աS `o,C94[jMjj*GpIfz:ȑ0aEV\eϣ x]&$$nݺaÆb͋mm$s~N:ufo.UWϋf@h籘 ~7&ܬ yc>I;Km-lDžvR2StCeN$uW^eMM uZFefӿ~ 'y"n?#&{a=iuЗMGypdGx $cO77Ij[<|M{$6<<@LL <ׯϿ@pruJpE[DIrօ] )4~7g`}IBٲu6+y1aG/Gq_ҚOXGUVX%OAYCv)8S<`<ÊU nk{{{ۇgo<7>?1?[qM\ ,%n8^0纚g}(ݾ}{u/CW$:l6Sr,欣,fcU8vns8=C a۶mtСbq'Cl۶=7Æ ~xjذ!FرcLhh(Z._-l[6ȱ-mbѣ>,7)&^$ Gw3/>z ;v`0pg9|0?C#4hPw?ocšT}^5{nQ_cº)׋z1al:7\w:v߿YI9҂}n!1xUz/9>KVn(SP(vniJ:u׸ϼf3*ɜfUqѹzvZM~|){x'Ɠč<օ-wfzxը[wIs9QV{SV]12v>Zn;del?ZUcVԓh{=k<ۮxX:]F0J GL<&OL=8x`ҍr4pjϢط5pƽGpU=ekkbMP #DkW3[[sty~ 7vΧB9/V}Ae,oI( iiUzܘL&BCCqssq<{N M9t`_i;„;;WG;36V3Ͷ u͛7~jP5mƶG\lL?p- zKF}2VwkK,h4%\;[gBDw=vbٲem۶mc۶m̝;77n܈jO6g3b/s_}$ѳgOT*{&&&Xc렠 R/]$t:fVkicX0FL&. N%ݩ 󗾙΢(p5IHHuhZ_q1{EYY>hݩ;UnlAZǚp2WaQN-Ӵys>fbyp$[.MS&'vm*KQE'=k7bC`ў_lŷb[-B2칺Ow|Eu%R8(8/7u\>מWoٰaXu'oufZu\kф>Fd (/-y7@rS 4#0hԙ*;Ǿ} hQP1l| Cbоb:G 8;3sgHo7V͇C>QK?>ݍWM}#H-N%¯/\~ek[XK}+?>[31_ہ}~}0ßg}Lś+as봚gLQ\>KYYM43fX`Ĉxzz2i$}PذMmu#-]ϒ3,y5]l" {`oo;}|r ׏:qw넏3gLeķ6_4:712ɾ:${*4%(6n':7]q` Ijck;a(bwdždYѾyp{`/xjԻ4=}}{)/u)'xRq IDATC*fs'uN#J?BRq ٓ9ٳg '=%ԫ֑|-[$,,Ç?M||@W j9||ҏ[ÁY3(=p'66Eӈ 8aaEJ4 fقRTT:}gLa$((5:#YBK)?EQqLJ,zWi5 u%b|s*P^{}'U`cY j["5 WW"Қ<34FO0sn2UCC0m-"?5HzVr/LSz^gg5ny3׸*U#Ϟ! ,wzj0bWbJčȋ\4f> z,s-m}%~[Vy`^Xs{CEblRޭS:f'~~q~?s#Ƃ!իoA#{tNS%/^bŊ,GuX,L&f~&85Ւŋ\&L余3UɎMq# ˒WSf--6jլG0nСC'V_?bSm6r$66_Ȑ!C$e˖9AaX㏋3i6ľNCQL&S{ⲧy3!Cжm[ׯC=;IIIx/$I"!!JoŞuq1}]mV w 5 E?&=}Cҵ_QyT!-%e /Ax@$L&_;>n6=! zh ,(XX(ĂQ| H @ }f7Kz|>̹̽3ss ]okGSw+?=E1tP:uDf8HNN۷o6ΜΠglX;SdE7Ɯ3#ºя'~:߶3o<O*۳眈#z̬G\"^Ǿ GƴlلfЄ NQտ^qi<t_TԞ+ $Uy868 u-{͵Ɓڨێ\|sKxb1P:P AnLׇ35:u UUX,tA\ەOW_ꏃ…\1thoٲm(rn ܫ Dk"ݝ>֭⚫rF4A ٝ]x+ǖNc-l氮-Ou(Vn}vm?hWy-g $=9_~c/0E`YJ^l.)\[לk<֔.*8{kZ&X~=_8-^QO9syhy?{u˖D|zϯd-qDUU|^}' Vs S;6hvh"."Ơ^ l4ZV1YZ[FfY{'+1sۗ_|>*\necc.bE=;aw' {P~[GGpNxķ5C0ݎNby9{jxkUCM<}vQEaϗ_~g}ƚ5kh׮f|[.x_YK.[o)89&jĻ fƨaa􋻎#bpiڵs2DS|>=Ctj`#uyZ$11Fd>яݻ4cJIގ~~~~ǥXl"őZ4NIM=<Eh-sv6]C)h]ո>jCT~QQ{>u~\%؟\2&0*hmb"C6qgrOwmߵk0fE T*lGè@^Z]v%۔/n$=]b1` h4( u&)c>ӓ*vlj:?7nH󉊊w^FTUu }2d;t&a{a*U%u+FTC,ً@} 4%:P >skp%^ U&/":XoFT`ٌMT;Nb3zjƬa0 Ujb6ӵdmy/lC@3fkim4l,V;*V=+*WW%eZWgҞmb}}MGYhaq8jxIrJNHѣܙ֭[be.q*ыݞ!Cؤg4H밡 Pm{#4krssXjCFnnN ݦzf~gFIVZFqX,7_CKa*9;C(f#vO*)KvO:[Tkʺ ֡|>cn:ΜɤADd^gQ>}ycPݎi=Yf( wu!ǹVU6y\~}1Kq|-ÇҤi3=ƍ#..%Kп|}a-Vez}mpgɒ%qt޽Fqvw 58d>)ұ?0rHP3g^#Gp/5%ZYޠAw5pQ+xŨJ\\z$O?DvGO\f3)BČ(Oq{GϜ9FkV~۟ysҨ;-rwƔӡ TЬr`k,6|`ʔ)7aw#i'VUa[C`w|ƛ)"Su=S0fРAm6,V ˖-EߍԔV+5!9:t PjCi;vqMG:5IK~5rC '#e%%W%L0t%V[˦&踟la_F}MظM,/.uXfM]zt -~Hn|{(ጘ~1A<w!־U{"o&MCt: L+nn~˯~Nmϧ#ˆyCVbٿǏ $4W^ˡ16f77?V-s-7vzOT09={U{@ZBGZW^ba/V(ukVo@^{5}]F#&'9 377W6l+/T%\Gy!, [YJR֝ ɬO5k,iǏW@RQUU]r 999s)Fv1i*e\{F#ݻwW^zSo:q63qDV+9995jyG]ñNW͝;XaIII@a?NƍٷoM4q8y$pԩAn5iz U[9;{gڰYm.j͚5׏{7Wڞ-,#:~Gl/FuAԊיUl\/]x dMt:FR EgXHZ|t&Ϧ9.xX|zMzzk>b @|Ba1-5膱(W@V7kK#2-l{ 7ʀL*=0KĪ?2{6a@Pޙ]O/vM!^7M%0cl0|p>Ād2;Z?ħ]dh:2i«!>wl+KzFWsyCmӉpÇ`4ғb5;o aήS2bu0zhfxS;l U{eF6m\=z{4iC.GҥKgKn^ciE˫\IغzhdM riԎ_Cj8~X ǵ_xޱ'N~w?Y~+0qD|||0 7oFs ^^3gdn5E}l>~ĴɏĽو\ZK*uv]{o?pE{~w짖E^*F~.SaO?q(6>4-ʴixq =ǹپ};'WMg6+aOE8z$eOyUdOqdeeU=Ӓo$"1V{,0cƍݻu-@yϿG_<wp}*Cdd$wt _C !::8Vh]i[˛揤9哸5 `6ݱXm-ΜdΫ+k@ǎԑ핶(Gls=1e=RbޡrJ٣JFV: @$buBծmhhXMfS C!'pH텍Ȑz!dfuG@(>- _TiM-;Ci2Y5Q~í h2ͅ|$&^ۙC )! O>qpYWlψG'p3)qp\Y|ÚoC[a(t'Cs? }1N*tEXX>NJykREr8t aż+d|Ddϡ tz13N1o0RU*N{-\P=O=]]Lm0A~B;)Ѳu7oD`/GPV7s*4McՊdhNX,޽Yfc~*W>@.T(,RQ4iK it>ocDe8y4h]nnM@ĞLئ#D7lTsxn̹C"̞".1L6v>}zXzK$)NnH0ZU5P}> cS w88>EwU{`gg*WCip)݉ +ޗ'~J\ }Vq߬vو5S\?()>Vr[~%m(' gJ ;3goy]|t/^k|_}K{/:E]<[ң+w˗C]W hG>߫e>6oT(>+wX7FuM44*(o6 ᎑r%T{vn/xQ<3D,}Ȥ8iQZ4ӏqqK<ӪԯbSa^nDT4W<9%oel޼6:3y*]c׼V5jDvv6l?.K{4hQ# b1a6[ NaW^nwAQ-Uɂb%ps~[ Ke:nW@իWӯ_?ϟ'~[eN3FuΝ;H%mQqCFA)dga-,~A! 7y걇+gOjdgfߩюEq٥Xx1DFFһw/z+ȐGZ:Bpԫ=_0 @ -.MS0궢t2Ul[}wZj3'ѽ[gX֭]"jerK >?TTĞuw #Ue T8Fa5f.l0_Iɹ|O\JqiҞ~HM-hV+,#O)i g0|zې|h!]LP1=1Fn1 B:1]lJ5vDFEu~^_Kvesɓ1[R6ǏV"@V*lyǬD7<F͛S$6&#HJa@S]|Q̝mtwząjaL[wcD-Elv;qL6 ;1T8*jύ#s)TYaڪXmVv3/o_IrƗL=vVrD39d^pUkxw;\XLvADF7/͊^k*JSH;!/\}phbW(Rm3#e#,RUʓ[U<֔z:tӬ]ܼ||ݟ0~] ϟoB#/ߧN1meIj+qc7̆]>/cɬ\uѧߥ]Ua=B;Rng}Ͻ|AR ;/e8>k׆{-v;"Tޘ#TumNƏѣG9tPԈNs3wyGE}iV+~~~DMƻhDu SLyfc fvL`;D^/X: %pJ$,,#ӧOWʖo`UZ,KYHDF1,p${EQ=@ IDATl ,K뜤W=gS9dDc9L'ҰٔOXx"lBUX|T U?ߚ5qZu?ܩ^"rrr7MYԜ3mU#ɧm.Л+ו\PA=P∳9q|sȗ_99v _~u cYr} ,?@ӦM9tPӶ)5MՓ7!֮[w>'Nv9}=1Ln1aѩh޸5#wlVsλS~HHRnHҨW#Q fEm68o盟Ir8FF7t uy#*8O܂Fhe~1%?8v W.СCg+ ^ =W꾕92?ܸJ*A9čQ]wxޠ~t۷n.QFnD4&{Ğʖsࡥ?{|!2o K *?ϭw=LG \rssyd\4 ȁ f۶UfV,>>T*c38o #*eӮ)U~T=Q I,G)o%iJ_Kʛ6@9_ObΥv۳q.ѓȨr&jrq6mH||<Rђmˊ 8zh~/?_~=P(i&Lg&S3ZgYǟ!ekOqd{}' KoߺwjQ8)ѧiwI8zZmK5:HrRYYYu$t:hӶ QQQlΨ;U**˯{'j2#8881"J@JJ %?/OUm@AAʊiX.o( FQQqDIVSXS~>w&Wj(nUYC=ǂ  Ԫ2>F_|:` T'q˩+BXx8I7$ g@h75bjij^.|SΌ*$"y#   'TxB^Oǎظڵk/jAڿVъ83Ӽy l6yBBHEίѥGqO+   T G qm?֑qĀXi]>KlEQya6&kA_ɂ   BPs,f3Mh5&^_1 |Mʴݦi*6 S#EQBTkK5AAA0T6\rsǺeJ8-rP!   u 8oߞ{#P_K7b kA]Ǟ>yAAAk\C&7j nݗӄmR߂  TV8kjekr*3zƑ!\?4 .P~5 p6r75(yyt%nǢ(@:oa2)7V+( z=*_U}^Mfl6|4 sחd ekGG{>/&a6Mâ(18b4ON*hf2Q`Ɲ|iR |rIN\̈́d ,qWAAlMN˿ս4@ӟUb^99Z*P|G2鐟_ظ4X,4XDN[h6YY46YvVL&.+`fF\6-0ʏ+n皌 B N4hhV&?FDxlod6s3{L48Pl:&efb#f# S3ׄAA[b]:u$4& 5 6qD#a6a%4z{䔘L<<CΜA)lpIv66`G $}FUOSڻff3CϜᯰb{pFJ}^SAnv2cqTkW$+3S*E!8$qmس;YUFEFEC: N0 +ffeZ,mm] u_UL~??6h9/yy9GHƸ*$((:M1*gf:kBC9f4V+}hP jX,ęL |n4.#LJS0T{13 PG܁n]D'=fDUUO`捨v;%*6oDq:XXkCBh}IzX .R{p֒n y 8j4^=l5!!*L&@<6Wî((v]\n7YRm^nǟn00/<g]gsppȹ  u^zZgMOߏ,]NN_mw`Or##A[̜"zpX^=n sŶj4QYRFZrӬ ᤏO7sR˨o^0>oGBPK)   : +Ym|W\W9-8Fs=X'#Ny5ClUX`Z,V+7:1aW|C ^H-&8' B-~0ցjP˞ "EQ4MUN BV}~bVU6~}9C[ýB#dVV,^S9yf p~WZp+" VYQs5 Kuƴf-4ӠM&丂#) ţ'×XkI ]QйmѴrIm  Ƞsč E%7V4s%Qj$+ęRz(VВ8__z-h<8{LbTt:#VGqf p?@SUTM}vKc\A־rOjPw/GgYO{s^^DKEQ] 8UE99=).tyR̬4Q' ^ZA@r/Y٠!cWeW,}%g@VYʶ8Y/AqMł^間S   (__TZEvK#^Gra5&!!\Z!/yyŦYLrϏHxem^->6 .v+3fӧLdbP}~ח75Ϛ8jAAHExS *GJC=w|֟7o6}Lt( ̷b!rN=~>>rG/11gS.ju$ ]mʐv{/0gӖDRl b}p0e_P[},416Oj E.ESK_'W3߸/~QԔw~Jbf ]x?3;HaMvIiJew#0kFIjq`;Wg燰:Oˁ&%"թ#>~/^sDZ/0 య/1gf@9ד΀rztͥDݎVȎn=y Cn.-fBl68☜6΀"n dhl6bc4Y#[縏y롲>t D\~> V|5 ׳ߟ^%P'0K.Aƃa) B{*?u=Ĭ[w}9>cyYq^{Mg;% y+3z_Wof8"屭4A]j#©}'7CzCpw}VW]e˾Bj );Ѵ4ǩQfX:)|]rz618B6Zu:6%`̖ ;ʬݵ( {[x=0OA^iK{K35n4w? ljw&k/`΄Hz/-ƿK1[㙰~:A_`t;^Z?KmĽ ^'tP7l. fm盩]s;S^}||O;ՠ |0;쯃c'h_0'/{e|6؟<䯖t!^z+Dm_4>).]ofOIyZc>೵r3rUc|\,!za=)7/N㮋ꡯ.  #{ɋpy@׺h| n՜zAAydfi+ᶔPW_"^;!+NP +{=,z:+~Ư,]_>Y8#q;^"!-̛47:!;K/ߖ3=~իXE+gAi(v}sK(nCBWsܽ9Mhs/G8chBiAqPyqDh 7Ɲ௹{1iɖ֣I`4_ڍi&OOٓ_xCQ(:SY0xnZラ ?.yb ˝o>[8-ǿ} ݚ+/A HPcA7nss $ƍ pA#.բ|cesCc.G/_>k+ &XIwM VjwUzZK؊hvVc a PmV*͊jC|h"řFl3pHkؼ˷۰kykv:k6Ќ1-} Ӈ%<{20+|]a/ӧ >s8Zδ3YK˒3 '+iцnikeCOE- [_4M?ahBRq B]@z .=DEwA"…/>j^<-7?LŠOTvcXQMg^ 4i'ejv;9()a)W#u if񞂸\=UȼVo~?>Kϰi=cjtD AɅ73[~eқaXʔ?%ij۷oc ջ1/[>,s?5M׉W%m˂ZS  # B)$$$A NQR]cF4]ưЁkx#}{&tQL+[Nǯ%\OA?iz'=A6hyCKKS/pQyyYI+!:c4Q|GChXcbm0Fw1Y'׵v"4BWJ>V}sS3_D _6s۬ $Yqm}h^L[|W<ٴ,=OSm,:/6BS`K>Z*/Aa5 B (`$By0.T~?u;OF'}C%G q&cJGن(D^$1jJD )x|0zDFs݃Z܏lyj]4Z$^IGבϽmhv˟!~kޭ:t_d~=Bt]q B]@4M0aWo q1Et!$$$R@)YmAj+V8kW½FaZj׵pinRed{b]#+]{Nʞ4{=F# qѵk<7( a 1!Q%cXRTb(c YYsD9"Ա{@ApG_ x#  4_ YDAAZA "  BfVڝsTX0)GB\4 5# B# \Ylm̝Rq)?t,jGA;Aj1M v:'9I!("  VU ї`DjGAA7'B q$:EQPUU"("  B1X,Za+ R_P tA8  …AH0VA4Ml R_B@AAAb %GA:H A(kPl R_B@A:# iJon % $ P!$( uQXԢe #S6׶E,.m5-i??3KR*Js8yaU>rk_e)Uue y+(mYf-k,ۇB`m1-w}gvձMA~E PuRӊ\ ES#!}[k)[Y~ݨ#i~v$2zW`Y3I:_F;{,*6"jW IDATUF.80ٌؗߎ`m}- : `~*մMAxLEfA(I)zAOqi,%k=@~|tgIkeasVt)Vb=Mt?!V'~IC|1"K,]5&P)"?Վ)75]͠ &(mvѵ>=_# &4[?fŎt͇LeI[aykLPP01=FS6ե,VU%8K{K8j@%h6*`l2SZ.INFa m B $" E)-7kp/LhڗbCvz=1lkxr][^Xmk>ϸul˫CxAax?kLP|oDُg_7,~TJ;Ƽ=yt4C. ;ޟAwY^ﺙGMCWqG >8{^J3XU_ي~ȎktIf\?ؤ8b@9}so_ctPtOfd"vhe   O/m~buG"nMjMNWgSK>Iȭ_dlIdgC olNd{} m &}̧4LQH(Kד0[Yj<;"&ޣyu.CR^a>!3tl׋۟KLzv#yvZ [_ z۴#W>6Ŭ'8"Pwđ@0֏g3{8NAAty6$LLƝV+!ɧaO.83אo81M1oM"x \:lF& q퍡,XU+Wj*~ \VKHh=3EIxQdẾH\ne vWdt/(7df!Tb/j!!tkTPB"Cؓ{EKy1DFir k*߮]z򓘻( N; .S iJ!)4RsUհMA8pXUpbX^Yjbшn|ձbT@f`#v 4RInSAD}؋׃z~ah۾=ix;4'@62 Rʭ_-)a{*VN6)GǹcA v1g:=F"NEK3X=|:;&Ē%syS_]  y"hĎ|Qtm>jA\} 톽o3ήVc;޿fRq,hvhJmlv 4I"tUr3Znў&ZES#ᘛ.\AKQ.j9s670{SVG C)L[~+r[_舯kh   Լ8RJSSn:(#'~i ٟLV3߸rxMJQ46{. ߥBjp!Ai;K!v2-5LƶI[vhub=i0OE_or+tl;Ũydxc}{vCXlڽ!7/WKآ[ZUoG`PcXBԺUpp77[m:]y*YCKTV m~j,j1^פ>ŷm5],"ſ^Yտ_S]#OO#Ůvх/(%ݧ(5;ܡuZ7'@p8|&83(|ߑ^2 (,G HQV8RIuoY;{NGRTBb-JjǪ=qh=F]Ū,{f5˒ono =!K9%zp!:Tq,+^Q$ɪփV"v8$ ?!mUXFe5 I2\r:ϫU 6sut׻:%uĞL[ ٢TגE':,ǵ7':"ryZp.,y,r\K p _se\GE+8WuPƅ` 58L QjjX9kS@U# Ah{}+5: #>S`.+B8}7(g@p!xZ ∼\*'<H*8r apXG@P#ApPd9)Tj Ngd*p5G@P#Ap5G@P#X3E)Δc~^#Yp^"Gu(1$Y }Pj#  *ّ\g&8 DFF8j#jp$L"P LX2F8j# F8j#jp$L"PF8'(3EG2j# F8j# F8VG*59TG@PSX#h9ŃT# F8ŦNfjaP 6@ F@s8$ @ E8'3E#Ap5ԘFP|9rr-GNFtL@ B8} Cڥ1钤ZFX,T+VKE@M~ P7ڥ=ʃ_nܘʑF`T-F gj..۳GYCeA}42#C=rrD'(# % ۷'3#YNg 'p>(#.EuD|~=G߾u&V:D/7lWݮ4]կW5bLÑLE@`Q7hmqdqXG]%r4^Æ=+RMP#yV=߸Egt7f@rpR>yGx#HOњRz(uuաY-~jŝzzWK )*$1I1V*m~dlQz5jE5SV3H5Y3}>MOLWGV2vԯ11>.N_%YYjXG6vjIt׭aaeU^N;p@rrԢ@unzͶٴ+4T##SL~QF #9͍՝-[*әp:*Qa(nך͍ӷ:!}u(?%nזP9slEDTw&֙s~{a(nְ0((>{'|8t:> 9+*: NyZR[U`lz*.SK]. Mv֕{*0n}^TMiH4hP_۰A<^c^NKӆ0=ܢVEV_AkhFnݹSyg߯nukmo4i): 8PGܡW /כ\rmI \}u)Zq8%|uMvIAumnzt6O ![ZCM @, Fv TVAU|Uz ͙:=/4ibt}~UL|{uiAטygǎJTRznv붝;@2W|گ^n;}?WǺouDDE)fS#[ct޽~o!izݺ^v˥sHz2"B"#53!A)QQ|U>ZkDF9ީ:u,+>__)QQam~&٣tqqr_VDD gtl6|jrc^FfdoVV}߯fVG"^]g{EbLLԆ0hHf.KMUhwzaC -q#^#R7n]ְ0ܸl6ݴkWi\$$1<\O5mйiiEm\J7]ޮ]96VDEiETvb-!:Vog7ևY<*JlW7lPbaOj,eWhM#"1"B&$hdz&mZ$@J:=3S1nYff[b-ZhV||e#"zDZML}P.Hә"p5*Yf{>p0RاuXA;P3iQK{X4^٦ʘXKH ?J6..D0Rʨ(Ms{֠.P8H P@l(4rtzf99jn"|R qO㤄Yg"ZUxUbQKRS#?H-gm=K޹SC23d(85+T1hSOK5LWԽ تNT#ACn^ٰA x?}K\?TT%-ar%LNP{qxƍ FP $-xJ' $l7L4ORAT;㆜h:#$ȾZ0vJR V& n5AK,[iiyd(f;܅řBa\%i?>"B͊443Sˢor`CxʲWdXvM6`϶d1 5pcm|jefXN h *`Zw57. I TjQC~CK6Ð0T֐Lqc JܼYT0u^Ʀutώ%ɷXc9Bbu9;>^łk]:/Oa>|>u}۶-[Jb#X$n߮323U0Tw"Ńh9vE!٣k - ,` P?-2IݺJ3ndSu%v>[W.0t=6qzN}%,LQR7ocHә"#@@ˑ &gH+"#JFb^kP*\zaR: p$_Wb6IoR,ԭ-[U :a^EZo֭-4~^=*u jb/iH[֞V?>N"$elzaC]ٶҏR_&$X~ Ex2ɮO ];=Jm֭SPs DFȃ 22Rtqii+Uժ=!!Zccsl^]k߯xGi!!.>^$&tU6w*fSݮՑZ`Kjv}joH~՛ >n׃-Z tfFzdg+@1$ݮaaZoaǦp]оLO./O^|UivϣL! '‘c\#ϱa~D IDAT2˰DfzY*/oyTԟUٰ,c<kJcTTժ/uǖ'X!-\xunnᇥ3?_q4qT%&J87Zpf>(=ԑOҐ!Yg-kfJi ߤ1kZzƎF a4eK+q GhF+@_#m,}gUi"`;ғO? x_\Zw5W_6-9Z.KM5_+VIfңJ&}!C@ğC/7N{,}t]Utyz$%*jfNS#޽jx1u`ğ3k[|)mt%}yU*VoZ29bPA#@måT~tuUw^nlR{LzuÇ-eÇtJ+Wy[噡Ky,FI]*pmf,ÇK-[V͹C喪9s.,v/ /_W.2+VHsV]P #HәBE@E⮻1[$~}m /˛5x.ans4}ٍ,a.Ҽ fۥ>+y~:F26ٷ,f`R4l=SN1ty7F~]}xf}*Yڦyss@žyF.\(-XPrJ'EE]ިt櫬q~AZn̠#.&MsΑƍ. J : hsgu`Stt飏/~}spbb̠wosn+W6С*FŜ9TS>9Eu*(Hf1mZ`<7(y-UR"}z5bǼ첒ˊ\Gs#@5!j{5-̮-b?\rwK!!;\֧Rj۶~U9(k"4[sp'$,k||u*Ku9S̳iiK/I]ZndwؾNq7m22EoN:䲽{K.۱leҪt_Hk֘c5Oyڱ#}QÆf gR.f  t_ήziܸQrY9%9mSՄpvۥ[wޑ 1,kرjc3CUSo_e s֛T@-G8vWHgAɽJaaE=ۜ2J.[eS7hlE-^\rYӦG?k?@8H76gqӒyr7 /T?>u9eJeyo.]nHn@E8.ls˾Vq g)~ں3 ᇥ,n#׭+~<骫]͖6"m\1r#7G**P7ǽpQG馒5sJ ;Ϝ2xsL\wK|?)_Mrr.*cr׸e4mيͺ(cmsscfdH)Գ4iR}z :_]ٵ1_;wxF_L@ ` C TγJVI+ѩ4f%׭\)v[Ŏ$]~yn?>LV3fHD)k3b1+ּD6i`4kH}}3Z@-tHe/@9WT/,իW3RV۷@:zGـ ;*Pa#@msg`+M͖#f8s6m*e-b?%uWZ˒w8GfܸSʖnnU9ࠕu)Ғ%Сm_N`۵jeN |xyZ?m~۶[7sى/Ҩ݋>}?o3sɎ"]iki 8n"gI))Ү]R~nf@%%{,@iiઉf0p9{K76Vz5s>~ds԰0J˖R҈g%}e:vs|Ds*Yi vfzd i~)" Mz2;֬G@aCTE#Ho*&4TBUYQ}{6_U}fͤ6_q}Ç2bb;0_eWO2|U=PntAp5GǕÑ,'8"ZL^*A8j# F8j# ;#YNg pj3EG2TG@P#Ap5G@P#ApP#8r:Ss#P8)r8 F8j# @-rp&ftF81d9icpjdG IBT ;UPh9j# F8Fq8tP8fGp:S`_ZF8j# F8q8tP8& Gl7ϋiѢsW{h:^gQ?ΞoߺU}s'3EG2T{~]B/YQ>y/װacuϲzT:jPf8MS?.ڨ^pTDH m~-3r櫯Yg8&PuYu|9TPPc$8SH@˜#}TB5Y0' 1J9=-GSFxh9*6pcjUk8 ̛3[a%ؠ[ n5Hf:cb~6k * (hϥJMM-s[ժDuSIIIT+\~we, @SG7oү ꌡCղekl6> Ð!C|ou٣_."hD%⸡[ jܹ:餓Ԯ]YVy=y<n].)//WZjР4o*N j3k@xX,& z52=[#\GnZA_?n[=:$F.]vE8rNS}i׽u飇VwM=޹]nV{vjޤjmnuyϻw?$շ[dXT`:GgbhUv=QDuVט.<|pu}vCՄ9 tyYOE_1sAk+orЀI2k #Un5F<>~ߧ"]_f=KcJG[9_gNܠ!oU)薨xCZ@|饣NzZ5.SS5ZO,xkuZ"+ur?Q#vI OS'\?5aڱ}[BX"@ [MVtf}[^A?G?SO=%YqWKO;|R'oyS,:pjwúzPiK<lM7sozWW>[߾QC _hr0uu}$խ}u':/s^ȡV,Hzzpmn#7tvȖC|I.wc:Z6 XdI\ s+ջQ#$Q~Yk.Qj֥zaNms֪ț55<"_+_D}[-r,y9<yɝ.lָlfǭuOwUH=J)tj#s΋yguyi`d);ũ I$ |M U)n5G҄{+*}WubQtv: R鯚wk0u?ް_@nmx,u,? oQigSn ->@t;xDz1Qn^A[#ZTgdeiݫ%gjw_rn[Y:4zYMRJ=u\Z gj3=9c.Ϭ0dSKNmmߡM/4qt Yoݛ>+cuM` >e?7kռǕ}"eK/M_.xWZy~9hk_!J7AwN״ՇCtY|}i0J689K7I]׾,NO΁ bPsoՐ?5sum=fozMJ2P]1ZI彳ZZpn`})Ǎj°jۢ:RHCmZcevxSPc5NM92_yw}G;*U^ϼ_!wO!:zc\~zFKdígCwKg}pK8wk z6j$D-z>Dur$h_:۲D-5[k iOg5tU΍u8rQ|xTWlܠohY~w.Zɑէ׹ ̮1F~̴HiIczR7ט ݴiZ/ <\},RZ^\6_-nƛ i3zdZ$5i0pڭUWٻc EOąլe:SMhٯc=ڒ!7(e3Mհ=5y.yMݒ8O}gW~{mؑesdvGWnoCGb8Qڜ$'遅b,M?<|y YdȕR|2-;eJiQBEƟRq(pТҔFuz)Xˏ ~n*@a*b=^M}O”4ht Fһ_mBkZo"exlG O6U4OXRʵGU{7}U:%k/ @8"VYSןYvi3JBV{k-Ż45A85I mKZ5gWK_Վ闪Їπzy+ab}uznS8cZKnvեbI.e?srn@뷮O\GW:,9ҥYo6.Z^N7*cԮs{/=[Y)'?szWZn.B%?#^F_mt %r,޳)a7ї?[rVOZvz_ai򌿴a9ggmϵD#ꮽܬ;Ո(s4Xͺ\gf]^VMYU )_)G߿>ep4@U xhМ+tex>\-0]=1uropjokv)tN|Em:֥S-zf@@҇ߟPub:/ipXG=2sj{Ni7RrDEj?N&>SW,G& ίn-դI[ wPI^Ιz_Jծo_\t|f#*Q&֫֊oVVŚ-дK\z{9 )OQJE?1kmjC^ÐkHP+ܐ<W#P:LQ2$FSzqɺu1]O\B=2Yy}|>? <0/?"ؘ)&jBh=@)5K #:V8 T[ E!21HQL *UIW'`Z”fyzY5.&QϮ2jMdXdY* "#P$; FܭH'-X_EQ0b5&QA+Ҏ IDAT٤'j觩2cx6OE׍ӌJ5Uш ǨkUl UՆ T 蔯[ES_*ВjYWQAX񦎺_vւ2V ֺ\'K *R'PkYWػxjY5-GtۭƞN=zW. Ez}rjm8IzZ&tkC#)>"[B]&%KԷuZ@=/"g?Ǧj6X4JO?_.~ ܯ)Ӵx o>N/wR"Kk|EյZdDлVW5_á2^7kq|j?07͘3cEG/n5_4yTkEX,=YNݠ+gޑ3uv\}.w92,vu{b\jRωV)?w⊴`jVwc:Z6 XdIX?~HAR5lxY_}G1Ɨ7 YiS#/ԇ=ƻv-՚gӫ{||rY;)Z5q%qF߳HGFpiꡯc4EZ] .Tɭ vUu}25>gW4}6e6I)yu 8TO6i?nIZJY" ʗfC#;DX¡N}təzh`d6ξC)w,C]-ιQR?)9}s_яQa@.z߿H^w;*uo9 x)Sc5zAוt- ;*uiCB4@r|*ضD-uGíBSZ&S^s e o̕{~:I;S)jxF_5S|ۖhRN9lZ_Cn*eito4Df=[*:-<ju>NyK7>W獗g(g35 -Gtۭ#uhv* [o+`wCuD,G2dȐ/+M;IF3Jƌj+OzEXe/H \e*]V6WASD)eEPNq .;+ouoGC%K4sEzmb*End+[Cb=G€4^G|U5TE9-V긼9MMPs 'u~jK*䮭jnu:ԋQJ2`~ @(]X^J2j#w_ƫ[ =JѪ#bWz?F^~~MQ]*Tް=ؤ?l_y<΋ּ+R+/Qc*U28ymJ<)i.HWrltPw59G/լ+l%nՔUGcBvSK}e(ge̜G/W;_~CR?Wz“c:YsK\um\!-onU-ௌOS5:׭ރHu:v|v*8>(US}Os yk/j}"|Ūz~m']sAq}ccԮs{/=쾁׮Fchi醯wS2Y\3K#vLf&hؔE8͸Oj]*)@0^o:Jkmpll%ګL>pjokD vTۻ\&4ƺ?SS۫uxMQ_|# e3{I!)ɢg?{Eշq-B衄$ R  Aǂ(6DU(k!gs]s%;̙s~{ծJ="og1<$jMngeЭ,C:=ʭay]"1At;.۞ N 5yI9zwMɭyGw[|ÅNHg2 >|8+WTn! 4ҳ' FB|Nysf3h0eV:%ab갭lxan6wvZB]/q) Ӿ{oOⓝazuvAʾY_f#>>KvF&N =FDDrccFe=F*k{g`h8㈎˜}4kY`)IB1~=R{3wͼŕPO@)F6Fm$B%)RVs 瞼dӘv]5;>d#qLWc~BL]g\r2?$T#]E6|=$ز:H9YEDDDSOC)vQV#"")8"""FDdDFFYE84zL[;q /t',/"RpDDI5""Oj2qce!\?$:J$Ki%&W>ʴ])g:9ӳԹլ&Lp]a^>hGGDDDD䒣n5縏-}-LvTOײu 8{$Gq}T\r ܭHbѱ ɄB[;@*[>LFiˬ5hdζplg|s_ڼ6Sd7FZ˙]_b)y۴|*}ٙ܏#I͏OQ%&*o1M??c{t 澄)/iJ Kv<8}csIӀ+W,'iѶ~|ƄUy?=RQO  S<4mKJ{QpDDDDD"Vn5}34~}G/g)-0Kٹ?[As\8lVoZB~dOXUSMLl8qޕaK1zV[O2]|{؝O kC[#Ԛv-wJNؓ>63x{?_O Sױk.L}JF:Fܿ,e7y֚~V#0} ѝOޅOc37y=>~k_'a1_rUEOD.%VeHx#E>(ODj\6{ijڄ/@c/}LW~zh&\IK2<Ȼm? 𳙌pF_agGUkbkҔ1?O}iQLYx~ a=U~acͳ4I{zs%3F-bwbzk7'm5ҺbLřGb^ZGSM+J`h S$mScmTedNmE[Ɏ݌k]dӜs4<I~=iz(‡j LX*S?$O盺k\c 1b˞h"Sk%ũVTXFPE+MU*P/!}YrDDDDD.9K/W8sT-8MK~a:z^}l ?OĿҤjRV ~LHL-~4dƚy㉁δ92<|yhCjgn鿾mO/i郚YC9VBbLφqr bf+>lvuq] +j9V=1Q!>zG;y7ͫ'iW L6Bsjg4Vҥ1?® 3aIY#fṤ.5/+?RR8ufePz\н)stcg*cXt2OakO=c%kN1Z Go Խ~۴b)v]̓2n5""Rn;""RXVF57ȑݬ9-ڄV eY4N }388V+}p߼W77၇"e=K譌ȼnװ~v=_|[yy`3OJS~k 1}q@{Z;3|rϳ]97]_z̕iqfDm;H ߰ p9kݛ$_mF}i" vձrprDDFD)LTD,pd(i@hh-_3-|SndC A[>???y({~9$ =''4 k@[H#gp8| |!x% vv4ߘ+Nfc)x~`Θ,rZw͂q\~vپ=Wi3Now5'nwaՆCxK(@5=C悥 CmOY՚|8+WTn(8|KŠߗg74d q8N͙͠!J7aoӖ`a%siɤAzXT$L nwfl6q,_72qDe )v٨.!""R\*k+0""#"".k@D qGDD'[{k ""Ո bODD.6>(儈HR)8""R" DEDDDDʔՈEO#zj\r\ԡ?#b+:|īE=Fc)w)eċ#F ?^];87X]2y[(>Z-lxQ*9\}iM‹8}o +JΗpryb¯NN[ή(mJY( 87yf-lZ/quO^ۂ^wUxO͋.|q7ue_^qw_I$u=5<n˅tt8HMM!%%dHJL$11A&eN|Xh޲%x ,O1ʖӹQ & LT3_sFlLB2)meqG˷v89.e7?>݃FLL&W.Hf83//K}cRN߇ɋ4Fn~|'MҎ[#_4lW,'iѶ~|ƄUy?=RAGDDR>H0N]7V?y);g7X;#WH}fo־O؍00b䪊ymUn(:4 ʰ=o+D'.=y5w$l#Wrž |ܓO^?҅;_c]( YǗnhEiXkv [ql[툈#""R4ol<٨k4dӸ@c/0g;ԯOcI3m)(\8Yޟ}W aF} u`x0I+ͳfmOrm=s%{^;-ۿѸ~c"{^u;LH*X c43k,bccϛwΝ{Y" H0ζt07'm5ҺLřɵ_m.\꾕H͸~i^L&L>yi[ 1NB_h^|٪6/-JL'JݿFKF`.~ ,ƾ}3$g;Xt2G#"""""E:~k\ə ?C}$\wSŒH:O)ґ/aklG>L次yYbig@, ֣7ˁiӆݺK/ѳW/y9U)%RSLP 8@>M!?юOy XL.NӮs3O.&{}e+i&~XqO dM~֙Ʀ^g) 0w ʺbсʴ~9RNI/W^y%&M+`„w1U/+)""b[:~)nmG׏bއTÌf˲h$mgCFk au\-'8u:S~d>\EXz+c2x5<dV;&{1+'% e=?{{!jW|20ά;i=v̧0rlw|=3WvġXa=*^ G)FD ((""%tlՊvs%ۘp|K7Fv<Ԉ4k!^~At;.۞ N 5y&XTuU9y<GYXw}$,xHba|BW`H٘l!.d%3 m~]|`}|Mlߞ+4xs'ڷQUt,@2p] 2a >+W*7DDJ(EĻ}yyACdޜ 2L%R@ӾCyx0<pyGlDZ|b֮ĉR*baXzUټy6M'eJDFFn{uv܉˕w,ٿ?̓2_g+ʺ?ӾB Raaas%Ξ'''fBBBrUԯ̓2H ؘƨ,TD}c?eHeID GDDJAOd޵FvDBnjrgekKعM %Pc O}8BODrjcƌIDDDk/FT{0e#**p"ic[REϞ=ٽ{7={TD#""@#sD~D@u$V_uoĎTHr͇ ZVCa?;/&nG 2{d;!f _Wc3`Q\ڱOϦop(J1mo獖h33F؜V|-0ϲ>=1s- Pb`$gIJwS`{wsiB+x1V<Ζob>}jokB%k<R [͆ ݻ7cƌ!..1cлwo6lؠn5#t HaOeci2?u3xzd|)8+DuP7_?(?bl_ w&L>'D@ 4G@"ۊκihDлsi߯}n9GNm~zlc͏|OyN-b <|Ul߁y4weMyy&v\緜;w3;䷾'_`Եl߽?0*-L;;Hv:q.}2eNv{ 딑,HP0Fśoc=F xx75j2GUY "O<\wK["RH~N[ruhjڄʈ=}>zzL8z^nOxoAs'r@L4㕯[ts2"0 l{Y1W8Ε+8Q"iU?<}u l@[dz=;<}x<:t =Y *4NwяV=|{q h'sYu,vgJf^ |!̞=СC;"j9""RHOI"IN\mD$?~yIq}=?Ul؄J B+ىI`?\: #Hkՙcٴ nF\swl'n߈aݼ~5r3}%jLUTb<9pq> Q%;0>6-mЊs 6p`{.@xx͹Um* FH~E* DD= twc{˯6Fm_˜a |\N'.W.:4酻W~[yqCCW2L 'v{jZ@c؝n \ToӎnI$_+`hU+nU^իIK`bwv\qdJwz.7\vT<@Ŏ<)KRjæ&BK' =X 9Nk?xޱʍsM~lJŞ. N;vO.Yerm7P$ Y[l,iF.IkYЧبUȥ&)-4~7ϼ% W}'⫍6dz],ٜqt5OѪi &DPu=:Bk/}iobTHmA9)}6HoUiyN{ħ[#Oĝ;}vi;m>0al>oS'~>1/zz|믹[5 (8oy}P^N:,**8y+ y/h$s[ETF/ύ$׹O~PP^(^1rZM? Vn9#{9ls֢E k:Y s K~{/Wb߿sxdM ztsIka3hv}cΙѺqY)tdy-cy#|Ϥ]͸؊/Z50v_^yOf;v~k}aF%|a4KeնphJFDYj7m+zSy#ʶ_~bАa$Ǒȏ?|?P y OFעug Ͽ޹v:w&uwS~A۶mW_|N~Z:q,_72qDe #"yu[PC2ad-KI )%\$ 0={UL_<mZfXVeQDNDDDD. qr,e nv;pğtG=l*B{vѮC'bOK2r$RSoۦ<[mϿKk9Rԭ&'g]LbՈd** rq8oKߗgSˑaCMsYǟy6Q6ILLc[5ia֐'O˾\.f#li@ֲg'd^,o?~xV>9\֏~ÔaŃ 8k e4֞xHl*{Ky%WF!\j!."U8k^d2cǨQ2Dh$E'g9xML-Yv ǎq&M#^Tw_6>ӧ} 2.J\,mu"PffqBReDňylFwoxg߃IbqUhLW?ĩxA5q˴ ϰ/!^fѠQcjr"%d2c)0"eN2aC̻*bs |8O~B>L=8l.Vj\֎zv#`#IЃ.M3ԓɣ)i'fwkԲNxAu#&e;a5-X;IxBntyٌfWCwO٢ "eN2zct/9os@Cg$<`x=ՎM:]tˊWo⮡8MkQ/C1E:`'T/)Н\Y忟r4sXҋErR{寵k8qxV@kbXLppZ ix$LT>ErVlPbƲQ#І6Lf71d싳8 2xghpgFLN`+f .ʭxgs߼9/1Sv/&Qv1pݻ˖ѳwj-ǔNyt:ٷo˖,r`nׯm*[>5(8r J4s^xx&#oL-%]T}Ĉ?)ċ?C R~VF]N'^py+J8jҤI30L^5VpDӻgj. =WDDDDy<N?NZKix<&Ñq:&X[ʧʧ\LTDDDDDJj V4mذVlذ!2ͦ )u)o[BTT}ݛsYvv{ !>щݲaKT+z)R)w s⹑o_СC/3gګn'*j=k׬w]w䆾+iv\jMKT)8ȉ|$GDDDDD )j5bN,#}q8)hHHj|w4VC&bG/3j@"?`Sk'O?aZ3^35yx "q jxMz;%b~{*oR},Zdjї'y:URd$}ڋ<Əl8k0M{>7Oc(ݮ~˖.p)|SŲ`#<1ȟQ)"%EW$҃+Eҡ<7neQG6`3h^\R~N`?)xPBVKj]>[&m:!62vԭFDDrSKTM/>C,5|*?xoëd6pF_K;&;5w]V<%Kywߞ/)HDe $lPr5$Juh֢NDJlj|n}P#.[nPSSSu{9 /bjDkځ q7yo7e>[*"ϑ3ڈ㈎РB~O刈H!e d1ԔBiu#/ x۟yAk`9K*5+aPS7[K4 vqp\j@ [gS@ w$hOub+WMB6{O(8"""""RҒ72eTryb+ >y-W3cu. 5m}ns0׾3og{a_ŠiGbqlxwx|C;2p<ԧ泭!=|;Ͻΰ`x6d|ʅ_i&S|HIS1f)=.S :ۿ5-E{rР߿?wu3k,~zf3)))^rr2K>iܗk&<&3oߎ<ڧ.}FMϨgqoaL?XHHLɱZXӌ|L1SALy"*pAs?Q CL=©'iǣwPpDDDDDbṈkM̩$Pn}\PBy ʗ^|]w2);vq8X,tة\\S}JuN<<Ӿλ2zawз?7_mvWr,1痿nw`[Λ@|O???DDOo݂LJkuvt:%^!""""r0a78z:˗tf䳡%[!-ŀyLc:@ZGɓIۯA|\EU".sOkKʴ|&2oB+vq:::oͽƅ_Ù{_Xme>x?Zb%XlRO__.#N C>fe3zھ6wԵ>.^z'sm6rJFZbql֭mt:8vOJj9"""""创#?+1DTZcSxKhl6SF :H{t) ){CY """"(Go|m)1^Tl\ޱK,f% }MSNv:vBDS$@DDDD3t`G +JmFocxZg)*UB۵4lА$KOͦ$eB)WA]x[A&4xqaŀ%ɄIl^c)x<#)GDDDD0swdnh ;*0tw|;`\f+ר*>&j5 'R Hr~!wx"=a_Pl6cZ-=Z)=YcEf6^?ޖv; ("#-yaѯDb΃7_ɀxwDT6 ˊS X&""""k*]_~1L) yŅ,>ÖƔv*]~tyVC>!3s<6Gs:>M-Պ@@-Z^nޖ$Μ#$&&jſbJԪJhݺxn]GD.$8g x/]|2v=_".EDDD$WN/}˛ek~nNhŽv9r0M=kR*gt:ٿo/+/0 4L;oJ"<@m Zcۉ9um[P9 uHOiJ( IDAT['z?O6`ˌ#ǛhUa1-ڲ?LDDDDf\۽ 7)cl64mdbͪ?];U1 ǎкmkvSOrbX^&`bԼ]fɳ:e'DJYY "x;r2i9)snOEM1GlnŋIJJ}dsp3I4ҙ5PA #bbw885]i֝7A7AƔEiFzrNޜnZW{)exݗ1fd"T4V޻ԯlbxTB3#مƮpglk G.>>6f׮]tڕN*%[6o`g7H-VY$%) "Rr)L~`7ox1 wLDn |-\1ēZFԩSi V᭘`]w- )Y"_s^+FDJ#"""""p9K`$]R™#1qc] gNTÇ ;tаʯa/zzQbӗX G.. ʼn'p\8Nn75k2$;ࡃ\.֭ȑhf35kߟ={vg;NA'yf~۲>?yfu1GDDDDDrDŚtcyq$ܱ7 ?me=NNo攞ĄB#&)) E`` u9Ixx8-Z`-T\sf6mBfMR% 6l@͋|k׹-X Sz2LJzZkWp{BTs1f6BfV3˯%%+)6v;5].@~R%ٻ7ۼaaܒH$z2ܖ7?.zeV UT9c(8"g}p8}>-n7>9'ΜW}^SCӿ``HXZV嘵dÇ%=6u\.0"6Z* 0{i琸8^!k 2Ŝw%W 8=˃qqpP㡪N{GΜa~Ŋ !ƪ!Јsrlǎ5\ަ nljwZ$d^ȥ^$ qښ6/=x7իSzu֬YCbb"_*lp$>>.[P^*u;oY))0UGrglΧ<KAۙ\LÇV-6`Kz5%Y<\{y_WQc蔚Zmz%%Z؛O+2 E;jʄ_'Gвu bNNj{>Knc0虶<>Zbh[ǝu%kE9x 'N񐘘6 && ^g_g@@ qgI02g{+iqg2)1ED0q[w.]lB;aѣ %[!E[GshM繴B|e^0|x#j|y7o +""^ɹ[D$ &3&s=VcNۯd'sرǏǀǁ/1gR`O`p(h Խ$0`MEmQbpac[QzuL&iiFX,aޭ&=ʏ*?YOִ涎QpD.@E$~}6{4S`$`jѭ5\r~? .L4EnjԬJLW|\%`InAtqLgᆪg7j:A:JQ~ 6RDADWQi*44- @lʶM)|8;wΝ;ޢprfo^OBBZUVaPB%rgg$ 4#׮^a.˜u-8bnBB!&1d 1J MJ%Zgd"= ѸY#_yzٙdg`F# :(F55@TeOix>0..%rO܎*ck;sLJ\\5&gp |Sz܋AyJ$`"$90r=`*ժc2Zbz'=&P2'tpՔq-U;ͳ\9r+W( +W.NnWIhD ;<&Ea7nkxm\&u}/GwEaOA"7aswuFܲM)\QbUW1(rFnOO&:z v (+`"!]ɣǂ))}oH?W+m3idx8|3#`H̙h={38st# P?/0(43T Af!GD1e^\>қdyLLYΕRnn9nZTZݽUaͥnq,F9 7Fnfv/eyk֫VXE)O)7U]]IpDFp'({8cW-mdL|do֙l7zc9\~;wB`goWc%rƷY5 ͥE/Gi5F .oIEb$Se9p#nI%c#߂& B$FCZ?$t:UAX,]uDQ4[eӓ))T(uB! ? 9y8)%6G"Fà %J#h>+2.!));x$9^^noVbKz6Khӓ^^XNgy`o]t ݄z ; ,ssc_{Պ;E!~RwؖMy#c!GlQFHmQgrNRX=4l@tjLYpUNaM..lsvFk%vqa-LVPj*22XpWUR4jvtd5<Ϣ'D]sݲf2`eFL&~uKh8׳ɉD `Gj*uM&B}(37ԕBQRHa5sD$U!BQn0B B!B!=M#B!BQXhBYȀ,!? !B!neKahR% U[7H]O4rtf5/QC4 ڴ;OםLђrӳ4 X!Bq y:m>ٿ1,f?W Vl92}m3x%q#Ψl!RJ%CNrQIpDCCsRWB!(y#vUmh}`aśsЍ~LU >ܔUBP'VFc*J!,YDvԜ_?ax4kF,0F24 YX/G![I^ 6ldH:ܑ4/xn3E맾ᴉr )YvO|0%«~B!B9m`sti7Z\ i׆рJy,{i伥xQ?_1Z3gtP65v= ii8֒f5fTL }}AwjǹZ/ }k{bߣr&Q9mXVǫr Q$8"B!CA|F?+_Dh;_'.8X#9}$"l`r Dװ|n;ߠiB!B6k?ù=8NSx2C^_Ǐm71f|8WoLR\o!v/>!G-}g4vXN ^͛ظt,Pʕ-tǺeN Iĩ5`Tf<8>gSmfzD놟SWӤ+8~i=ҰaC6[Ҳ>5k(ЮmFG;#M!BQ'f0>x1W7c4 Cxș=#lXV d>/ }P=`F~S5գ~\Z2jŸ|b{|ΖU㭵ed l3ukaxl IDAT *>N6cxM?-q49 s+#8ׇWtqmWR9d) PߥiB!BRf7~$ժW4(B&s,\;nlyh;3͟.1י/䰹/=߅h W|\ jڍrZEr ^<ѓaͦy_ɂCB,8TQKoV|U536N.❙G-L~ GlN88]&3i ֜Z ~7~SwT]-amʳ~eLXnͨdb ]\+Ugc·_fh|6cώSϼ0ۧ_'J sxiHUc@=zIǿ?l8[3'>8O CZ2h G9wF i *,Sl6mԨQْtT{{T =~1޺KiOưr|o' %3FS{CqhZӄB!B!Bt5>ys8ZP4&v,HPCz{f\Î2"ƨ\#M!B!gXB!w'V#D'jDI%aB!BqOB!B!i2B!PV,|۵<ͣǻшlFUeؒhppttB!Jg==W۰`66\5`KZj*碢ػoW$@R[ߏf͚jJQIpD!BZOPB3;Ŗm֭AAAh42jXVټy3T h4JQNIpD!Bj@vWC{=Z=@YVN?n\x}W`x+ &la|~Oc|\8ߟeͦ'ZQ3}=?K~fӭ[77o.TT ???֮]K:u$8"D9&!B!n0qt"ޘ㼉Z8F}_y7ofh8m^e`@1,ya03冗h9ׯfAC{`;V*e+E\|^;~$\/rΝ;8qx:ԪU-[J%QIpD!Bxt;@1ou2ޡkh/3m[\\ܭ5gΝjZnQ~RB!B\bP\zo*JS1FmtYvaF][.$ZS՝ols֖Zu9΁v*7Pˊ=ko$8B#B!BE?p br-7ճ}yڽr4r MݝJ\ؼzPOo3%|ɜ?ۍhD~zPz7e[C_Sh d\y0+Th;7mé(Ո*l˧gL+9SFnz85) ]W>͛njTTU*b0g^ E!..ή|;oS2%>jUQP7oqv]uťeZwB!Lxyc$9FW.HB?o_N`jȐ7vGzz<c)\ْ2e*`IlE ~Æ…^!&4 8+k#UÇ@Ɓ̜{pr~c&3giyPbֳ;}N_sW&/Mê66Y>?}j3ry+qmF_`6:tc񆟘LF@b1H(qvZ41GB!(5ymE^ys3>:rT vnTb?zscѷo,ǧA:h$,Dk(mIqLӌUpށ!j砬/67 !ņx#LrX!DiWV%4FA/ߌ ޣG;Ӏ^L}Gfl '~1Ccl]۽ĶhXŕ}2~je?R6 :2qDZD:O8cB{a' scsp t 78( Vc21Wө[? 1Gb7QLzgeiiQ#B!r@Gc!G5W^KFЮ/'OQ͕aese2ӆu[|‘_1bbcyawhz-%Ӯ|l5rֆ'v$Ao :? 7Q=I5Q@rwy;T4g( |o*4=t8x;濯?X?}TYsX7zSB#Bq7Jy|,(x_G&58q;p@`5̀S}\j gKոUv݊@qs6N;lډm2#L- 47zvUs &s6_j16lj=:|p¯mi ~9Tq9ҫmZNT v!dWP^zUZ[P^۶[# _H #=ʩ,܄_?fut懃ƒ~gҳwW\xj\NCvq W>ZyZ0 ~ɲ5m>Ɋ(zaJ|{;,V/jc̞{i+tN\G?"h|i7|4=BE~ Rm/2{"?L]iKκu rW ^|5G>MWkaOڜ"j9SkXswp-0{S^|“YgAt=&Ǿ51z~9K;?v+L !9w!>>>j4|}}iڬ9UVEݥtMz5pґ+u㥿93.g=6K&>Om=`s5a@Z|-f l̷x=jZ@DeϡȰ(v'Nab Uc~fx(O W;S4ZjMa!/-zYFBILYFf2sf*dLAԲDhի9@.ص&Q\N1u NP#Y&ơR_SX|'ӿD-DџfxeKcL`HF{ʡ Cڙ7ܹV#-x>Mz90`L%քz<8r8<1H3/W} &a.WsG܆5 ;śTWl^$e|Q]A0|\HW6cC!xI\~}r+1 LXs㨚0~]vl$=ʯq09mN$KVPS8zsqBb~^X{ c%WJXʥ-Vo <&X *X[5~wdv 6jg%@F)-S0a~N9SX^>bŽn7捣(mWvޏˋp]5jMu%%'Sn}^a…8AXs`yX w'b!P_%*W \!(}ׯ-UVłͦb4oOT3!Nh.X1\N \FRk2@FBP86yWOK .hPċW1_/W+5'%4| Sl S\6J!Ѡa#tz=ZvQDU6,/WEwMKM\T{#>HJFߏf͚jJQjLljISK܃hNgsN1ؗby?`CfWl*T{$B!  5LF#Ee6uFPP,[V+l޼ʪABc-8bI"V:8x}X\SP#PACƣ2U\ҹYO7'7Ej~|c C0K!QIљGt>wP0ؑ51DQK(3r.BQ>27F@ؔAS~`Z`2βrX>q3#2}[5Yk DMgHqm}ϴGDk>jޡޯX=q9fӭ[77o.TT ???֮]K:u$8"q A<xh-Gѕ''nDeWH>al"O^] 1&J&sbO =0:}cU5w᩾ZEkk_Ϋ[Ӄe[V+ЪWsZ<&.z@_`_jNN/C:z,ʣ 3)_-8}" OvO``\R˗@•ҽU;wpbuU6-ZKn hpD QYc z-Saqz/x|WFQS8zsqBb~^X3`QSQ*g`OMo`[:Veats21Gv$.Yw{s&~潳/f _sPS# 6d=tfɪyz?d7wK!B]fXcs:,ڇ0#[Z 0 z ? Υp/O!${iӦ Æ #66/R+ɓ'xgpwwV]ޣaa=8rQp iBk,S*zwZ%!Bq36m3fGa``.CL2aj9#=y*huB-* ʞ:3{lzt܅ &_!6).Szz:nnث!%lbbF8{5A4ۋ[y.Ɋ3B!J-aD>榴P~cſ OgVРՀͦ,LIl_q|#ץjQ=3Q$lFt0 :?^I/1vvƦN7͈Yubpl&|VC;H˳N\PnT.HܺM? (tžzMhpW'D7 Lʁ3q?:%4iQo5_?~6 ̦Bn,(X,i>Nբu&B!DzR/ 7HI`g: ,ũ:MZcrIW vn{yMO]d+dfÜnDu9+b>ǂfR~-;pE-R<^!ҭ)gZbc%=-i|:}*;0dT^T3'-1˧9l~ҽa_;#bAUUn_$8"}pDE_ա|J} ?A؎OjG"2Ej sq0.ҳe+k$FtuGc>Srxo }雚J{ xZj&n  `"WLCe+Dn&N/qLonKO0C6`i`ș$ ?dTNc D _غﰢ`5&s50zd(D#/o9(Xsmy=տ.gě6/پϘj!$ AMpgTBsYD4seܲ49W N-ďǍw߽U@JhP_{YO&z {pT臫2k? ո'|`Luu40^lZIz†ŀ*aI?fw Fm|g t+4}a]X dV]{-eiț| Oè/ۓTH8/`;w<{][Lamx"iWN~LϠ-s.6 4| zeUUŜF#d1}||NJuӦb͚sĚ7ȸoܼؖ8|$)ZZmR:ַѿYy-x;a7+, vGRDX3ɭ~?lľ\٩x5ȓػ+\vyi3 Sa_H:'ӣWiF}psCrkf癗%uՄ3MPh欼 ذWŬs'C̊yo|xῗ} /4`|F`$<32^`ѯ/R%3vb=@۪Ly#Lp?<׿vt[Ej׭Ѹ=ǜmџ'R2wq@`5Pn#TBL2i)="Y-bLrR"...Gj ))xwipDK~_u~0Wb[T3ItY)N0vfX3#p!p2.3=~Yk0C*F|l_&_}J3p S霺7M|@aKbkHl0UuvϞ-3fՆMKM:Ѱl+)7|РwrjJ#X#(P;ZjmVZhk1ǩڵ5j!5 (xx(}+Ws6ف&=yGq{k'6(8UiNM,KΥc+`}-px<&7d/MdLl߉?يT|b`HK:oݘ z ~,h[KK`sP z=9l4~ctWGcFv=eY|Ӌ'z٪FH26xd[Um>(cя˾DCMG+0qr'1EO9?mDtA4'~;]GWN{*6L.b35;88Z/n}qAgtU,&G^s‰U󿗏ѯwm-KN]GG2+,;0z+NAT%)9Q5/n/ǫNv$PߕLxVQ!gãI*MHE^ułs]ŵu]ݵZWAi" *BBBz2}c$3$p~3sv{ιd||/"k4~S]8]T7Lwn4\u#>5uҲT# F+Ѽ-hv0k@05V82Ʇ<O4lfRڵ%!i_~]4ZFvv3\o_OՂRӂjOs|y}gN|Ktxo4я\ɫzS"y3{>ŧ9yPKlr=S1;ڵkXpVfڵkO.]0rcrDёݜhx^KҮrqjdcHtKdQn..M-{BASmT;Vo%V۝1dRBQIa℠X饞[X:bmNIe3sW 6g*+-xťt݊(!*Nr,%pk:o8H⿏:X|8bmI*ذ6:/lαWيt%WaKYKh Qę4ҟ_f\If섻P{` >6<\Fv 7uیE=:~M2;sy=Pw!eSX@Y9%GrHч9d?-mO=&5!QbRR+Dсr~lڨd")#?45^QN/;}(s06m,GknjDUTܪ(UPӍ*zF+| Ky嶥rZl-ܔmWsv,S/aޣqLX*4d4^~'^ˊ!_yhq yjׯ?;w9u0S8#Vh XcBѺkcW\zr9lhI6oZ&6 ʣn*N,;N=7yp}Ȋ9%`W()k. 1cg \=BɡAU[9Vw)GU/*20jEKxAՏ)7 t$WM N8 Q([S\zĂۑw:eQmV۪Ӳs\~XE3IN{^=X:/n<'̨Ĩs_4)eYi6:eW&=ۻ-2(s fHWXJ6\B^?\ {@žjO.YI;CLPǴzJ~ ;Uv3>tp4cj69Bgg؄6^tr8X#@fvfԫ~17.r0N'AU)ja42lfVh(NMü/nn'&RQ(jXfsLpbaF>kdPvq=9WnsVTEWp(Ozl7ra13BNrj>UHiojH{;p{?rwtN£D?"""dU,9Q1x E(Mk5N|aΣ)U%;w WQ) yV?S4awV}û_zo'>j+ I^u, G9s[ bC6_9N'{VPcɧMG:kGaw1|@75Eڔ_\ IU[)DV;c5f`+Xs'iv(ͣ4 AW鷌NUL[3z* 2V+c,.-@IBP˹?66:v+)W۝(( =Ntʣee_^ ^Qw8ף0 Pj]J8J2]ZֺqaphEޭQדe |yb['s<.fa)\J4 ?ۚ鉩jA>)t:B0סpRIm‰ CB>*,0[g#ʨ96KYZߛkT[|xb~QW@zJhb4$%*,A ِs*<]C&I'DقV t7p[OȒyidwWzdj̴ۙ5l)}}mϟ>m<h>cG"Z.AvfH4T]MZR"A7 [V+-ᴉPWAFf#GgӘw6m U^  z[}HOdEE~#޴vQP@38+*QPW1PoZ#u]~PANrD)lw+û6'>:3J cg}.=bT[ ?Ђ*\٧%q x-RX.{wbeFI&D'wbս9(y tj Ƅ`1MwݚDҕW$X%ϴܷ;IF cwiBzcd$$cؑ4-]ϊ,;'eѺ`)&MHI G1ؔ&M6iwO<X9N}C-i܄{2I̶^Eaz^[rGiʑ>ZUoa!)._A}^VgFb0k\lkPN;hպ'WTJI vHSuNo-EQr(Q^PA٬ZI1S˯b51(pEǟq ȹZZ6zȠ\!sdJ:@[.Uݥl~ 1׏~"2> fJ6fКygڰgvNZ U3sݵ\@t\6W`'鬸z&={9Ѿ`|ʕ^\WY~8n[ ,X>WÛiaa\e1wtiKhv;/v'VT0$<]RS jq"BՊ ';%"<VKߔjz'+%&GDq!Ls:y`qq~뙑O_VNXJ8ESNxa] \;f$eg:Rd IDAT^ҷLPYQ9?ЧO/޽{?>C /; 3`jړ^S 'I_7*p.K2ea.?oc,11l1j4 ;*|kS5JN[Be0nOH/-0FL?WT`ZlXmCU~uBw zQ,/FE1>!e!!شZZ-L&c9]lUl63$)"#b4Rh)4XQ? zɒ%KسgFw.rسgfbȐ8 mZS2]͸=`prA<-@')DZ>lfN-Fν^v;7-aa\bsx i<5ZUϫF?չbjRYI:Imix8kqo܏"lCNBp0H.]dvBZ- =SYk 9Qf[;ݶc?^UsN~}DfX̄dǟi0h9sOE ǭS9tr:phD`HpFǷ,ÈR뱱FAqiچ$Y@Q(F !W Y_@AwQ,݌Zp&y-oXxhj~Λg!gq=)~ʏ[ιx=os8|.yMv;vĥA8D, ԋ? ~[cVgf߹2ka`hvz-Z,U;h,lꮇG>z |koccĕ2l8q y(}*ZD͑tIhwWANERD~ (?EE'?.N'ϗ%R|"NQx~7~2)arL\2 …q,>yYԛy.-7.+DD'S5O݋+#CA 5s&*-fN|Zb)`n@FޖWf@DEi-wF&%8)yXmC5ڿCOgB_Fr0j0%`UC*J2W1_9&{E={C8Іf8FMy-U'P](`` ATrdjr4&^z# iK%f9SIC؅|8{cJfhe.źTN GZZyķs:YE߂H=CZU PU8w3 }ZX_׋;|usfj_QORS_JZBnBB΋| ѭFA8ʑ *hKhYn@O\h*Doyy!*`Lc20}2K,uaiBIO಺=od'sЇ~ 48 5g'ߒ匡ˈv;}7l6-\ϼJST+5 nh)ѽ;^d/fv}Rmlo*f'#Usֳ]m˂*blp1vEZ/*66ql_F-~ƿךLh 'U>ciIx8ϗY^:ӹoVy@ﮨ8ișϷ9QÅ5\5h?ܽi~rNe]q:SY@բ70ļ[ryt/}h,%娺Z6| :l~;PCoeg@gt|Esf6r8baLײ0,hr gndϱn7.4m#zWEW ?olc?0*+5De۰0-+3pd5&f[I GE52L\26J2dƍV$$ӫW/trFǮ#ꐄqS9Sbٱ{)}]*a$čb,mi²2(}0&bK%UVZb/|#>V0ZBC]CYdĬTߖ6*s~$BWM%X89zȶv픋\N5sf&P}e-(ࣈ~ `C;+++.Dy?ZVMaVh(%Z-)n7WX|8ɲvs!.]Jdd$mȮA(#E/nGaG:`ϡ|\"a'M;=vz"XIAA23veVP+{VٗN,nw[NX]R}9맯3 (6clj?'괦3y{k7t$CRVwc Cqf-MY|:AR;.2ֆtK@r v8JZ)JBCU2lDkJMmZpĿ2Bv2펖"J+Mf~6Qq M$RrDN9:_s6fE7_S>yHG/{i)` z=ݼOT>.*jP}cVXgUx8~=XzL\g`>ƣw+ӱ;&iЫg˖֝t^NFQn UӟۭmwXt$SXԝS n daVչe'o4]}\4~ mėa#)osUfǼz : HA8-P w[-?n"ǘm fdQQا]X h- '19:>7v*[zG'8;:6Eca'=[ع`?}H7تveBu]U[t~ML6w[dzE ^'"vΟfW^x(BAAIII`:iӆv=yG jEa͚ٽ:f;Чo_\ruC8"(X hTNːae؋]L,c5#yC[b^uO@q$/;y}kce#_m.{ kHBJvTo* *['J\sCgk׬?M `O?v/'pRIh?|-M.J2sp9ږtځ(uSP  &@lwV~^>`//MDZ_ZؖQrEG՘iٷ3kRg|OhHW Mx]$)PYdA 3ruy\^_1nrRNic&8E]\݀ƘD[admM1S+U4Z 걗2TR9!}Y6%.tsURBCR|%*1<}} vPo\e||<zWF[\P:c84xwaHuQRf3ӹibĕ, rg x k 4Նc]3$NESk.ibׂݔeyEZ,d[i1~"qpAU-D)D#5FODJKg{/gĉ :3y֟{0rDDD*Jf3#F$=}/A ¹HͧRA־̃Ҷb=cOa ǾGqo"i0)d=7בj#szJǀ櫘$ܻ,"XGLt&EX JX_C.oGŁ;хvH*qzBJ~I?HzRʘmLw);ogmcs?2$]| XN]wky;K%NNn r_qU6-?19tb$MKcEus˽nvbII TNW2PϷ8ƑI4JK#` E16K>%` VvDZ,IuVղ_gĬ02j#>˭VXpWUHczf3_g002)k++e$ZQVXސU ;$D( Z-Z-یFVLd6SieZXwUs(""rEDԈ1< l tmnM౨Ym2mXkxa͡N[TTT #;`ym1[o%8.jٝ!o`&CDjUOFhWOP|geCy9r$هl0&OL\|B]l9*XV+a(w9 bPQ^~ktl,n|ζr7b1w/IJ"hOn-N\<.65 +ƭ7֙kZ\AXtp'Z4Rc G` a_T@q`oc(_ la' rlZDOf?gmRWS|sZ49?ŏc{D/E,:UV{"} J_=Rۊ?%u87C"s|MX.jgRqG*#~? <Ԅ2|x$X˲'u;TTSs7U %ɲlfQ͛5w6o&K@ VLEp>0&Sէߛ&8o[h~93ujC^$܂A=Uk<.dŻ8`ѓVs\ڟEܟ_aǿnsEB7-O߮iM&psEQPk5ɼ<\>_{R}qryWc'S T ^[pn(h8u#z(GQ  p6qXmt|]M$S9`ʃ̺ [fsVQb۟7Z({yo͢~n܇ho`Kk;UqϔﹹP?wJ඗[dc.8Unݲ9zZm{O?y(GAAAsߧP5j!w6kӞ-HXD4蚥2 =]W4D%fN}Os lw%(O__Oqk>jΔ7'10ʶY \725`sܝX`g*fG}Yv .j6Lvܥ &QnJ(GAAl( V( `5|WWUKBx' |o[9Hq5|vU|F$X;N=-q`qY1_h{# yɷD.^:p_PK|PS\Ʒ_88?7)hZO]NA]:L&T#rDAAs]<+V_f(IN+gdC:3 |7ϔ-:44__eQ\h9яfЊAAAA  (8s|Ǐc8NpA:  (?š+q9=-N'[7oRih2AZ#  BE=ضe+/%!)#LnQZR(Ǭjq2g˅V+ߜἿAAƀFs׋KH$4E=(++m/Z- dggE222HLLvBƍX餏v;}`AB|/hԹ+ 4i֜F Փ%KKjjlzq\߿ٳg3z(x7 gTA5a*OqCe%x32F#9ٖK AA< Ĥ$*+8F#]t!"" 擟_ gVKbbGCTV!A8O"**r(9(FADX,Lrh,4 M6t&q/t8D1"QT'JQ).cVUV2"9Y" )`0fJ ,2۷o`04*9vv '(G-Xp NI  .v9+/#ISZj EQ8r r PQTcFh(|&i,)! pZiժ+W,ǹsꌷjҁiٲOAqȖ͛p8 3Ao>{^U?U&'i%[  5Z-zx2bc>|D BujXX*+=AjIHW^t]6&&Ν>-  58qdfdي >dx-4W_Iz衇AA8+n:ҥK$M_*t:9p`?HZ4o![ 5WhlشÇӻw K !Ck*zk3n)Sয়`F8zHH=`x0'Ζ-P\ p5Э[hK`!&MO'A>~vJ_?>Ǝ <_ïŽUs0hr tZw=|>u*z2SzuM"w3AӼys?>iiiU#z=uUêUIMmf EQ((( ))Íi੧AÇ=y`t޽~mt5JK=LςQM w>ڴovysxyӟr/o Xd:x}2|+Wz:<^ -ZL?x|зVZUUQX\.z}Zh?u C'(t:JCPfoa`ض~mx+Fj_woXr% kСCppu_^4?<=z<[i?_œVZZXcM)ן pOĨGARqF:wƍ} ~٢q{Yb.=7'30y7}Rx5XɓKdA EUk6lF~5 7# d6;颰yD }oxܕ䓵?,sa_qKꛞ|݁,Ə;M_J{.|@כ_~1nvL%{ĭzIAAt3 ?D9r&7j=m~o"{LT6oͿ^oڎ}!lY{l+ 5Pݭ&OA!ʑ3I\o[7ߴvKй"ߴN|ӎ \.*U_}9FAbrDA8"3)3 ߴ2?f4?nO+VԞr% (& U~ZǏF ."hdS[ l?ى[11u/?E7^{°aݳAtq 7}f ~ƦMu/>};LY׮Mk));w_uUl{޷RsAAc9bZ[ Pil|oZ~f:T{YY'7gMt?}zz7_]yaHO\޽cAAsDQ1FAiL.Lw;loSnȋq_W߼+8'GmL|cL_~~m {9A l۲ù9lߺF#AND9r>qMop0\۷àA#ef-|^|z˳c{,Oԅ=3fxvq={:k`rO o^j;$1Gs[߆ !!~l>ͽm# (Bj۳#yq9dg}ۣ:NA!Y'Oέ:<ɾ֭ N=  `; ^}fZ|GLZ=_  t+c=S5c”x)AA8f9B=زi#w`42a(n7κbz b9Ҹ^zބ c/Y= ЩSS?^y-MY0&O=nRO pV9*&Sݺ$<< rp8$ Q4VXNxT_p_kkh4#mAjjn~BC}Voyr%! c.AAQQQԚ?D>}q\v|AAG V}sN/KK!)c?FN0qb{gV<` ~٣`hj?eKXRq'.$.O[媙q|t@.'pkY B{̥4t8p>nDYAAzNKBB4o\$s.x'Ȋ%DB‰`7??""<ʁ5kzQUU}G_bF˗`ܸqKֳ73~b * b̞=ѣGڦ 6֓RTTڵwTT6X@*VlP쉨hwc%DcIQcb11Mc}!b%D%V+ HgJ9gq{3w9γ3s_zg,$pϖfxsYҿ1#?bb`DigXJ‚:u`kk{UZf>!XR*-[N^Fy>m^z*F !F#\._GKL (%9P(0 TuBZ5e>Ӗٳ$.B!(}#I<% |B>e ty<u4o&!! SŶH B!D()4PdU((JT?''Zuyj_݋,*"B!#III ٳDDD|J;;;jZ FC$ТX|!A I!t`^ԅtB< !!;vySrT*UΙĵWؿw&Uܤ(}˖7Gu$QDD#ƍ0V Eǎs LJSC)/I`/+L9/YILR/Oz5T߼M}0VVk4mj~dZY!㢮ɟss3?ullOlʖ5rZH6) ={J9ƒB!DN^% !OFsx7UV+sj4UBȡ8::aI^̈́éӧDz:/FuTƣzWзn"g}ۨ_#[LU^ ?!!>Dw{c[FY{udO^׹s0jT^!66H+ jԀ7 •q*0RĘ{"x{C?7Oٳе+ԭ[kDE5KZqq9ǸUWBN)Bd"\*Ud2ŵ֭t\̂E }z*Od޿?Sb$]/(ՙ3#[/sb$#޽N(zYڕ)1{׮?gNdt"+v v킛7 wc֭#c ?3'Frq1##GQI'jxFju-ƚL`4,V(.L&>7[R m[mѶu&[[}}1YYa};KWL⫯Bv5e9 řG/5jԲ%&{{>zɓ=U["t:tP6;c4`yI4&&¦M %OT J ?0OrΜ ^ҥsgprؼ>L?*:e֬_s9b.CN8u'ٕ+wlhSKrD!y%jx1%HWҽ_g{}6wl'6>1&uR PMTh3$G_m۠>59#ԩt#F$zwޞ!>U``J']}$6[ 5ZA+' :͏E2gg:u۩SKg^IX.upp+޿yGG!5^- CTb\ ɴ!dZuYnc͚ӰYl߿ tJexY0F`V J;Ӱ,IL|?1͜hFnS2Th6_@ nYRPK?#o3*Wba0?k8mþ,Y5?i:&;xR(G2$UMdH^uz֧Rt-6oF<#+.jyEW")HGRERP#ʘiMnҿ_nݒ?B! =t6-xe\3`_?S[[YbmedW﷨L 7gz~+uG^pf7Cy;㠪FF 2M>&n;ı> 8s҅\zilS~ޙ/ go_"F|~4jԨ#M6(=e_xE)_wkBx;3W 7ibN#Ҝg?qϦTP2U1n 3݉KC5}$z|waנ0VuLE<(*m[N-͛h@bݽ[z?M9IRPi\ 6վ=XƘ`X¼HkV1B!JQ\-+W|dI{JL}eƑ%/fAVEBZFM͏`rC|KC/ѯY5*ĔSq{VEdv>opkjF}vr]j6yG"!Bs^~o_8Q1 $$͞$$szV1Fg>Prt c[jH9Ў3<[ߛ蕤Ĝ;2>%[doOB1/nCTy1Ҵ \)pXɒ0pyN/˵c\ ɚ#B!DF#Sjx: oW%}&yTIq`zx9}sT-z`EUFֶMYCSޓ#9I974PQVU‘n6Qc[El819NOxH.j'%eu$]~ jW+e.cb0h,x+/x U]ru?dRk`2M Kv^ױɵHϠPfL0a2,'2X9թSYvGo ZOڵL@1B!BV#ċGTJU5 (2I+a嫗QUG`2Iz􄓒%~v,h™NL*Msbԩ*m,~f'/CnUP<5)J:`\} oB}KPDERFg ]2'I"Xhy`k CVI//pwG0OQʼpkŊ/nbO<6?&,?^ɏGpgZ(Sj>ňd2ƌ_Q!"dZ/DO4lD ܷoSl %$؏ٍFɜ!!11kׯ2o̞9jժK5sl/^,OB! HbQ*rw#l֌eˢPl2{.GꡱHWZOcLPX[[kn+S)Q$H($9"BLŢRqqhԉDEEzqUR=uգ\XYYKW\\\HF}vTJ%SNN$BcB!FCe8:UDOT *J4I<ߏVVV JI̜5[!B!x%/ZJۇ/XXXSlTL;&B!ؐFӴmۖɓ'3rHlmm%A"D1#jB!Œ`@VS~}oNdРAHDbD#B!BbIףRh4xzzOΝwAIDbB#B!Bb`0Q(lْ7obii/֒ 5GB!J#VZ?z`Ԩ*xyEUJDB!BK!ejJJI3t`{ai)OE&B!BKj_ix׼{$dн/.V#B!(VܸqV򲰰~\r0J(Y҆ >q* / V#B!(Gܺuү_?̊+xe# tFTjb0%9" JF!B!#AAADDDo˗ٙ>}1zX[kQ(h4*Z JL ċHB!Bb͛tꫯ`̘1++m'FCxfKj.˽pM Fg :WS8<ʷٴq߃=.!B!DԳgO3f SN嫯"**crr4'OGb⣮!wT<,P'6_{v[y*-%Զ~q;nl]ilCU.=귁_y|ZU(B!(n ۷wߥA4lؐ%K`kkG} ..)vLh;},:E]$_,:eӂ&J6 _ fbpTU, !B!R̙3;v,קdɒX[[0|^yΞ=cz)1g݆hFCbDzw8?u0:tr0ּĒ>)A \%mN烮DOFPB?PL&okNCkfL9y%FlN׍]V-8u1~o{^w׍ eӔS&sC YQu=\ԊVyke|t|SFIcH\ǼOgכ9C҉.m!<ދyKy.xL nbl8{)lG1B!B|`ܸqcooRDTҼys~j׮Mδk23 fs=b 'bz7cښ39xOLxZXGN5ʎFfu_G2Ï3A?rLT E X F~ /2lz ܹkggA8&ve9:7Ѿ?[98 M~cD+tw@5j $8yY%u=V6/P=_/vggt3;Ϙd=鹲x$rs[)MB!B+* -jfƌjUʕˣ<٣@ƶGjiϋ*]a,ِ5uJ?ګ>xɻ:n>MqhukwTvF2ۍZM2e+a{r3 !rKG0EEy1mUjғCx~ WTh'_LGH F-WθVϤQLt `dJ$FjThʘRﮃ]RNFX|].eѲs\6ڱa7l%gdy[١BME>LoW%hi KG7y:cW?j(B!(n E۴Z-Z[[ \SjBGUWc&蒋&kYò73} x{f܍N~gW ;4f%*PFkK\W815)ũWhN\!tJg&ֶJVCji-KiҕJv!I8 '&,,~>||,D,oM R"dsk(5h %E-׺{,@TH)*[jַ%T(F3c ^ 2rD!B! dԣ2arJUNEv簮̤]?qltŮQ0 @@L [ PH[gM{+.NLיi#*KF$dZ$PSk>Xfxutzo&*ڟ;{Jrl]cC4kyrݓ{I3e2aH2k~_<B!BǟFA1QP˻dZϚERBs \M! HOt;i qn m Lrh؎z~#|:GubDK3u$+b S\fGǫ[0#?cƆ ݸCh9+.$eqEkVS\/2MTRm_Is#cű1}uèʑ}(37+yP&oeo I!B!ij0bS,X5yoc]7wnfX[1o̫}*2N'`:.Ş~516Fsrq `gtؕۗs$Gokk9{Uv0hQLЏ?!:BK=ofqs\ȱl~3;[Zšb4Y2/u.?zN-Q_ seHt7`lZMF¶ZࠪH&vr!oeW I!B!ijpm2zik*LDV~@a8E=^S_4njږ)] {I9ė\:-*VϋiP;=9 ]7G\FZVX8U]\h?q,d+Ekdl+{4 v{v (tt3vzjXZZQp;i0r-?m, OB!B"-KiҕxHBu;bήviT *JTС" !pbLGAS2ќ#Yf`F{O;xVwlyLgTQu%W;ӎ3D$ֹ9jJ-ҊVI\1JO`wcz+3#*brv)oeO$GB!Yc2bR`ť.Ȫ^+UiΩ 0Ȱ%K `L>ܺ2vDq?>Ƨs^#hk=)dQdKaD/fgOl_ܜ\;忍ES!B!(ZT*%ëu)] $:{6ԩCA1 )eՇ$8&M6D):2t|ik2GM$ʐy_c\2T갆akF^bɀ8uvl ~4%"tEg~(L6U^CAW' {m#(gB4CA8C4?']coCbe "FOq˱Ώxm>ئECʭLB!B<~z.| ycФn,{jw56;G9l(;/F~5[Zqr>gb O ]X F~ /2lz ܹkggA8&ve9'sxɻ:n>Mq#|QCx33iX٥Gph~X׹WIHce|2=UN;ur[(E< IDATtt3vi*25GB!W2 APVdJ5m6+԰ܙM%gU48$q-Ɣr87v:]),JӨգdt(b5@B%+S:%WhN\,kx3?#'xn6 p)6#M27qn7LgTQu%W;ӎ3Drscs*Z>c8vEJK$GB!(Ӎ)*(0aTs(2&L&2lɒBa~cc֕'=4> E3+.-wyE[Y{\(2͵_[rUbHe*?2F!B!mj8B؟R)]9^`2moW\WSRIulS]RibRvofRچZ7]G&~֜՝\cߑNxrC=GlԹ%eޖkQX硝D!B!㧮،w[ǰrwl<}g{;1) t:sEġ/r9y1_F.7zjeiRǹc9c29V-ǁwu+&0T-[ǎ1cn!QV}-So*Z:2`Ps)K1b9WFB7.]8W4GAc"I!B!ώ1_6HW;ѩѪҁzZ/_\W޸]?1E+w ӼR<[uf^ԳkDFYounT?M󲍹JK5@Vϼ0^_615Ԡ{|֍?vW|>e\K,&׳ȵ,Y(m9H|eLzUZFDa*1YB!B۹-Zد M1䗼[Y%)YB!0u|[R=Dz1` *JbDI#B!Boy;jn]*5!&I!B!P^*# !B!XB!B!5I!B!4[z}VwYՕ oh1,;ҟ;112@e-DB!BG}m}|(46o̾0Szf]j[?v 4gCZn^y6ڳ+MzF:O,bSutoAdB!BŠʜ>flWE Ǿ׹k\WvXڦI< zjC*ahE z7X [ETFDc 8HI!B! Iz?l]8SVeҞ J1*%}SNi !Ow0=dSn]NTPl$O]Ņɭm$.&)B4C5qj2/>Pe4?']coz؋ҕNi>Ao3P83]>K,0gtnm284?fti$/GRwOc؝b)~icbud!y*).MSS%1fKz \:M6婹3ߛ^Sr'y=·fh4'_cs+uX0$GB!y*QhNL>u9l(;/F~5[f*lysf .eۏOnXμA?1!i{2>d)Ĕ-8+VtdX+db-}&r,p)k 喖ckLc׵ gtؕۗs$Gokk9Zov9#prn)6Jn=WDWn{3ƵKjϯl"vv1 wӘr F|BXsz3+ٻ3mL0t9ޱtH^i5a[ BXFk"ء2OCxڍbQ?^dzHs98_]ςB!B<_g^=wZ8Z >F1)Pph~X׹WS9ė\:-*VϋiP;Ly߁meFaAnN&J,h錝*5.Φ_0EsS^z<Ϯ xFTөQ.$X5)]H}v j_ ZwUFkL̟s1./5ZK ,--PbcϝδVrg^mq{P׽6}Ja7i$ݨմ-SF,71B!BM[w84G9ε_ϟsGtɽ_-N6I\1a9߿/nHb:3AθP4 n<$ [QQu%W;ӎ3D$ !pbL딄;^e9q)CcsU mJrVI\1JO`wze[Gnԏ;w{[S 5~*@MGr6Y$R۷&ө4]0=Vrg2GUD 8Sؑ0jRYQ5GǨ(QA"6|W#B!B QٔOt`|d,_9*ުY S!% gzYc3ꭜB_Wգ2|^z&#&U V\Z@YgHxe~j3;%ϢQ:Z$3:^RZOy l _ӕ틛S^xtPVfĿ2F!B!Kq$I_OJ ; l\*Ѡ* TA}GJh,qf~ߠ*a⻺Skl;2ShRF_A7 zب-qALxBu7 1p۴5p5+ u>IҭgC" <_/ͿcwgPH7~pS":S#1>D u.T:ޒB!Bލ[+GҸ9~Ư :ZOEK'qVEY7R x+_}W{e!˨p25z}7,+fMdVZF9Ф;/(rF54i6VY#9o0Q:taL(mпu/6stXwoYo uk楡s&Vڱ5KATpM~'h0971.u@ʕFG+WOUat)L0G!B!a;wvE$B2rD!B!Ś$GB!BQIrD!B!Ś$GB!BQIrD!B!Ś$GB!qK u->,JÅ7Izq_택]OʝE U'u͟D{<3D#B!BQnaFTA ڼ3LyYRKtmf,Ҝ ky}XkϮ4Y|u<3'P {i3D-„B!2絏YyUQou6NmmJv~} VoD6C#8!B!($}iw`NZiJ{&/*sƨ N9:1ZDH?MvE8{gSACw8?u0:tr0&|Mi1c9_e&aϺ {ZOr=t .zƲ=u0?']coz؋ҕNi>Ao3Px|1Kf[>]<3 ۴mGٶ9y~No^Sy ^*yZvi9؟Ѧ<ŢPc}5G<\%9"B!EMU4ќulGx X# eEO#f K-o̞%lb [ ˙W[U<'&d<-tE7Ѿ?[98 Mb^'wl\tiE`aWfn_ΙÓL\j9مv\ɹ<*\Y<]9.@vænm!,l aa obPeV&_4Ӂq~;˒8*oxڍbQ?^dzHs98_]ςB!B-E)Oޒq)ݛF?a;f"ivhdC׬F)h ϋT]bwiJ:Eڮ~a1d|=Z%-@E'xVX8U]\h?q,d+:Njz.qŴU}KO@ƶGjiNKG7y:c!"ϩ=PZZ`ii"xG}t5Pe*K_Tu)Kט?b\^9^hukwTvF2ۍZM2e+a{r3I#B!BQ,>h*毟ѷ W~EGdԁ&{[l1a97uBl Sy΍lGؠޝeU7D`,2h$Y-ܗThBMˍ,ԮQ˟5niVM(0 f+}rN9~އu}|ιnRrǝhQ6v3pLU2cxӸ'iqKRR9fm= !?ξ l[LJur89Mk]-qq1]p b TW;}/^hB6*w3rx틓o lǕ{RGj_>2/iJH0ш<і8we'&Lfsčb[T^Y33*BD3|9NJTLtb4 fe߂002RX>~9q1 ؆8Q,q#{w2lϠ07a'u; {Q7HD 4!;;QP3o5:%4r18WF%˷*doLDDDDDdiE0Xv#6OOLlJ8SJ!1h\kb5wA>D%8gW0/ b*1XWJ؈3|jo^ _>a4CzLEOpJ!ƈF~8߆GLŅOtI7Bݡ뵓z?pVZK?mnmIݛ"- ޟG}<ĊwPnɎ*a^wf2 _$"""""RlЖ͇I'÷\ͬ$\q#Nd7cM5F#""""""0{hۉ?E5gFw9"""""""2DRFDDDDDDD45GDDDDDDDLSsDDDDDD^yi8}Wշ+}5~e#""""""1.i鲇F=QtqL~av4uR͔RɓchJ1/J """""RD9)6xFN󊠼0,?EN0m\3VkVxZ~N7JɔvKky9(?Nն%~0D5rXak?냏5kp.;IqR45UjvMG3y3c|+ !cCcnEwS5bOjӋ@g6wN5d]We0m5]F\Gűi _ jx8Ҍ)=_HD$^q6%AFNKTLGrNN"v}ΘAnq;,%.+>315_Ć+E6fWe<ӞƲv򎮧O.%.cxq190r+rQb]txS73/N쾜"ăv+13!Sr.c-)#inX}%]w#"""""":&q/M.,,I|u0tٺ5Tcmۈr:l7P?.=_}j mLG[=4-{9 y)ۙB9~sdkz[ˠ,hB5 _:RXnFEwsqQ1rj\jN˓TEI#~rDDDDDDL̨ nO``e>BZNIt5LɄMG}dNf\+ ;6}AFQ1?a>߈!l:MqũOqQr2گۛ$5ҡ戈]95r#jUplqzQ~\kd13"'rԫe-&/JŰN$r8+*B ryn;N붌fWh۷|svss\#G9̹+M/G9Ą V ˜}Q)5ߤAM:8Q,q#{w2lϠ07,TbŻJ(}G7cfo"?o#uBP˧Y>tY4}g"C\}ysز*x6|2>j>nX/'.X8#;}n j[<~/FMJH;1?52+UWO%ll>x@^Q1h\kbaTN=r`}W&CkpDDDDDDnaC[7BL6uEh~^e9)6yg _Tzݕ} ~l}O܂ A[jm.jƨ%}[ijH戈Ƚ+;p ro3{ە>eߑY9u˷F{2sv'PCA/nґƳN^k.MBV{ɓchJQ%)ˏt4}`:3}W:ujjۏؒf?`g}F` e'ɾ.LaoʧѯQbVkJk |"J`N5ɗl<Ȉi6Z=t<gMG3y3XXJnkHT/7;5UG1Su+/^o2.mٖQȘ/ױф3S?"(?!/Qok<.t8[Q#""""""Q>>ٿ5bWTCٝ0:VcZ]gB,\Z2Rbm#iȘH`##U4ՇφfYNaK{|V/n96=1VocTv1Cէvxy%CH۲"/e;6}ְU<3/Co+]ݛT:c0?a/1lg: ] wCÉY؃I_2›#NqR'["ϛ&f{kEZ>fe߂002RX>~9qF?+n51ӿ8y2 !כJw2ja< G>($,V$o#!ig+`PvRWC疗]fHzdg琝C(ޖV?@DDDDD>`ˢq;ŖUSyZxes5ZeW)M䇓m{ĴY\j4Kt# c:-*}oU:_+?E7 '|Dq`[@]ύaXm ʪ8ͧOP;VZts <؊Gm\N=r`u_ڮ$.ػ7SEl[HK[e22 H!6lhKU2BՈH戈ijH戈ijH戈$Ww綟O›lJwࠜbNNQbv2r"+O_^.،|ڦ-!{D*[#鲕#SN̛^AedD_{݁{Zyn?ú[bɸ'~jDX|QtFN󊠼0,?Pӿ*]y\莱q,&S6۹/؄&jE{%nȦ O;gsXj6|r'3֒1:4kUaK91 sl3xq< kKx1L\[@\3wu15t}f)qYH`##U4ՇφfYNaKR)("-mii+_osm{bXnq6%AFNKTLGrNNrp周9""""""r!jc*ݺNfVKQ_v2y+C{@U'bFdz`9J^>L_Uftc`тå6KuDudOWνz g 5:=K+K_p i@<͘-)y#&dXbc+$y R>Es`Σ-DږЋV*=.&"""""rea9Y|ROZqca]gC {*ƓzӾZqڢ*^W$oqORIJ.$rϰzWJΕC(=׷WTpH89Mwc\L`bxp2 e8[`#?<O &f{9򠒝ԕЪ'ɺyA%mfgAvv9 ,Ķ87=3a} k#v \W汎ZmXhӍf%{`+rov*Oow.F1{ %\Vk'vXĚh¿ߩ~^:dfb;5]s^l5m?bKqY|*0o؈ȷ8vv?M!Yqx`Sؒ~zZƁVCؼ+؆9.Q <2u*A#;9)"FF3h?ѫ{xiMc 9#""""""#CX wL8޳c[WɌj)9J+^^J97-81l;j1Eژ`V.>އQJڌbl,ZpG ,t3 'fa&}ɜPD=ɅXp i@<̑g؀!uST;_< fh`<ڒAmˉ\+]ݛT:c0?a/4',*ȽgK=yjō-u eHiXRrǝzu{p}jm^8ǿ\.G3LU'{7xJݿT*{)%c&rBCI(,u|:Z'w*r$Ŧ;_F~x\郧KyQfVLȴAx3Sٸ#So:Sb.$';w_5eaw.Sqķ#c?;"tA4=7U|L00(ɻ L7'h( c噊Z Mq3O|*[Nʝs*kůll~k=ny!sfaDŽ XTAX j8_78sX BkzR̳$Ж·In:ن96NG=ʹ8ˮ.G9V8sc&ygNaC[&HacώQL݂BQVn=Xj#*J """""RBΕx<1OFa;֭Q, 3ɛHH=G,áXr  ^ʟOq|v%q^CF-Pr2ߙ{է:i} &>Ha(`D㟃d ɤ~_=yNwhdN}~͉yIO.|Pa5GDDDDDD ևhTL୞ch0ZVL6VFs Up(dc~2?ߝ'\-L g_k!d_j*aȳ_ >a4CY~/HSOb0 KfjR hI7k GBgh1[©%U*w.=ǥuyd9@=IV`YF)J2P!IENDB`PKHolh"&fityk-v1.3.0/_images/mode_add_icon.pngPNG  IHDRשsRGB pHYsttfxtIME 78uPLTE1Nk555D-FFFXXXXt^DcccdT5g^Fnnlnnnxwuyyy{{{}Z"`X_i~6iǯẕv‘‹ōʖ̙HtRNS@fbKGD70GIDAT(m0rй@;TŲV}`m'< wل.a߇mnTaR86 LXnwa9GaA$_W aQF r%IpM,l ,PQFB$120BIP&R>'PҙY_Ewp$IENDB`PKH=%fityk-v1.3.0/_images/mode_bg_icon.pngPNG  IHDRY sRGB pHYsttfxtIME:.yPLTEwє (ntRNS@fbKGDhQ`IDAT 0 0?#1rʟm@:)(!!C!J!@!-(耂 =Ͻhg{h}/ .?w +rIENDB`PKHe  "fityk-v1.3.0/_images/lock_icon.pngPNG  IHDR sRGBbKGD pHYsttfxtIME IDAT(ύnQT:!Մ%] +%GpE3 tј Mù.hm]xv_r~>Vjv:60L}1 |ߧj> fIۥT*rznG;ޡX,Xur[-u⦴})=&S%D'8 iy,6TvZy}PJ"bߒ1"RR$#!;8 "%k&HJ$ӏQƶAE|<&PqMAȁ=7Lt- )`5HXFkVIENDB`PKH%tii!fityk-v1.3.0/_images/fit_icon.pngPNG  IHDRשsRGB pHYsttfxtIME 7)PLTE ! ('$640@?:@?;BA#|!` &OFrBI ;C!! 5oM*[ 'Gndo"qp&1hl $y W}AK!"B-ʒdw}ce$`30F,zڋPNE)'T~zi_1kIuIcw$LĜy9fLMi+J, +<q7w-G2 Iˆ 1<jXz|*[.j-RBEu^\鮹1*yd !&-eIMe q Rم-〥¶VZZ-m՘7) X*-<J5,J, qq]u5IR!k۴E0mݲ>AKf^ rcxe/k=(gD^`\IENDB`PKH=y**%fityk-v1.3.0/_images/eq_fwhm_icon.pngPNG  IHDRY sRGB pHYsttfxtIME".gPLTE1Nk1Nk1NkqtRNSLfЂbKGDf |drIDATc````.@8 )hAY8R U u 06H L (G"HIIL)PN D0 1A`yV DƎFG"V+?PIENDB`PKH`\.GG&fityk-v1.3.0/_images/eq_shape_icon.pngPNG  IHDRשsRGB pHYsttfxtIME76'PLTE *+ ,#1&4!4H"6J#.:%1=%2>)AY)BZ*7E1Nk9AJCYpHS^Ld~WsXtl{{B'tRNS +-.;<=@ABGILSk~bKGD%IDAT(ϭ1 MmZա0K-B,9$$G"͗-Ux/\5DFl"W "x"V|/c;IU:] es\O:)U8sB |y1eYD)_YaѤ6"VRj> kY?gxDiDͶiZs\ Ѭ`IENDB`PKH& U&fityk-v1.3.0/_images/strip_bg_icon.pngPNG  IHDRשsRGB pHYsttfxtIME8cPLTE %( L wNNNTQM}zrg¼ĻĸŹŻє tRNS@fbKGD k=IDAT(ϭG1 C% z/,2a^G47J$aI,3-H:JA:q'sND >opv:d<$Haqlt!ۄ0!ILW{#\lOdK # JhތIENDB`PKHIГ(fityk-v1.3.0/_images/mode_range_icon.pngPNG  IHDRY sRGB pHYsttfxtIME 9 ʍ9PLTEwNNNgggtRNS@fbKGDhQeIDAT} 0Df.Pw돶;h/<=)6@|1A28AC pJ]7pYVDrXϴڸ{3pot)YIENDB`PKHH3qq&fityk-v1.3.0/_images/add_peak_icon.pngPNG  IHDRשsRGB pHYsttfxtIME 7#amEPLTE1NkXtˤԫ&1=QTXfuw؋ܗޜ9ntRNS@fbKGD|ѨIDAT(}I0 C-[0CK*$ޗ-'ߧ r8>9|tkYqK)xG3 WgpdAU^X*9"HR)2"B5H&~gT`l@K63[|UۑJM+$גΎ#IENDB`PKHRvkk&fityk-v1.3.0/_images/copyfunc_icon.pngPNG  IHDRשsRGB pHYsttfxtIMETPLTE  '''(((222333444>>>@@@IIIJJJSSSTTTXRtYSv[[[^^^b[yr}}}}w~~~ɯ .tRNS@fbKGDQi|*IDAT(υZQ 4 V3DPAGb :#caf\8y2&c)`^zÌgKtZo0dLTCD~_xP[`~'N&gAaFbvw ֖R颇1'V_߇DG>k1F"B+ dpW6y]ƿёڝ"՚óNwB!ؤ25(7>oY"{%IENDB`PKHN'fityk-v1.3.0/_images/load_data_icon.pngPNG  IHDRw=sRGBbKGD pHYsttfxtIME 5o:>AIDATHk\U~d&dlb"-Xŝ`ƕ( /pхQ?Ѕ-B+UVMӚ0!jt$343w|T';_8s=<>s΁]'"|N7WqGloItK~{sދ'3Q O~rl# n@}]pnZ#0 (V'<y wBd< :C&Rnw":s0 6!PRaC6Jiĩgp%r}ұ `L+ ~ , yP) [&0~% LJ< ƭ &$Pr`jk)_ ߈5T &0*?v7宋 ( :`)A܀0/j(6g\DM,ĊhVBL$ꃞDkjjmJ_JSkP=I->$:ݹ_`7f3-@baf)L4[scX57O1ZfTX ~t0"Z:Cbt͖X-ĔTbaSqx-6uBa;nS4lA4K M9jp›Ygg[՚&[RN2[q|Qˆ셾g+tl8}3~4>D/5@٢yv{Xl.Mzf],/ȁ {XX5 cID-?\s~m>| S^ׂWF5hz&sӄb@6=1) `;I،:y l/l`8 J;7 y] | j-WA7*lQ2RZ@3G!,k >ۣoA؄282jHP Ԭn!"\9q4ldJ l(~Ek6[j&"wuI4G ~94<?>k-^aN%E{|,gl)f] ~RpGIzF "lio|]~sVrv?v2X{EMX [%_ҥݟT/7a6rD" /w 3ŀNBt6S!dC趑PBYv!626lfodoTI]|yPXƚilπidbȫeJlJSDGzٔ&hljb˘Rnd ~ٌ%Y$b:RlB!bލ 11NH})P9 #!X$lbkkjC$QHXnd!:0ۇˇHuԠݛHR++b(okj{jDdK{+ \; &kl]Ö"lJϸ6}p 1=([ĽKD3lije? Qq )T.y,#v`yl<)%J,]]6-fX/ljѻo6W]gH*pBiӜ !p6f tEXtDepth'# IENDB`PK3H;3Ffityk-v1.3.0/_images/math/52b1bfbee8b72201dd7b677f47d1e44c19f76083.pngPNG  IHDR ? 0PLTE/ݠ{tRNSt`H2O pHYs+^IDATc```2b 0eMp', X ĭ~ZfEsNν3ˢ&]jU Twg<~f@>5aS< svRsG|1?0qROM*d,zkm.@N&^r[uxנvDZSq+cUJC2{ `!npr/E>P9~|^y'q?| u@ԩM^gݳInP)z6Zl;*$:ӞLXxX2 .j[dZ ?7:a׉[& 6Wdvm,r}!Ux ?}wgxR&ί yAx&O[`BMw2F0,iImhۙE X";tIKpĒ &R֐ҡ_`<: L"A(bV8|D(*xX"7xX" &xX")x?n *!Ddq~ *);]'(@$0#%SYc#hV0~Battt`&s2WC-̼ājH5ĶVebCäd0 PR\zB' kÜls b7"V $.aN4u"+yV[<8b$ƕ ix$*\ v g>g:{@NFU0d#?Β)VKɻRR_I%'0DVX!Wl_d o,iXeZ 76"EdrbH8Uf"Ό$h-tH$aC0.FQZŊlȭ0jFGhF~j*55!q_R*&,UZ кgH@_αL tEXtDepthP$њIENDB`PK7HیFfityk-v1.3.0/_images/math/0b5ab196db6349cacd88ce55321b1830f9ebc304.pngPNG  IHDR3C3PLTEL tRNSHt`2{ pHYs+@IDATxٖ*E!Q `'\yU%3KQOÈ슙U Q>;v d!dg7:s3f#qaA0 kMiP-I&^عg"SrJX qa.XY`f݁=dfjk? *CV#iEfJku1[:5e`+KfA)y0 1\HMbW*Z`_Zc˪ fʸ~Q<0iC6B9 𪓆r600*3'`Uw077p `UIz L 9 \I lXcf``o`Ѕ : o@+63`ڻdiؗUlE | ]FTD$%4a%jyˁ `0#Q̠i2!`6R:p&>7/:y|@[ tEXtDepthWIIENDB`PK6H~ +Ffityk-v1.3.0/_images/math/661356660ddc5c3f606be9296d3c031df7eeceb4.pngPNG  IHDR 2g0PLTE/ݠ{tRNSHt`2z. pHYs+0IDATXX]U>g2I"JcUڀ5B[YXlobvۈTѦRB Xtq("($-?}ˢ͋HmVY>mZ;Ln̝ܙo=ιª`lY3Ϯ b&p s4'8\Z[lWIv3GT x8q{ z`+ojXwnGƺ$>Ę8 8ަo g\q:÷87wBw@HFyjhoP5GUzЩ3wZyӈc͌ѽTI9-P/ _%[gUTn/ U?:~q3}t\=d?kc<~(6K0 *I: ӍٗX rb&)Q<"Γ}1uD, ß'!Qqir% ;`1sI6ߝ(!d NǯA 0e_?bI)iߓ(n 6ckug}(pI,>v8K <ݬ])g~$~>+ YRq t׺H|20WF/q`SeXd36 ndL:mLzvwjD49;F?t2F<\Nſ6,c%M|K8Z  &>,#%mfV4-o0So4L*+0?~` B̪56 ::Il+* IA-n4;C _b gbX,֠)@@MFek-27+^E'Tդ4(B| 6KC <ˤJ֐lOA:j8D3Ñp FL_u&ADm;lbأҺ@hqn4ym})ju+0~Rx\IZiA@ <~`(k\ pR#+zK+礣I!HD1eTmG]asI's|e +q`"l'pcVRů}K|Fy_yk{8: %p{-Vzb?1X]}tib]88e]AewT6,B8'6gQdR_ݥN4<z@c'yBg)ߎvbѕ)mOX&O[\oaY!jo1>GЭkaaUUCz'nP0 tEXtDepth'# IENDB`PK7H{  Ffityk-v1.3.0/_images/math/038845d9aa07cc91c645af2e46b7958765ddf3e1.pngPNG  IHDR6V0PLTE/ݠ{tRNStH2` pHYs+OIDAThKlE?;&>*9P)DKU *p(F\"UᕚG`.ĤBU_z@-CUD#UN$)mefwv3`x|VC+P awkB >֪1 áhvZ5^Pt3V Ά,$7_?+cOVLkLq'J_E0Iۚ(ѼI/n=u[{8!~+"ȃz>HrnI!Z?ӢX].]E}y4؈\h$׀ʸdhABV| Op/}k8UU+U-"YA¬N6ŝRk*M=j IfD68xY\ܩu}6ŝ XZ&Woi+mUBu[ ڠPԛ? A6ǝʩih BW8 b2jB;>vr¬Ud;Oӝ&e4Ѯy0R 7sRt^Nn7疻dJyӟ,=ةu{e<I%j sAMY֭iMѦ-4)v޹̒\O uz7SLR5 c$zB :K[= Ӕ`Ӹ*3*;N gz[к>H:IB_0)ՆiغCI<'3-9ם=4Z T'H~ ,%Q+`P1h􃂯D4ʂQΝ[:S^v=93 ] yB!n#v$Y39$8P-6999999Y:7ubo)s>~2g;x> OKs;F&t9AO؟Y鵋cA~JVźqTח,YD2'Z!9gyBaRp80v=I$`^ЫPjA&p~sog[y4p7M|a{yܔ6qg9_Ki;,bײghgʂ/*83nމ3e3T1KċG K7Bʤ`q8'Y (u7YhH#9ϩQ6g}5ZLN4 !sjs-ʻ)L_RĜ!0s]L-/o8ӽ,ycpfY: ^O&hS Pܟ)VZXݕeِ M2vJTIf"<"sןg33\.^L'ҹ4)q7Z]6.yU7,C/X$q&+-6RէU1DHzuB+* 5FcK@6v\#i{\h֛z<ύو^5/UM-g G'޸I}Z28ȭo]'L{ tZ7b\]mQ/fr85sPN"UM#Bbe٣֪P{ %S'ŏ҈4ݝD܉/Yep0K951 qiDWMĉ,_W5EsrK_v4KsocԺgłӝĥG4B5?7;0~wWL%Ï< wys \ZgOwT9oA1+G}\aE9a[ɗWQ5<9G\c eSBΩ @}qi4лW 2Tva$Q./Si>>tr;T+;Tv-ڍ%rdx/?WRSHsVXٓp- C>qM)=hلbB tEXtDepthP$њIENDB`PK5H{:Ffityk-v1.3.0/_images/math/2982f903778fcd237a209f4ce59a67472f237e86.pngPNG  IHDR7#_0PLTE/ݠ{tRNSHt`2z. pHYs+IDAThݚ[hU%3d 6kZ֬+M"vԆv-T%$ ԩhK7mؒMɋA(tChf339gΟ6>FQkѸXsZS:7Ss9I:*$ItVv$ ;~k@ u=%&ݮG'3,TLRuWO>JIfv9:yCS(m!vK/ARNfi;3""AސGJ-z{%/;|}9"J'ɟHnEˀga&A%}tᾁyh΅];^X *|YRI1)>DJ· Bx" 6$<F@d &9= ZPB/S5u~PJILy"7r̀ށ4U&KB]Da)o&Io'}FnY/ RI)WRJ*+ 蠻$0`=|@rqee$&p 5W%qtR w ([KҪJ+2Z;r?T'H H%wa妭si5ILI$XLsZxQUTIܱ$$c7ܤ{bAWF%8{Arr $l.PI15/ܢa?M3:/WEZK?oT`kg裮=[>|K݃"Z[fb6穤VNw}:xaH0upݧ&*2)#j"w@H uEyG'NǡI|E$Vqu(*H]O0~ft&‹Ht Qw &&*e%;䈺h3_WdTR]LU\O4E%fz2G] <"׶yQ~"ҢԀ~TxGhLjEoFZ* @OZ2**hNxgtiy@RCOZ2*]̤%R/s[}bh% :v^$}BfK+u9dq†w9ҢjRoqQnŅکnk|Mb%%m= QRyVliQMnZ_u[IoiQUnb &746&U5k" pHYs+^IDAT8c` 8c 2@+9n;J7g SD\`+@Lb -vH\0`2q(nbo2!! CY!گYzE$gAbi pta`Ѭ`] C \`Vp+@l-V20㱿.& ,(`01`AO ~Y t 7l9jjUH @!CQd10L2?jڷxan~ x".šYI^=y tEXtDepth*IENDB`PK8H7Ffityk-v1.3.0/_images/math/8bdf5aac2c765e5132ecd6cdf51b546dbc541583.pngPNG  IHDR|(/l0PLTE/ݠ{tRNSϿݯ`Ht2LN pHYs+IDATHc` x/2"|, Ϊ@ `7@~^Tk+ uXp tt2MP#Ψa<Ѓ]qv.Αсr$cK܃{|6\ V{M4,AU; DjN O@kZ -U0B^. 'hyq cAA̋Pu!3fHV]nFYe411ñ)10b`Pq x6ac88V$n\~-U͹g(!(*-Y`jg|Vٚ-!< ϑmf]"bI` ~@j{9( kR}c8DBJTVCOE")͹> S0i tEXtDepth'# IENDB`PK4HŪFfityk-v1.3.0/_images/math/3262f202c4f0920b5279a9135b168a48b3168e94.pngPNG  IHDR&0PLTE/ݠ{tRNSݝϋt`2HP pHYs+,IDATXVMhAv&f[W6V1"'VAxQ`D/ dAA\J1DҮ=EBmiKu$JI6 oo{oSҿJmaIϖ\g8kYM9?)hXҷf'xaFp q =n>3.5c4A( E!]eaf'[.h儲]nS̍ϳVvnu>*(h>YVDpLc4@ w.c`d@X,=!jc!:-츔_G0IԇGB3`|MH-Rò$3Ci4F0҄+w 9!5!2]1f@lCXl+\cFWgo{R)y{3 ^u_Fr)q]CL;͓Ksԓ9}Tc +oYR,_'pNT m}?&:  ֝,7\t1#qk z&K.:%dBܒ 2jr_3oR(Q[D߶ VbX9 tEXtDepthyIENDB`PK7HFfityk-v1.3.0/_images/math/2888ce7a171931fa57c7e903fdbfeeab84cb1370.pngPNG  IHDR  0PLTE/ݠ{tRNSt2H`@/y pHYs+YIDATc``2a 40İU]mCg/Pj>c-̧7l9`-u"qs Bh 63 `)_ tEXtDepth*IENDB`PK5H0Ffityk-v1.3.0/_images/math/02b58ca52e39107dfad1e16baf9996d993246743.pngPNG  IHDR@0,0PLTE/ݠ{tRNSHt`2z. pHYs+IDATXKhA$ =xPlh (,bA%(*hBRj=ARb=(>PЃXh)XEAkzS,!-T qgggfzX2˿f7,X$}%SNh/W -QgO:]c}QVIJ26sRqxփݪt+! n ؄d_CY ٦a-̲{;Sz G>GN"G.69_`]:Ed*AHSCq m.]rCz^z{=2;Nɏ؏+RsΰХ=G “m:ll7q vP.>$]I) {&oG*w?ߕLI]] 6Hp^4HoXsgfM>4mAV-H/?5DY_7*I6_ c+~.POG7QW>琢j)vygyjϺ{㸨4WKR;KP}^4slu,B Cҭ2d{Y,aeJSXnȖ⑁m'5հތrbYbȖ0KFSa|ISF%r4@5lLiFfe&' F,ljYWiV6|:Q&e,N*7ۇ`FMdB|LyͅTյBLdE;jcM[.bbn xZQ_=-O8M:#C8$8w:\Pؐ=;OY`փCao:W7)&Ng"<@yXZ蔒"Q<:z:kv H2/3H]d NƯ@>R:Z}}vX_r =;7$@U3FNeg.Ԋ@`~{iiAKRWAs{byNtDoq7X ^(jq3OϜN C6N5q<P!;?5ضӝX3d{X@UDSp++4 @q :OQx~ b aHq[*!q:/tBXuWBF\JD"\<ݟ*bZP3P`8.>dG 4UBӠxJŖzԘSՋ qv{Mj-SOTAϿ._T+ dssPza*3%O(ͼQڳ9IʯIXGGO AϺHٜrPzSܠ8L=iȐݖ4߁C'_3m?%+7$KB^F'\#?%iξ6>.|SF'[o %5;_H>(I&]U+u'B灈l.:;&E'&Q<}D)Ĭɼ֘^-߀i5opB~N3d42əsQB^+ ƛ|>ѕoutA 1չb3ڑj7y,b5)DvCլN3Mn*JY:t. ~N7tiF""7>Al/6ێڀ_5[B!N'"vw~zTaY`(_ >5 .DwC]=APFYXkEmp7Kd4%~q?CiA@]עd ΛYkrK#TP&]Qk|+^L` Cw:i4%#$}7?\KFe7:5kc)?(xXڴqai57LXS-6ҙyAf˛/~ܭvNX&QkvצRgS~(qH{.~$}qR:5&qLZ_ '|yAP|#eZϕZez-,Zޟ@iM(jCkZɠWQ}ş!ә1P]#yPf鏇ΣL8ot:S<Jǡ=7Jc0yԡ 5*{{}}:=S<JӞ׃KS3gZcqXjP&?>V:mto@y8 ֚~(%iOS+LX95kyt:SR_懯=E6= 0|QY"ۖ|yS|m#iOw)|EtPV{ʯmC 1:k.thtPn~ZG~ @~`Afc-tڻ۟ԇqq @(}UӨq~QUi tEXtDepthP$њIENDB`PK3H{盁Ffityk-v1.3.0/_images/math/0b23b73caaa669a28e9f47f32b6216eeefb508cc.pngPNG  IHDRX0PLTE/ݠ{tRNSHϝݿ2`tOy pHYs+IDAT8˝UMhQMMjRRPꢭ/[J`U$A҃ UD؃=E "P=THA ^쩾}?ٷ ~f};;oY?EzZX^9:B-{6)~mYoAMy =Z^~!ʡ,}>{;^}b&{QpЇ0,\Ωӗ' # Hf}x ArGx. Hc;zWbF,ꂚEuEsٚKJ JBց "Tңea }%uƘU:fXψyYPD1nzXY]iQM)Uô(D hFޝsxpJ yh??l/u¹;:C8&t^lg S*^w {ڻB{}=Yֻ"!pǩ0Rwc1T淿E)w|Ks.GQJ?lnRwsq!܍?u6vT[_5ihW tEXtDepth ^YIENDB`PK3HRFfityk-v1.3.0/_images/math/b1f0d5d72728cf446219bc974304b25c15129ab3.pngPNG  IHDRO#VM90PLTE/ݠ{tRNSݝϋt`2HP pHYs+"IDAT(c`|CCC7` >HQ(~CR`|aB`ƏZ00AQbΏz/TO@q3W L"UaxN:/\0ã r`TđKCh PIdc(o41m  LA 'A DZ nY@M U F Xր#|_A|P;Mgcc?4 ge<XDY* tEXtDepthWIIENDB`PK6H:Ffityk-v1.3.0/_images/math/0f8217c3de937a264da9f5ba8318686c9737eabc.pngPNG  IHDR.40PLTE/ݠ{tRNSt`H2󿽁? pHYs+IIDATHŕMhA'vPRP-nV멡SуDAzXVl)B< =EJJ*(jmуxZJBZP1~ٍbb:$/;73@u_T) ^?/S:[Ery/i$蛩Zξ2 ~ab?h\0G:U=izك* LQo[S3lsRͨPαS @ >7HqYLBujCa>D7 c>ꖱd75Vm EM)GnXq2s"h`1`a8V[#/JV=F(}H=sJ-^Z{# )KMqTfF6< 95jr#{a 2>t8Jᬢq B!ђб1(:%^G$[-M?~HǹŠ戀=ʙj&E-Vw%NL)`ڔ.?_+@ќakkh>J7\b;n&ẀepH%N+$MS9?h34\1$iy3To gZϗԃW9e'Mo;k3rŨm nx=W+A3I3)P,agy_- yBJEO1z6eA6ɠU:VIeHpG Q RL 0sz@p[ Z3S.OsLP*uS}TNWtDI⃫+䏾9YV#POC%s-c^n5u W?3SX? (?j3$}p tEXtDepth@D9IENDB`PK5Hr7Ffityk-v1.3.0/_images/math/b29037954c6d55c57424cecc75b29768894afd1c.pngPNG  IHDR?0PLTE/ݠ{tRNSHt`2z. pHYs+IDATXMhAMIv=IC=Y ADtbEfGmJ0S⁾ڎ3 :1MBiWtSPPm :d;ݩ6awM˱u0zhDlf_+ZOz='ڼpD'd)٣ɻò^9:^@X%fEы}Gَboqq5ڐBQFPsAlb"whXB8=U՟|v;hq ]HS+h5uKM.jEH꩚UգViͬ^:LZS=:RE5`|UofJOj֣7z7]`u~I [XIhf1 _ST-[)1tin8N( tEXtDepthP$њIENDB`PK8HjEFfityk-v1.3.0/_images/math/57c598b9054d3369f6280a1e55ee58f9c05234b1.pngPNG  IHDR0"2G$0PLTE/ݠ{tRNSt`2H c pHYs+-IDAT(c`h*d!`K" o V,δ䘀cŀC"D(TLg߽{w w`]M;WJ>a<h4121(H@yA9 " xb,``/ " xe.? żquHh+RZ(07kqd3IC3R2>>)ǼW╧ABcX)'\ LTt tEXtDepth 8IENDB`PK6H\+sFfityk-v1.3.0/_images/math/90709990f459b243f1d3c1dc2256439f8e783f8e.pngPNG  IHDR Ol0PLTE/ݠ{tRNSt2H`@/y pHYs+LIDATc``2a 40PsȘ`߰u"qu`K?A- tEXtDepth N%IENDB`PK3H3`iFfityk-v1.3.0/_images/math/d6e0133ff7c9f0fc1c1e78c65ece0d278b21384b.pngPNG  IHDR 0PLTE/ݠ{tRNS`ϯt2HZ pHYs+AIDATc``TvUc``/`;O +00`fey ' tEXtDepth'# IENDB`PK8HH#l++Ffityk-v1.3.0/_images/math/143d854fb865a4fb5147d3cc5d31426941bdd597.pngPNG  IHDRz-PLTE:tRNStH2`ϩ]K pHYs+tIDATc`C}Ŷ*}b 01  m``e7`H>eU, 3OƊSCxPC a`hb&g tEXtDepth N%IENDB`PK4H?yYFfityk-v1.3.0/_images/math/c63785bdaf96c8a0319ad19bb2b258c425153bc1.pngPNG  IHDR7&0PLTE/ݠ{tRNSHt`2z. pHYs+IDATX_hGǿ{ݽKNM H}h E `Bp/S>AB^Ym@"WmzO[Ň]gv7#;Bfs3]˚G@%bYn(y)(fCjߊ@hu#xj0·e!þj/&M.ls \ر4/\SRa֡`=Zy{&Ys<{3PǦǾgp{7ΪG8O^hP2^w[ؚ++NKu|:ɍ376}tT<mpٲ {Q<99cŕ,6vy䟶=axbrٽ]iٝV^簤c!dJ\pl2Amx-wػc خc-%=6Wc[73!vˈ9==+)Su Sܲooц|A>9>|+);j-Q9r70ʮqE<ͺgqtGU9&@b'" pHYs+YIDATc``2a 40`,PUUa>c-7l9u"qT'0K'b sY$ tEXtDepth*IENDB`PK4H /?Ffityk-v1.3.0/_images/math/8af5886035cc95e80157ce63ce8c4366d6e4d058.pngPNG  IHDR Rio0PLTE/ݠ{tRNSt2H`c>" pHYs+YIDATc``2a 40İ,C@  x7l9u"qC 0c> tEXtDepth*IENDB`PK4H'RFfityk-v1.3.0/_images/math/db8af72c921cf238131d291006ee474a2b2e17bd.pngPNG  IHDRT0PLTE/ݠ{tRNS`Ht2կc pHYs+IDAT8c` EiQ $PSbI 6t0+6bŚ!Ʉe]q v] V 0$Qx@|uM`^rC7# fsp2Tkco-ܺ81%9c?7mʟ$fȂ(gH#Ap֣7tH$ŽS|˔c3{=0*Fh4G D^ ?Sv`8ݬCԀ; hRdH30lءM ˁ9C)E:sFF˜C箂Da$mn`vc7C^s ]H;&MRRgD"} h :|/(!sizNKK;rjك`8P`83 Wja(isK`E;0V&ham4 &E(Ds~}3k V!2sM1u ߴL&;֛= ag+J·)ܵ Bz&CG' ׁa0\`0tN)#!-_ƺtB4Lg?\ph rIқKvS!?:ʂj0mxȽEޗ0+&zBJ? C;sl o#hD0͖<@xo Lʴ!q WK; %?E# B#&;'`D)SlZ'~F!(̪vhÅK; ʒvHJ dOcEtK`RX? &e PqrjokgInAYv s&P.Q]nˏL@/L&0[;K!1XD  FخrjΔ55SN:(`(u~Ma蚫c|+XhNXm`s*Gv, W T11|G )Wfp3-vLU881K.KdS1⩗d=L?0iE.y6*0wݕ%J8f"Zژ0:2a_ʩ, G Svv!Z-c^0k ʔ=>5?Aw.irG+<W\t_`/K{'"5Ņkفa;ʩ3\ Z:΅2)|baign`a*l{YW߄X!?~ˆ?l%x $ZGKbQh|1 w0? * 6J/Kz"gI&0Z3[Zε^?$`r€, 4|~UPsqgU9 י'j0"Xvq@ Ү>Li0t$)B0e/B}ɅA9aJqg`p!.ǢԊXJ #8ar] |1M[yo(t2s\`[D !s&X%!4G#L5vE4 m.W RʴrjYt 0פMյ?{J69RvFΊ,hr!И992/2fb o`&1)$C ԓ?e_O귧.0E 0EoloDK^HA]X9ap`pc*T _kiao_HA]askDG ͏7o?٣$sPb%۹0Ѫc~Q70bģ܆~!-uu=-`ץw6#}QQ{^:ENC9n`xp['\ 0Ѫ}~,$ߐIm>^A=I5`( %DdmJF}~Er:-i>~MZ{եF-_)\`BwVݷ-o?ǯHw&ۘ=Nx.=G6avӎn] Ej?[&rt.̩ \N4`@ƪUKB߸Nv{8QH!܄UDϢwa , tEXtDepthP$њIENDB`PK6HUFfityk-v1.3.0/_images/math/df16deb31a39f5c975de8621f5496cbe2f9e0c0f.pngPNG  IHDR0PLTE/ݠ{tRNS`tϿH2,J pHYs+IDATHVMhQ6ivk ()E//ZiX+-[/<,ԃi{0 z7EAA+"뼗M^ܒlvvf޼Y^\3pso{ʺYWm@)D*W.x>W?[eO#m5%2M_@̑nM{l)hgWOx77 qoJ.A-6Ɖ|a΍zI?&3EMǪړgDp~pwQ2qJMPblS4xO$#F0ɼxY&/ xt!r-r<1ajl#]_-8ƑV DxA#E8VC|n^OsyH ?b18SOF8*M; f}>`ؑ8}#CRqT:g?gFo$Of:CY7vO_' Nl#dAoԨԦbhZBs ~rq},'x~(3]{&~=T1e$XN93VL1t[NMYIXN0E; aKCXN2D2b'I8-FXח\fD3uhnNrA*VIRf"0I,N"g'ƕGfE4W}pR-fÍh6%܌&9b7ELs"R-VGDthߴZS:h`8Z'S:hh:"^йgDSڡ8[3twEbh:Iν"%V4tz{D4V"(1KhI8SCRg@o_ ig YjփM2Fw tEXtDepth'# IENDB`PK6H<Ffityk-v1.3.0/_images/math/1256d2d6adf6ef0f67fbb11cad708596bbb2e516.pngPNG  IHDR7$0PLTE/ݠ{tRNSHt`2z. pHYs+IDAThZUAw`%Ff"qt7lu YD"ץ䯨$ ;K/T?DTBW,+{y;3޼gg*U;cw?kxY/rrOtiWʼ\b' omMֿ̲PN~Qtr7帤{Yœ /,i֯4}=34OwW2ïMue #0w%.N}=iAE-P^:cP[uPpťG \'R*{nyHʫ9{flYՆ>a , C$264`t4l҇Xp?ڗi@7˴v#,9{Jx (`QF-j4jqmMؚ V 1h=ҡJLxނ*]/FH?Wngv}ay3i?Tv&_qˆ)S ]>1?XEd|‚l?) 7xL%)Ot  p*޽"ߺ!X.xhvY: E9?#bZpgn,{`х j}0Oژu52khхk~jYsa ‚v.$tSlwy9i+Y]7}4٘54ʀ2kKXЎ(j?A{(¹l&/lԘ2k.,XGDҴߺ܏\G6fm j̘),rDܓI _?ޮ ߆YSac _0k",ָSa%C}kIQ̚ N)|YSaau`S!d6wڬ~άʼn;<.2zFt W\h){/A^˥u/C7>ku =e\ tEXtDepth'# IENDB`PK6H/&Ffityk-v1.3.0/_images/math/e0923e2fa1c752337bd1b563184846570beac113.pngPNG  IHDR Rio-PLTE:tRNSt2H`7$b pHYs+ZIDATc``2a 400Q]I+3ac [(^xu"qcgX'WW,Fc@Ϯ tEXtDepth*IENDB`PK4Hu*,Ffityk-v1.3.0/_images/math/d44a0f74bdcf5d93e31fe8a0dd1fa0625156dd63.pngPNG  IHDR VR-PLTE:tRNSt2H`HNos pHYs+LIDATc``2a 4|]$>1`-x70pHm}?) tEXtDepth*IENDB`PK7HE_Ffityk-v1.3.0/_images/math/de15ac785d9c2ad7268d61a4cef43e1bbb58440e.pngPNG  IHDR.Q0PLTE/ݠ{tRNSHt`2z. pHYs+IDATXŘheǿw]59a2\ZC"bNV9d qLmE4J-a&` KTdQBcle8 D %qI@Ws>|mz }5We9?Bz<[O_Qn "0-FG-/'11͜q7}(P^-݄6e-o b BxFkv|oJx{1`4bt_>V0vh, Fjb@iCD͇'y@qW'k̬!"z~M,8W;QIBÏ1!d߄&v9Vbb.;g>ބu-<[ ݞgĄ8BL/eC}dpxSV"~PrEhF/M%U0!|)J͂!XlAV0BV1Lk&DHC(B1*!>*Ĵx yJ!EVCJ!ʄ>^ RQ_p_+]ZZ \CeFRJ 9$ }R!sG2D`6pQwa/Q$ K6dB{D<,9; {/8O#~ǖomMh9ĩ Gw`2EaUE?"[ߔひC08VqwUqu"E냠qI'B|Z Dvŗ:FG p08ƯeM7="v(_3a'5u?ں} tEXtDepthIENDB`PK6HfFfityk-v1.3.0/_images/math/54ed1090d96367d1aabf7d11c87324d0850494db.pngPNG  IHDR7`'/0PLTE/ݠ{tRNSHt`2z. pHYs+6IDAThݚ_hEǿ{^6 mE%iZK<Җb}(6YԔS VRU'"`H_|5V^lJ+-ZufvggoxZ۝_v3;} c;Fom^bn6 f /JK|W ڳKMb9kfL/3dAk`/k 9/I p=tu M =9?{gzXojhMT y4[}M'hUO[Dz*nMzčoi%/f͘zVX$ckw˾X%J7 4ӥ9 w0[#q&-t"#[1M?̉]yN2]hK̽8"mEMߐ{2Rdq8\F~ط'pmH~hN$:Z oXI_'@Ukx4KH\<%af~vƱlVǝ{+Nj;W|n9&$$6;C9Z>7yL$$=2Cb_BcV4ͷ6W;+S#-^Jq_,$Bb{#ԅT_8%cF[I{tG5Ġ़F뵤$Ǒ-(0Ad1ri>'yS-t]Jv6&7!D4$K)aڑ'C氠ZX9Dյ~'fr?ǏQ 5G)0ȡW$ 4Ŷ_x Dk0Q VGgUE"(&*Lno)JH@ҜB5Mu"1"HA P"aBH NTO kå ;P/ԎIJO'^YEsó喲ItiZgѦA12*\Dsp<);PLI y|qF#bSaT\nxXBJ M qÓj1 K 7%+17^_'1X9~7M_ u+[lßŝDRt]IqJ:-Ec@NSLSb7~7$oQf 熿 ߧ;IKv^ʤ8%Fqi膧6 7Hԏ`n/}Z;Cɗ}r&dNĩ덷}}/B>awէ[pm̱j45}h_@`{Y kaaeD'=}zqdt""xĄ&᱄Gҁ$íu= &`{Bdރ] H#7!Z,N^qt#Ɣ7͞r|xN>`ғ)[\o+S?MO}ĝ[WuL j5Ȱ+#-X 6_b122n ^U2nȵ'![`p Vm8DB9W![m|[unQVǵ~z&oSe(N]O(Ls)Q]~6y7:7n:w(a)WּnQVV7[|#Poۨ9ҁ$Bus9sjEɿE`fj;B忰I}{@;qKDq W+ۦWɭ(c#na)+:qrtM.P$8x}?Y:!me=~kk^ 9?ğA)rt@:!-Y==="}UgGG|% plYH`̏츱\u}$68pcan^e s d\Gw苎>PJl{+7IlY,SA\ 9GJ=H`Bp6c>ʍx}p-Đ-8QO *{RYX[vIWnTۯ2e=RPhXV _W=uv^ƨ]Nj4RMxز)wD8mH~:5QINT:Yf3ܴNhfj'P*z $Vga4yZL,~` ,5-X jٴ4ȵx:9TF3]ձ鈸oR#5R q"[ [9h;`r|}K\{vD:;$- iP%JkQ%yM1~G-觽 ]G3?l~7D\ߵG|9:UеOsnL[Q4=lOLu _ę ];Al v383SvᴿK_'t!*L*aN>Nun;z)qv:+CqT陓2>W4h$@픈<>EqT].tM(i{ϒh$KM/rɏvԊ*^sDpW|S/RsQe&z"v$.bHȾռiKth'7e٩t1[ܱ'qA'$I8h 0!EHSNm\Z5&"K!vǤV?q"2r,XN'99eH~t'+?{"}NE*gA~tCaHf b0vD^.j_ԼEU<`RY}ڢG6o;x o,Qf` ^n#EBLU4"!|z Q813V Q5<>\K#OKEX%wg < J X>sɝemAϙ\roRɝoLUxe*ayPk{e*v:]@-Cl|PZ2X>f]Y>O0&` X>kd|Poex@  mt8L tEXtDepth'# IENDB`PK8H[[Ffityk-v1.3.0/_images/math/8156e69d4f3368818ee3af462d23aa79d3b7a170.pngPNG  IHDR3)x0PLTE/ݠ{tRNSϿ2H`tVc pHYs+IDATHǭWAHafםuv I%14!yJY(=O}q{$A.PyzGw4҇cz@! ڢ(bjf?lE zѪ"pMU<m~s%%z–ՖI~2u[QQGa̿ P!gsJ2'b#`DR#.Ʃ v 33/AdpÆDY)jD%f#A #S Q b0]Mknq]htx" z9v2K*Db?@AnM>~q%nMPmC^ :1I7HI}z8n|PCtD=9?>#ҝ؟/ S* yio%J8!;8.$y :O)k* V ZCDzG?DZ4BzI[l(! ^fy'md 8Z !-oBlY ܍yBzrO$hp5{?)qOQ>_;f ussY|СtIo }&@Λ  ϛWF=yM>י tEXtDepthIENDB`PK8H%%Ffityk-v1.3.0/_images/math/91ce46aafba3854e5e4822cb08a8712660646044.pngPNG  IHDR In0PLTE/ݠ{tRNSϿݯ`Ht2LN pHYs+jIDATc``\*lYppe8ƌ)r};0a3'%Z5tڄ@82~ 2,gԗ/XI ܔ0?^ƾy<I4W4MM&+l2k^#7 Byyo 0ݤ?}5SĔJL$;W,< 0 okC_Q\WTyE\3]QWQRf@#"[0CwO/ |;'{:? yA#/;m;ӏn=|$Ih^H`0fATo`>~ҮPLLsX H"u#heR^ AfGaP* U䋐@7Ax\mGPiE>7XԐ؄+Fsu߀Cx&7l/Pѿs$1@#lT{RRWic^`6đ߭`. Ʈ&`BRJ {r)ujQʃQw*VUyBl@7'TGM#ڕSj 1Sjn8)W|˰VP: tEXtDepth'# IENDB`PK7H! nnFfityk-v1.3.0/_images/math/c363846f9114db62ae575be00f981772ccd53f96.pngPNG  IHDRB #)0PLTE/ݠ{tRNSHt`2z. pHYs+IDATH͖?hQǿw&w5A ]ê`tP"B[TQ!I%CT\P*8f)h){w ^^ &z=df$c4~l'&9D)$3q82$as4~:֫T3Q@SN&߹+X4$^G f6IxZ+)~0IVuUZHސl>:ҕ)>^bJj>}>ڟ0.Σ˞mR}PLk"vL]zi=I_&`'KO[z^ޫ;M X8>$VĄdQKt"c3X81PbL%=66ׂ;,eգd_ފu+_^U`J,T*[|I$lxJ(+%aY &Kkaw{Y9 U"gJbQJfMHJ.Ox3nJ 7qUB1qt٪7+OWfGHJvШ݈ИJNY1P|DJFXMvDbXKDwIr?g<]#B?JJ> K͎X%R9F*H%U*5qlk  { tEXtDepthIENDB`PK8HTCŬII#fityk-v1.3.0/_static/searchtools.js/* * searchtools.js_t * ~~~~~~~~~~~~~~~~ * * Sphinx JavaScript utilties for the full-text search. * * :copyright: Copyright 2007-2016 by the Sphinx team, see AUTHORS. * :license: BSD, see LICENSE for details. * */ /* Non-minified version JS is _stemmer.js if file is provided */ /** * Porter Stemmer */ var Stemmer = function() { var step2list = { ational: 'ate', tional: 'tion', enci: 'ence', anci: 'ance', izer: 'ize', bli: 'ble', alli: 'al', entli: 'ent', eli: 'e', ousli: 'ous', ization: 'ize', ation: 'ate', ator: 'ate', alism: 'al', iveness: 'ive', fulness: 'ful', ousness: 'ous', aliti: 'al', iviti: 'ive', biliti: 'ble', logi: 'log' }; var step3list = { icate: 'ic', ative: '', alize: 'al', iciti: 'ic', ical: 'ic', ful: '', ness: '' }; var c = "[^aeiou]"; // consonant var v = "[aeiouy]"; // vowel var C = c + "[^aeiouy]*"; // consonant sequence var V = v + "[aeiou]*"; // vowel sequence var mgr0 = "^(" + C + ")?" + V + C; // [C]VC... is m>0 var meq1 = "^(" + C + ")?" + V + C + "(" + V + ")?$"; // [C]VC[V] is m=1 var mgr1 = "^(" + C + ")?" + V + C + V + C; // [C]VCVC... is m>1 var s_v = "^(" + C + ")?" + v; // vowel in stem this.stemWord = function (w) { var stem; var suffix; var firstch; var origword = w; if (w.length < 3) return w; var re; var re2; var re3; var re4; firstch = w.substr(0,1); if (firstch == "y") w = firstch.toUpperCase() + w.substr(1); // Step 1a re = /^(.+?)(ss|i)es$/; re2 = /^(.+?)([^s])s$/; if (re.test(w)) w = w.replace(re,"$1$2"); else if (re2.test(w)) w = w.replace(re2,"$1$2"); // Step 1b re = /^(.+?)eed$/; re2 = /^(.+?)(ed|ing)$/; if (re.test(w)) { var fp = re.exec(w); re = new RegExp(mgr0); if (re.test(fp[1])) { re = /.$/; w = w.replace(re,""); } } else if (re2.test(w)) { var fp = re2.exec(w); stem = fp[1]; re2 = new RegExp(s_v); if (re2.test(stem)) { w = stem; re2 = /(at|bl|iz)$/; re3 = new RegExp("([^aeiouylsz])\\1$"); re4 = new RegExp("^" + C + v + "[^aeiouwxy]$"); if (re2.test(w)) w = w + "e"; else if (re3.test(w)) { re = /.$/; w = w.replace(re,""); } else if (re4.test(w)) w = w + "e"; } } // Step 1c re = /^(.+?)y$/; if (re.test(w)) { var fp = re.exec(w); stem = fp[1]; re = new RegExp(s_v); if (re.test(stem)) w = stem + "i"; } // Step 2 re = /^(.+?)(ational|tional|enci|anci|izer|bli|alli|entli|eli|ousli|ization|ation|ator|alism|iveness|fulness|ousness|aliti|iviti|biliti|logi)$/; if (re.test(w)) { var fp = re.exec(w); stem = fp[1]; suffix = fp[2]; re = new RegExp(mgr0); if (re.test(stem)) w = stem + step2list[suffix]; } // Step 3 re = /^(.+?)(icate|ative|alize|iciti|ical|ful|ness)$/; if (re.test(w)) { var fp = re.exec(w); stem = fp[1]; suffix = fp[2]; re = new RegExp(mgr0); if (re.test(stem)) w = stem + step3list[suffix]; } // Step 4 re = /^(.+?)(al|ance|ence|er|ic|able|ible|ant|ement|ment|ent|ou|ism|ate|iti|ous|ive|ize)$/; re2 = /^(.+?)(s|t)(ion)$/; if (re.test(w)) { var fp = re.exec(w); stem = fp[1]; re = new RegExp(mgr1); if (re.test(stem)) w = stem; } else if (re2.test(w)) { var fp = re2.exec(w); stem = fp[1] + fp[2]; re2 = new RegExp(mgr1); if (re2.test(stem)) w = stem; } // Step 5 re = /^(.+?)e$/; if (re.test(w)) { var fp = re.exec(w); stem = fp[1]; re = new RegExp(mgr1); re2 = new RegExp(meq1); re3 = new RegExp("^" + C + v + "[^aeiouwxy]$"); if (re.test(stem) || (re2.test(stem) && !(re3.test(stem)))) w = stem; } re = /ll$/; re2 = new RegExp(mgr1); if (re.test(w) && re2.test(w)) { re = /.$/; w = w.replace(re,""); } // and turn initial Y back to y if (firstch == "y") w = firstch.toLowerCase() + w.substr(1); return w; } } /** * Simple result scoring code. */ var Scorer = { // Implement the following function to further tweak the score for each result // The function takes a result array [filename, title, anchor, descr, score] // and returns the new score. /* score: function(result) { return result[4]; }, */ // query matches the full name of an object objNameMatch: 11, // or matches in the last dotted part of the object name objPartialMatch: 6, // Additive scores depending on the priority of the object objPrio: {0: 15, // used to be importantResults 1: 5, // used to be objectResults 2: -5}, // used to be unimportantResults // Used when the priority is not in the mapping. objPrioDefault: 0, // query found in title title: 15, // query found in terms term: 5 }; /** * Search Module */ var Search = { _index : null, _queued_query : null, _pulse_status : -1, init : function() { var params = $.getQueryParameters(); if (params.q) { var query = params.q[0]; $('input[name="q"]')[0].value = query; this.performSearch(query); } }, loadIndex : function(url) { $.ajax({type: "GET", url: url, data: null, dataType: "script", cache: true, complete: function(jqxhr, textstatus) { if (textstatus != "success") { document.getElementById("searchindexloader").src = url; } }}); }, setIndex : function(index) { var q; this._index = index; if ((q = this._queued_query) !== null) { this._queued_query = null; Search.query(q); } }, hasIndex : function() { return this._index !== null; }, deferQuery : function(query) { this._queued_query = query; }, stopPulse : function() { this._pulse_status = 0; }, startPulse : function() { if (this._pulse_status >= 0) return; function pulse() { var i; Search._pulse_status = (Search._pulse_status + 1) % 4; var dotString = ''; for (i = 0; i < Search._pulse_status; i++) dotString += '.'; Search.dots.text(dotString); if (Search._pulse_status > -1) window.setTimeout(pulse, 500); } pulse(); }, /** * perform a search for something (or wait until index is loaded) */ performSearch : function(query) { // create the required interface elements this.out = $('#search-results'); this.title = $('

' + _('Searching') + '

').appendTo(this.out); this.dots = $('').appendTo(this.title); this.status = $('

').appendTo(this.out); this.output = $('