# Basic Usage¶

In this section, we will illustrate how to use XAD to compute first order derivatives in both forward and adjoint mode.

As an example, we choose a simple function with 4 inputs and 1 output variable, defined as:

``````double f(double x0, double x1, double x2, double x3)
{
double a = sin(x0) * cos(x1);
double b = x2 * x3 - tan(x1 - x2);
double c = a + 2* c;
return c*c;
}
``````

We will compute derivatives of this function at the point:

``````double x0 = 1.0;
double x1 = 1.5;
double x2 = 1.3;
double x3 = 1.2;
``````

## Prerequisite: Replace Active Variables¶

In order to use XAD to differentiate this function, we first must replace all independent data types and all values that depend on them with an active data type provided by XAD. In the above function, all variables depend on the inputs and thus all occurrences of `double` must be replaced.

This can be done in one of two ways:

1. The variables can be replaced directly, given the desired mode of differentiation. For example, for forward mode `double` is replaced by the type `FReal` and for adjoint mode the type `AReal`.

2. The function is made a template, so that it can be called with any data type, including the original `double`.

We choose the second approach for this tutorial, thus the function becomes:

``````template <class T>
T f(T x0, T x1, T x2, T x3)
{
T a = sin(x0) * cos(x1);
T b = x2 * x3 - tan(x1 - x2);
T c = a + 2* b;
return c*c;
}
``````

This means we can use the same definition with both forward and adjoint modes.

## Forward Mode¶

As illustrated in Forward Mode, when applied to a function with a single output, the forward mode of algorithmic differentiation can compute one derivative at a time. For illustration, we choose to derive the function with respect to the input variable `x0`.

To initiate the forward mode, we must first declare active variables with the appropriate type. XAD provides convenience typedefs to select the mode of differentiation, illustrated in detail in AD Mode Interface. For forward mode, we can declare the types needed as:

``````typedef xad::fwd<double> mode;
``````

We can then use the `AD` typedef for our variables.

The next step is to initialize the dependent variables, which is simply done by assigning the input values to new variables of type `AD`:

``````AD x0_ad = x0;
``````

For forward mode, we must now seed the initial derivative for the variable we are interested in with the value 1 (as described in Forward Mode), as:

``````derivative(x0_ad) = 1.0;
``````

The global function `derivative()` is a convenience function that works on any active data type. Alternatively, we could have used the member function `FReal::setDerivative()`.

At this point we are ready to call our function and it will compute the function value as well as the derivative we are interested in:

``````AD y = f(x0_ad, x1_ad, x2_ad, x3_ad);
``````

We can now access the results using the `value()` and `derivative()` functions on the output (or the member functions `FReal::getDerivative()` and `FReal::getValue()`). For example, the following code outputs them to the console:

``````std::cout << "y = " << value(y) << "\n"
<< "dy/dx0 = " << derivative(y) << "\n";
``````

This example is included with XAD (fwd_1st).

The adjoint mode of automatic differentiation is the natural choice for the function at hand, as it has a single output and multiple inputs. We can get all four derivatives in one execution.

Adjoint mode needs a tape to record the operations and their values during the valuation. After setting the adjoints of the outputs, this tape can then be rolled back to compute the adjoints of the inputs.

Both the active data type and the tape type can be obtained from the interface structure `adj` (see AD Mode Interface for details):

``````typedef xad::adj<double> mode;
typedef mode::tape_type tape_type;
``````

The first step for computing adjoints is to initialise the tape:

``````tape_type tape;
``````

This calls the default constructor `Tape::Tape()`, which creates the tape and activates it. Before creating active variables (type `AD`), the current threads needs to have an active tape. To ensure thread-safety, every thread of the application can have its own active tape.

Next, we create the input variables and register them with the tape:

``````AD x0_ad = x0;
tape.registerInput(x0);
tape.registerInput(x1);
tape.registerInput(x2);
tape.registerInput(x3);
``````

Note that only variables registered as inputs with the tape and all variables dependent on them are recorded.

Once the independent variables are set, we can start recording derivatives on tape and run the algorithm:

``````tape.newRecording();

``````

At this stage, we have all operations recorded and have the value computed. We now need to register the outputs with the tape as well, before we can seed the initial adjoint of the output wit 1 as explained in Adjoint Mode:

``````tape.registerOutput(y);
derivative(y) = 1.0;
``````

This uses the free function `derivative()`, which returns a reference to the stored derivative (or adjoint) of the given parameter. Alternatively the member functions `AReal::setAdjoint()` or `AReal::setDerivative()` can be used for the same purpose.

What is left is interpreting the tape to compute the adjoints of the independent variables:

``````tape.computeAdjoints();
``````

We can now access the adjoints of the inputs, which are the derivatives we are interested in, via the global `derivative()` function or the member function `AReal::getDerivative()`:

``````std::cout << "y     = " << value(y) << "\n"
<< "dy/dx0 = " << derivative(x0_ad) << "\n"
<< "dy/dx1 = " << derivative(x1_ad) << "\n"
<< "dy/dx2 = " << derivative(x2_ad) << "\n"
<< "dy/dx3 = " << derivative(x3_ad) << "\n";
``````