My stuff
Imitation-Projected Programmatic Reinforcement Learning

We study the problem of programmatic reinforcement learning, in which policies are represented as short programs in a symbolic language. Programmatic policies can be more interpretable, generalizable, and amenable to formal verification than neural policies; however, designing rigorous learning approaches for such policies remains a challenge. Our approach to this challenge — a meta-algorithm called PROPEL— is based on three insights. First, we view our learning task as optimization in policy space, modulo the constraint that the desired policy has a programmatic representation, and solve this optimization problem using a form of mirror descent that takes a gradient step into the unconstrained policy space and then projects back onto the constrained space. Second, we view the unconstrained policy space as mixing neural and programmatic representations, which enables employing state-of-the-art deep policy gradient approaches. Third, we cast the projection step as program synthesis via imitation learning, and exploit contemporary combinatorial methods for this task. We present theoretical convergence results for PROPEL and empirically evaluate the approach in three continuous control domains. The experiments show that PROPEL can significantly outperform state-of-the-art approaches for learning programmatic policies.

A growing body of work [58, 8, 60] investigates reinforcement learning (RL) approaches that represent policies as programs in a symbolic language, e.g., a domain-specific language for composing control modules such as PID controllers [5]. Short programmatic policies offer many advantages over neural policies discovered through deep RL, including greater interpretability, better generalization to unseen environments, and greater amenability to formal verification. These benefits motivate developing effective approaches for learning such programmatic policies.

However, programmatic reinforcement learning (PRL) remains a challenging problem, owing to the highly structured nature of the policy space. Recent state-of-the-art approaches employ program synthesis methods to imitate or distill a pre-trained neural policy into short programs [58, 8]. However, such a distillation process can yield a highly suboptimal programmatic policy — i.e., a large distillation gap — and the issue of direct policy search for programmatic policies also remains open.

In this paper, we develop PROPEL (Imitation-Projected Programmatic Reinforcement Learning), a new learning meta-algorithm for PRL, as a response to this challenge. The design of PROPEL is based on three insights that enables integrating and building upon state-of-the-art approaches for policy gradients and program synthesis. First, we view programmatic policy learning as a constrained policy optimization problem, in which the desired policies are constrained to be those that have a programmatic representation. This insight motivates utilizing constrained mirror descent approaches, which take a gradient step into the unconstrained policy space and then project back onto the constrained space. Second, by allowing the unconstrained policy space to have a mix of neural

π(s) ::= a | Op(π1(s), . . . , πk(s)) | if b then π1(s) else π2(s) | ⊕θ(π1(s), . . . , πk(s))b ::= φ(s) | BOp(b1, . . . , bk)

Figure 1: A high-level syntax for programmatic policies, inspired by [58]. A policy  π(s)takes a state s as input and produces an action a as output. b represents boolean expressions;  φis a boolean-valued operator on states; Op is an operator that combines multiple policies into one policy; BOp is a standard boolean operator; and  ⊕θis a “library function" parameterized by  θ.


Figure 2: A programmatic policy for acceleration in TORCS [59], automatically discovered by PROPEL. s[TrackPos] represents the most recent reading from sensor TrackPos.

and programmatic representations, we can employ well-developed deep policy gradient approaches [55, 36, 47, 48, 19] to compute the unconstrained gradient step. Third, we define the projection operator using program synthesis via imitation learning [58, 8], in order to recover a programmatic policy from the unconstrained policy space. Our contributions can be summarized as:

We present PROPEL, a novel meta-algorithm that is based on mirror descent, program synthesis, and imitation learning, for PRL.

On the theoretical side, we show how to cast PROPEL as a form of constrained mirror descent. We provide a thorough theoretical analysis characterizing the impact of approximate gradients and projections. Further, we prove results that provide expected regret bounds and finite-sample guarantees under reasonable assumptions.

On the practical side, we provide a concrete instantiation of PROPEL and evaluate it in three continuous control domains, including the challenging car-racing domain TORCS [59]. The experiments show significant improvements over state-of-the-art approaches for learning programmatic policies.

The problem of programmatic reinforcement learning (PRL) consists of a Markov Decision Process (MDP) M and a programmatic policy class  Π. The definition of  M = (S, A, P, c, p0, γ)is standard [54], with S being the state space, A the action space,  P(s′|s, a)the probability density function of transitioning from a state-action pair to a new state, c(s, a) the state-action cost function,  p0(s)a distribution over starting states, and  γ ∈ (0, 1)the discount factor. A policy  π : S → A(stochastically) maps states to actions. We focus on continuous control problems, so S and A are assumed to be continuous spaces. The goal is to find a programmatic policy  π∗ ∈ Πsuch that:


with the expectation taken over the initial state distribution  s0 ∼ p0, the policy decisions, and the transition dynamics P. One can also use rewards, in which case (1) becomes a maximization problem.

Programmatic Policy Class. A programmatic policy class  Πconsists of policies that can be represented parsimoniously by a (domain-specific) programming language. Recent work [58, 8, 60] indicates that such policies can be easier to interpret and formally verify than neural policies, and can also be more robust to changes in the environment.

In this paper, we consider two concrete classes of programmatic policies. The first, a simplification of the class considered in Verma et al. [58], is defined by the modular, high-level language in Figure 1. This language assumes a library of parameterized functions  ⊕θrepresenting standard controllers, for instance Proportional-Integral-Derivative (PID) [6] or bang-bang controllers [11]. Programs in the language take states s as inputs and produce actions a as output, and can invoke fully instantiated library controllers along with predefined arithmetic, boolean and relational operators. The second, “lower-level" class, from Bastani et al. [8], consists of decision trees that map states to actions.

Example. Consider the problem of learning a programmatic policy, in the language of Figure 1, that controls a car’s accelerator in the TORCS car-racing environment [59]. Figure 2 shows a program in our language for this task. The program invokes PID controllers  PID⟨j,θP ,θI,θD⟩, where j identifies


the sensor (out of 29, in our experiments) that provides inputs to the controller, and  θP , θI, and  θDare respectively the real-valued coefficients of the proportional, integral, and derivative terms in the controller. We note that the program only uses the sensors TrackPos and RPM. While TrackPos (for the position of the car relative to the track axis) is used to decide which controller to use, only the RPM sensor is needed to calculate the acceleration.

Learning Challenges. Learning programmatic policies in the continuous RL setting is challenging, as the best performing methods utilize policy gradient approaches [55, 36, 47, 48, 19], but policy gradients are hard to compute in programmatic representations. In many cases,  Πmay not even be differentiable. For our approach, we only assume access to program synthesis methods that can select a programmatic policy  π ∈ Πthat minimizes imitation disagreement with demonstrations provided by a teaching oracle. Because imitation learning tends to be easier than general RL in long-horizon tasks [53], the task of imitating a neural policy with a program is, intuitively, significantly simpler than the full programmatic RL problem. This intuition is corroborated by past work on programmatic RL [58], which shows that direct search over programs often fails to meet basic performance objectives.

To develop our approach, we take the viewpoint of (1) being a constrained optimization problem, where  Π ⊂ Hresides within a larger space of policies H. In particular, we will represent  H ≡ Π⊕Fusing a mixing of programmatic policies  Πand neural polices F. Any mixed policy  h ≡ π + f canbe invoked as  h(s) = π(s) + f(s). In general, we assume that F is a good approximation of  Π (i.e.,for each  π ∈ Πthere is some  f ∈ Fthat approximates it well), which we formalize in Section 4.

We can now frame our constrained learning problem as minimizing (1) over  Π ⊂ H, that alternate between taking a gradient step in the general space H and projecting back down onto  Π. This “lift-and-project” perspective motivates viewing our problem via the lens of mirror descent [40]. In standard mirror descent, the unconstrained gradient step can be written as  h ← hprev −η∇HJ(hprev) for stepsize  η, and the projection can be written as  π ← argminπ′∈Π D(π′, h)for divergence measure D.

Our approach, Imitation-Projected Programmatic Reinforcement Learning (PROPEL), is outlined in Algorithm 1 (also see Figure 3). PROPEL is a meta-algorithm that requires instantiating two subroutines, UPDATE and PROJECT, which correspond to the standard update and projection steps, respectively. PROPEL can be viewed as a form of functional mirror descent with some notable deviations from vanilla mirror descent.


Figure 3: Depicting the PROPEL meta-algorithm.


UPDATEF.Since policy gradient methods are well-developed for neural policy classes F (e.g., [36, 47, 48, 30, 24, 19]) and non-existent for programmatic policy classes  Π,PROPEL is designed to leverage policy gradients in F and avoid policy gradients in  Π. Algorithm 2 shows one instantiation of UPDATEF. Note that standard mirror descent takes unconstrained gradient steps in H rather than F, and we discuss this discrepancy between UPDATEH and UPDATEFin Section 4.

PROJECTΠ.Projecting onto  Πcan be implemented using program synthesis via imitation learning, i.e., by synthesizing a  π ∈ Πto best imitate demonstrations provided by a teaching oracle  h ∈ H. Recent work [58, 8, 60] has given practical heuristics for this task for various


programmatic policy classes. Algorithm 3 shows one instantiation of PROJECTΠ(based on DAgger [46]). One complication that arises is that finite-sample runs of such imitation learning approaches only return approximate solutions and so the projection is not exact. We characterize the impact of approximate projections in Section 4.

Practical Considerations. In practice, we often employ multiple gradient steps before taking a projection step (as also described in Algorithm 2), because the step size of individual (stochastic) gradient updates can be quite small. Another issue that arises in virtually all policy gradient approaches is that the gradient estimates can have very high variance [55, 33, 30]. We utilize low-variance policy gradient updates by using the reference  πas a proximal regularizer in function space [19].

For the projection step (Algorithm 3), in practice we often retain all previous roll-outs  τfrom all previous projection steps. It is straightforward to query the current oracle h to provide demonstrations on the states  s ∈ τfrom previous roll-outs, which can lead to substantial savings in sample complexity with regards to executing roll-outs on the environment, while not harming convergence.

We start by viewing PROPEL through the lens of online learning in function space, independent of the specific parametric representation. This start point yields a convergence analysis of Alg. 1 in Section 4.1 under generic approximation errors. We then analyze the issues of policy class representation in Sections 4.2 and 4.3, and connect Algorithms 2 and 3 with the overall performance, under some simplifying conditions. In particular, Section 4.3 characterizes the update error in a possibly non-differentiable setting; to our knowledge, this is the first such analysis of its kind for reinforcement learning.

Preliminaries. We consider  Πand F to be subspaces of an ambient policy space U, which is a vector space equipped with inner product  ⟨·, ·⟩, induced norm  ∥u∥ =�⟨u, u⟩, dual norm  ∥v∥∗ =sup{⟨v, u⟩| ∥u∥ ≤ 1}, and standard scaling & addition:  (au + bv)(s) = au(s) + bv(s) for a, b ∈ Rand  u, v ∈ U. The cost functional of a policy u is  J(u) =�S c(s, u(s))dµu(s), where  µuis the distribution of states induced by u. The joint policy class is  H = Π ⊕ F, by  H = {π + f|∀π ∈Π, f ∈ F}.2 Note that H is a subspace of U, and inherits its vector space properties. Without affecting the analysis, we simply equate  U ≡ Hfor the remainder of the paper.

We assume that J is convex in H, which implies that subgradient  ∂J(h)exists (with respect to H) [9]. Where J is differentiable, we utilize the notion of a Fréchet gradient. Recall that a bounded linear operator  ∇ : H �→ His called a Fréchet functional gradient of J at  h ∈ Hif


respect to  H, whereas ∇Fdefines the gradient in the restricted subspace F.

4.1 PROPEL as (Approximate) Functional Mirror Descent

For our analysis, PROPEL can be viewed as approximating mirror descent in (infinite-dimensional) function space over a convex set  Π ⊂ H.3 Similar to the finite-dimensional setting [40], we choose a strongly convex and smooth functional regularizer R to be the mirror map. From the approximate mirror descent perspective, for each iteration t:


DR(u, v) = R(u) − R(v) − ⟨∇R(u), u − v⟩is a Bregman divergence. Taking  R(h) = 12 ∥h∥2will recover projected functional gradient descent in L2-space. Here UPDATE becomes  ht =πt−1 − η �∇J(πt−1), and PROJECT solves for  argminπ∈Π ∥π − ht∥2. While we mainly focus on this choice of R in our experiments, note that other selections of R lead to different UPDATE and PROJECT operators (e.g., minimizing KL divergence if R is negative entropy).

The functional mirror descent scheme above may encounter two additional sources of error compared to standard mirror descent [40]. First, in the stochastic setting (also called bandit feedback [28]), the gradient estimate �∇tmay be biased, in addition to having high variance. One potential source of bias is the gap between UPDATEHand UPDATEF. Second, the PROJECT step may be inexact. We start by analyzing the behavior of PROPEL under generic bias, variance, and projection errors, before discussing the implications of approximating UPDATEH and PROJECTΠby Algs. 2 & 3, respectively. Let the bias be bounded by  β, i.e.,���E[�∇t|πt] − ∇J(πt)���∗ ≤ βalmost surely. Similarly let the variance of the gradient estimate be bounded by  σ2, and the projection error norm  ∥πt − π∗t ∥ ≤ ϵ. We state the expected regret bound below; more details and a proof appear in Appendix A.2.

Theorem 4.1 (Expected regret bound under gradient estimation and projection errors). Let  π1, . . . , πTbe a sequence of programmatic policies returned by Algorithm 1, and  π∗ be the optimal programmatic policy. Choosing learning rate  η =�


The result shows that error  ϵ from PROJECTand the bias  βdo not accumulate and simply contribute an additive term on the expected regret.4 The effect of variance of gradient estimate decreases at a �1/Trate. Note that this regret bound is agnostic to the specific UPDATE and PROJECT operations, and can be applied more generically beyond the specific algorithmic choices used in our paper.

4.2 Finite-Sample Analysis under Vanilla Policy Gradient Update and DAgger Projection

Next, we show how certain instantiations of UPDATE and PROJECT affect the magnitude of errors and influence end-to-end learning performance from finite samples, under some simplifying assumptions on the UPDATE step. For this analysis, we simplify Alg. 2 into the case UPDATEF ≡ UPDATEH. In particular, we assume programmatic policies in  Πto be parameterized by a vector  θ ∈ Rk, and πis differentiable in  θ(e.g., we can view  Π ⊂ Fwhere F is parameterized in  Rk). We further assume the trajectory roll-out is performed in an exploratory manner, where action is taken uniformly random over finite set of A actions, thus enabling the bound on the bias of gradient estimates via Bernstein’s inequality. The PROJECT step is consistent with Alg. 3, i.e., using DAgger [45] under convex imitation loss, such as  ℓ2loss. We have the following high-probability guarantee:

Theorem 4.2 (Finite-sample guarantee). At each iteration, we perform vanilla policy gradient estimate of  π (over H) using mtrajectories and, use DAgger algorithm to collect M roll-outs for the


imitation learning projection. Setting the learning rate  η =

rounds of the algorithm, we have that:


holds with probability at least  1 − δ, with Hbeing the task horizon, A the cardinality of action space, σ2 the variance of policy gradient estimates, and k the dimension  Π’s parameterization.

The expanded result and proof are included in Appendix A.3. The proof leverages previous analysis from DAgger [46] and the finite sample analysis of vanilla policy gradient algorithm [32]. The finite-sample regret bound scales linearly with the standard deviation  σof the gradient estimate, while the bias, which is the very last component of the RHS, scales linearly with the task horizon H. Note that the standard deviation  σcan be exponential in task horizon H in the worst case [32], and so it is important to have practical implementation strategies to reduce the variance of the UPDATE operation. While conducted in a stylized setting, this analysis provides insight in the relative trade-offs of spending effort in obtaining more accurate projections versus more reliable gradient estimates.

4.3 Closing the gap between UPDATEH and UPDATEF

Our functional mirror descent analysis rests on taking gradients in H: UPDATEH(π)involves estimating  ∇HJ(π) in the Hspace. On the other hand, Algorithm 2 performs UPDATEF(π) only inthe neural policy space F. In either case, although  J(π)may be differentiable in the non-parametric ambient policy space, it may not be possible to obtain a differentiable parametric programmatic representation in  Π. In this section, we discuss theoretical motivations to addressing a practical issue: How do we define and approximate the gradient  ∇HJ(π)under a parametric representation? To our knowledge, we are the first to consider such a theoretical question for reinforcement learning.

Defining a consistent approximation of  ∇HJ(π).The idea in UPDATEF(π)(Line 8 of Alg. 1) is to approximate  ∇HJ(π) by ∇FJ(f), which has a differentiable representation, at some  f close to π(under the norm). Under appropriate conditions on F, we show that this approximation is valid.

Proposition 4.3. Assume that (i) J is Fréchet differentiable on H, (ii) J is also differentiable on the restricted subspace F, and (iii) F is dense in H (i.e., the closure F = H). Then for any fixed policy  π ∈ Π, define a sequence of policies  fk ∈ F, k = 1, 2, . . .), that converges to  π: limk→∞ ∥fk − π∥ = 0. We then have  limk→∞ ∥∇FJ(fk) − ∇HJ(π)∥∗ = 0.

Since the Fréchet gradient is unique in the ambient space  H, ∀k we have ∇HJ(fk) = ∇FJ(fk) →∇HJ(π) as k → ∞(by Proposition 4.3). We thus have an asymptotically unbiased approximation of ∇HJ(π)via differentiable space  F as: ∇FJ(π) ≜ ∇HJ(π) ≜ limk→∞ ∇FJ(fk).5 Connecting to the result from Theorem 4.1, let  σ2 be an upper bound on the policy gradient estimates in the neural policy class F, under an asymptotically unbiased approximation of  ∇HJ(π), the expected regret bound becomes  E�1T�Tt=1 J(πt)�− J(π∗) = O�σ�

Bias-variance considerations of UPDATEF(π)To further theoretically motivate a practical strategy for UPDATEF(π)in Algorithm 2, we utilize an equivalent proximal perspective of mirror descent [10], where UPDATEH(π)is equivalent to solving for  h′ = argminh∈H η⟨∇HJ(π), h⟩ + DR(h, π).

Proposition 4.4 (Minimizing a relaxed objective). For a fixed programmatic policy  π, with sufficiently small constant  λ ∈ (0, 1), we have that


Thus, a relaxed UPDATEHstep is obtained by minimizing the RHS of (3), i.e., minimizing  J(π + λf)over  f ∈ F. Each gradient descent update step is now  f ′ = f − ηλ∇FJ(πt + λf), corresponding to Line 5 of Algorithm 2. For fixed  πand small  λ, this relaxed optimization problem becomes regularized policy optimization over F, which is significantly easier. Functional regularization in policy space around a fixed prior controller  πhas demonstrated significant reduction in the variance of gradient estimate [19], at the expense of some bias. The below expected regret bound summarizes the impact of this increased bias and reduced variance, with details included in Appendix A.5.

Proposition 4.5 (Bias-variance characterization of UPDATEF). Assuming J(h) is L-strongly smooth over  H, i.e., ∇HJ(h) is L-Lipschitz continuous, approximating UPDATEH by UPDATEFper Alg. 2 leads to the expected regret bound:  E�1T�Tt=1 J(πt)�− J(π∗) = O�λσ�

Compared to the idealized unbiased approximation in Proposition 4.3, the introduced bias here is related to the inherent smoothness property of cost functional J(h) over the joint policy class H, i.e., how close  J(π + λf)is to its linear under-approximation  J(π) + ⟨∇HJ(π), λf⟩ around π.

We demonstrate the effectiveness of PROPEL in synthesizing programmatic controllers in three continuous control environments. For brevity and focus, this section primarily focuses on TORCS6, a challenging race car simulator environment [59]. Empirical results on two additional classic control tasks, Mountain-Car and Pendulum, are provided in Appendix B; those results follow similar trends as the ones described for TORCS below, and further validate the convergence analysis of PROPEL.


Experimental Setup. We evaluate over five distinct tracks in the TORCS simulator. The difficulty of a track can be characterized by three properties; track length, track width, and number of turns. Our suite of tracks provides environments with varying levels of difficulty for the learning algorithm. The performance of a policy in the TORCS simulator is measured by the lap time achieved on the track. To calculate the lap time, the policies are allowed to complete a three-lap race, and we record the best lap time during this race. We perform the experiments with twenty-five random seeds and report the median lap time over these twenty-five trials. Some of the policies crash the car before completing a lap on certain tracks, even after training for 600 episodes. Such crashes are recorded as a lap

time of infinity while calculating the median. If the policy crashes for more than half the seeds, this is reported as CR in Tables 1 & 2. We choose to report the median because taking the crash timing as infinity, or an arbitrarily large constant, heavily skews other common measures such as the mean.


Baselines. Among recent state-of-the-art approaches to learning programmatic policies are NDPS [58] for high-level language policies, and VIPER [8] for learning tree-based policies. Both NDPS and VIPER rely on imitating a fixed (pre-trained) neural policy oracle, and can be viewed as degenerate versions of PROPEL that only run Lines 4-6 in Algorithm 1. We present two PROPEL analogues to NDPS and VIPER: (i) PRO- PELPROG: PROPEL using the high-level language of Figure 1 as the class of programmatic policies, similar to NDPS. (ii) PROPELTREE: PROPEL using regression trees, similar to VIPER. We also report results for PRIOR, which is a (sub-optimal) PID controller that is also used as the initial policy in PROPEL. In addition, to study generalization ability as well as safety behavior

during training, we also include DDPG, a neural policy learned using the Deep Deterministic Policy Gradients [36] algorithm, with 600 episodes of training for each track. In principle, PROPEL and its analysis can accommodate different policy gradient subroutines. However, in the TORCS domain, other policy gradient algorithms such as PPO and TRPO failed to learn policies that are able to complete the considered tracks. We thus focus on DDPG as our main policy gradient component.

Table 1: Performance results in TORCS over 25 random seeds. Each entry is formatted as Lap-time / Crash-ratio, reporting median lap time in seconds over all the seeds (lower is better) and ratio of seeds that result in crashes (lower is better). A lap time of CR indicates the agent crashed and could not complete a lap for more than half the seeds.


Table 2: Generalization results in TORCS, where rows are training and columns are testing tracks. Each entry is formatted as PROPELPROG / DDPG, and the number reported is the median lap time in seconds over all the seeds (lower is better). CR indicates the agent crashed and could not complete a lap for more than half the seeds.


Evaluating Performance. Table 1 shows the performance on the considered TORCS tracks. We see that PROPELPROG and PROPELTREE consistently outperform the NDPS [58] and VIPER [8] baselines, respectively. While DDPG outperforms PROPEL on some tracks, its volatility causes it to be unable to learn in some environments, and hence to crash the majority of the time. Figure 4 shows the consistent improvements made over the prior by PROPELPROG, over the iterations of the PROPEL algorithm. Appendix B contains similar results achieved on the two classic control tasks, MountainCar and Pendulum. Figure 5 shows that, compared to DDPG, our approach suffers far fewer crashes while training in TORCS.

Evaluating Generalization. To compare the ability of the controllers to perform on tracks not seen during training, we executed the learned policies on all the other tracks (Table 2). We observe that DDPG crashes significantly more often than PROPELPROG. This demonstrates the generalizability of the policies returned by PROPEL. Generalization results for the PROPELTREE policy are given in the appendix. In general, PROPELTREE policies are more generalizable than DDPG but less than PROPELPROG. On an absolute level, the generalization ability of PROPEL still leaves much room for improvement, which is an interesting direction for future work.

Verifiability of Policies. As shown in prior work [8, 58], parsimonious programmatic policies are more amenable to formal verification than neural policies. Unsurprisingly, the policies generated by PROPELTREE and PROPELPROG are easier to verify than DDPG policies. As a concrete example, we verified a smoothness property of the PROPELPROG policy using the Z3 SMT-solver [21] (more details in Appendix B). The verification terminated in 0.49 seconds.

Initialization. In principle, PROPEL can be initialized with a random program, or a random policy trained using DDPG. In practice, the performance of PROPEL depends to a certain degree on the stability of the policy gradient procedure, which is DDPG in our experiments. Unfortunately, DDPG often exhibits high variance across trials and fares poorly in challenging RL domains. Specifically, in our TORCS experiments, DDPG fails on a number of tracks (similar phenomena have been reported in previous work that experiments on similar continuous control domains [30, 19, 58]). Agents obtained by initializing PROPEL with neural policies obtained via DDPG also fail on multiple tracks. Their performance over the five tracks is reported in Appendix B. In contrast, PROPEL can often finish the challenging tracks when initialized with a very simple hand-crafted programmatic prior.

Program Synthesis. Program synthesis is the problem of automatically searching for a program within a language that fits a given specification [29]. Recent approaches to the problem have leveraged symbolic knowledge about program structure [27], satisfiability solvers [50, 31], and meta-learning techniques [39, 41, 22, 7] to generate interesting programs in many domains [3, 42, 4]. In most prior work, the specification is a logical constraint on the input/output behavior of the target program. However, there is also a growing body of work that considers program synthesis modulo optimality objectives [13, 15, 43], often motivated by machine learning tasks [39, 57, 26, 23, 58, 8, 60]. Synthesis of programs that imitates an oracle has been considered in both the logical [31] and the optimization [58, 8, 60] settings. The projection step in PROPEL builds on this prior work. While our current implementation of this step is entirely symbolic, in principle, the operation can also utilize contemporary techniques for learning policies that guide the synthesis process [39, 7, 49].

Constrained Policy Learning. Constrained policy learning has seen increased interest in recent years, largely due to the desire to impose side guarantees such as stability and safety on the policy’s behavior. Broadly, there are two approaches to imposing constraints: specifying constraints as an additional cost function [1, 35], and explicitly encoding constraints into the policy class [2, 34, 19, 20, 12]. In some cases, these two approaches can be viewed as duals of each other. For instance, recent work that uses control-theoretic policies as a functional regularizer [34, 19] can be viewed from the perspective of both regularization (additional cost) and an explicitly constrained policy class (a specific mix of neural and control-theoretic policies). We build upon this perspective to develop the gradient update step in our approach.

RL using Imitation Learning. There are two ways to utilize imitation learning subroutines within RL. First, one can leverage limited-access or sub-optimal experts to speed up learning [44, 18, 14, 51]. Second, one can learn over two policy classes (or one policy and one model class) to achieve accelerated learning compared to using only one policy class [38, 17, 52, 16]. Our approach has some stylistic similarities to previous efforts [38, 52] that use a richer policy space to search for improvements before re-training the primary policy to imitate the richer policy. One key difference is that our primary policy is programmatic and potentially non-differentiable. A second key difference is that our theoretical framework takes a functional gradient descent perspective — it would be interesting to carefully compare with previous analysis techniques to find a unifying framework.

RL with Mirror Descent. The mirror descent framework has previously used to analyze and design RL algorithms. For example, Thomas et al. [56] and Mahadevan and Liu [37] use composite objective mirror descent, or COMID [25], which allows incorporating adaptive regularizers into gradient updates, thus offering connections to either natural gradient RL [56] or sparsity inducing RL algorithms [37]. Unlike in our work, these prior approaches perform projection into the same native, differentiable representation. Also, the analyses in these papers do not consider errors introduced by hybrid representations and approximate projection operators. However, one can potentially extend our approach with versions of mirror descent, e.g., COMID, that were considered in these efforts.

We have presented PROPEL, a meta-algorithm based on mirror descent, program synthesis, and imitation learning, for programmatic reinforcement learning (PRL). We have presented theoretical convergence results for PROPEL, developing novel analyses to characterize approximate projections and biased gradients within the mirror descent framework. We also validated PROPEL empirically, and show that it can discover interpretable, verifiable, generalizable, performant policies and significantly outperform the state of the art in PRL.

The central idea of PROPEL is the use of imitation learning and combinatorial methods in implementing a projection operation for mirror descent, with the goal of optimization in a functional space that lacks gradients. While we have developed PROPEL in an RL setting, this idea is not restricted to RL or even sequential decision making. Future work will seek to exploit this insight in other machine learning and program synthesis settings.

Acknowledgements. This work was supported in part by United States Air Force Contract # FA8750-19-C-0092, NSF Award # 1645832, NSF Award # CCF-1704883, the Okawa Foundation, Raytheon, PIMCO, and Intel.

[1] Joshua Achiam, David Held, Aviv Tamar, and Pieter Abbeel. Constrained policy optimization. In Proceedings of the 34th International Conference on Machine Learning-Volume 70, pages 22–31. JMLR. org, 2017.

[2] Mohammed Alshiekh, Roderick Bloem, Rüdiger Ehlers, Bettina Könighofer, Scott Niekum, and Ufuk Topcu. Safe reinforcement learning via shielding. In Thirty-Second AAAI Conference on Artificial Intelligence, 2018.

[3] Rajeev Alur, Rastislav Bodík, Eric Dallal, Dana Fisman, Pranav Garg, Garvit Juniwal, Hadas Kress-Gazit, P. Madhusudan, Milo M. K. Martin, Mukund Raghothaman, Shambwaditya Saha, Sanjit A. Seshia, Rishabh Singh, Armando Solar-Lezama, Emina Torlak, and Abhishek Udupa. Syntax-guided synthesis. In Dependable Software Systems Engineering, pages 1–25. 2015.

[4] Rajeev Alur, Arjun Radhakrishna, and Abhishek Udupa. Scaling enumerative program synthesis via divide and conquer. In Tools and Algorithms for the Construction and Analysis of Systems -23rd International Conference, TACAS 2017, Held as Part of the European Joint Conferences on Theory and Practice of Software, ETAPS 2017, Uppsala, Sweden, April 22-29, 2017, Proceedings, Part I, pages 319–336, 2017.

[5] Kiam Heong Ang, Gregory Chong, and Yun Li. Pid control system analysis, design, and technology. IEEE transactions on control systems technology, 13(4):559–576, 2005.

[6] Karl Johan Åström and Tore Hägglund. Automatic tuning of simple regulators with specifica-tions on phase and amplitude margins. Automatica, 20(5):645–651, 1984.

[7] Matej Balog, Alexander L. Gaunt, Marc Brockschmidt, Sebastian Nowozin, and Daniel Tarlow. Deepcoder: Learning to write programs. In 5th International Conference on Learning Representations, ICLR 2017, Toulon, France, April 24-26, 2017, Conference Track Proceedings, 2017.

[8] Osbert Bastani, Yewen Pu, and Armando Solar-Lezama. Verifiable reinforcement learning via policy extraction. In Advances in Neural Information Processing Systems, pages 2494–2504, 2018.

[9] Heinz H Bauschke, Patrick L Combettes, et al. Convex analysis and monotone operator theory in Hilbert spaces, volume 408. Springer, 2011.

[10] Amir Beck and Marc Teboulle. Mirror descent and nonlinear projected subgradient methods for convex optimization. Operations Research Letters, 31(3):167–175, 2003.

[11] Richard Bellman, Irving Glicksberg, and Oliver Gross. On the “bang-bang” control problem. Quarterly of Applied Mathematics, 14(1):11–18, 1956.

[12] Felix Berkenkamp, Matteo Turchetta, Angela Schoellig, and Andreas Krause. Safe model-based reinforcement learning with stability guarantees. In Advances in neural information processing systems, pages 908–918, 2017.

[13] Roderick Bloem, Krishnendu Chatterjee, Thomas A. Henzinger, and Barbara Jobstmann. Better quality in synthesis through quantitative objectives. In Computer Aided Verification, 21st International Conference, CAV 2009, Grenoble, France, June 26 - July 2, 2009. Proceedings, pages 140–156, 2009.

[14] Kai-Wei Chang, Akshay Krishnamurthy, Alekh Agarwal, Hal Daumé III, and John Langford. Learning to search better than your teacher. In International Conference on Machine Learning (ICML), 2015.

[15] Swarat Chaudhuri, Martin Clochard, and Armando Solar-Lezama. Bridging boolean and quantitative synthesis using smoothed proof search. In POPL, pages 207–220, 2014.

[16] Ching-An Cheng, Xinyan Yan, Nathan Ratliff, and Byron Boots. Predictor-corrector policy optimization. In International Conference on Machine Learning (ICML), 2019.

[17] Ching-An Cheng, Xinyan Yan, Evangelos Theodorou, and Byron Boots. Accelerating imitation learning with predictive models. In International Conference on Artificial Intelligence and Statistics (AISTATS), 2019.

[18] Ching-An Cheng, Xinyan Yan, Nolan Wagener, and Byron Boots. Fast policy learning through imitation and reinforcement. In Uncertainty in artificial intelligence, 2019.

[19] Richard Cheng, Abhinav Verma, Gabor Orosz, Swarat Chaudhuri, Yisong Yue, and Joel Burdick. Control regularization for reduced variance reinforcement learning. In International Conference on Machine Learning (ICML), 2019.

[20] Gal Dalal, Krishnamurthy Dvijotham, Matej Vecerik, Todd Hester, Cosmin Paduraru, and Yuval Tassa. Safe exploration in continuous action spaces. arXiv preprint arXiv:1801.08757, 2018.

[21] Leonardo Mendonça de Moura and Nikolaj Bjørner. Z3: An Efficient SMT Solver. In TACAS, pages 337–340, 2008.

[22] Jacob Devlin, Jonathan Uesato, Surya Bhupatiraju, Rishabh Singh, Abdel-rahman Mohamed, and Pushmeet Kohli. Robustfill: Neural program learning under noisy i/o. In Proceedings of the 34th International Conference on Machine Learning-Volume 70, pages 990–998. JMLR. org, 2017.

[23] Tao Du, Jeevana Priya Inala, Yewen Pu, Andrew Spielberg, Adriana Schulz, Daniela Rus, Armando Solar-Lezama, and Wojciech Matusik. Inversecsg: automatic conversion of 3d models to CSG trees. ACM Trans. Graph., 37(6):213:1–213:16, 2018.

[24] Yan Duan, Xi Chen, Rein Houthooft, John Schulman, and Pieter Abbeel. Benchmarking deep reinforcement learning for continuous control. In International Conference on Machine Learning, pages 1329–1338, 2016.

[25] John C Duchi, Shai Shalev-Shwartz, Yoram Singer, and Ambuj Tewari. Composite objective mirror descent. In COLT, pages 14–26, 2010.

[26] Kevin Ellis, Daniel Ritchie, Armando Solar-Lezama, and Josh Tenenbaum. Learning to infer graphics programs from hand-drawn images. In Advances in Neural Information Processing Systems, pages 6059–6068, 2018.

[27] John K. Feser, Swarat Chaudhuri, and Isil Dillig. Synthesizing data structure transformations from input-output examples. In Proceedings of the 36th ACM SIGPLAN Conference on Programming Language Design and Implementation, Portland, OR, USA, June 15-17, 2015, pages 229–239, 2015.

[28] Abraham D Flaxman, Adam Tauman Kalai, and H Brendan McMahan. Online convex optimization in the bandit setting: gradient descent without a gradient. In Proceedings of the sixteenth annual ACM-SIAM symposium on Discrete algorithms, pages 385–394. Society for Industrial and Applied Mathematics, 2005.

[29] Sumit Gulwani, Oleksandr Polozov, and Rishabh Singh. Program synthesis. Foundations and Trends in Programming Languages, 4(1-2):1–119, 2017.

[30] Peter Henderson, Riashat Islam, Philip Bachman, Joelle Pineau, Doina Precup, and David Meger. Deep reinforcement learning that matters. In Thirty-Second AAAI Conference on Artificial Intelligence, 2018.

[31] Susmit Jha, Sumit Gulwani, Sanjit A Seshia, and Ashish Tiwari. Oracle-guided componentbased program synthesis. In Proceedings of the 32nd ACM/IEEE International Conference on Software Engineering-Volume 1, pages 215–224. ACM, 2010.

[32] Sham Machandranath Kakade et al. On the sample complexity of reinforcement learning. PhD thesis, University of London London, England, 2003.

[33] Vijay R Konda and John N Tsitsiklis. Actor-critic algorithms. In Advances in neural information processing systems, pages 1008–1014, 2000.

[34] Hoang M. Le, Andrew Kang, Yisong Yue, and Peter Carr. Smooth imitation learning for online sequence prediction. In International Conference on Machine Learning (ICML), 2016.

[35] Hoang M Le, Cameron Voloshin, and Yisong Yue. Batch policy learning under constraints. In International Conference on Machine Learning (ICML), 2019.

[36] Timothy P Lillicrap, Jonathan J Hunt, Alexander Pritzel, Nicolas Heess, Tom Erez, Yuval Tassa, David Silver, and Daan Wierstra. Continuous control with deep reinforcement learning. arXiv preprint arXiv:1509.02971, 2015.

[37] Sridhar Mahadevan and Bo Liu. Sparse q-learning with mirror descent. In Proceedings of the Twenty-Eighth Conference on Uncertainty in Artificial Intelligence, pages 564–573. AUAI Press, 2012.

[38] William H Montgomery and Sergey Levine. Guided policy search via approximate mirror descent. In Advances in Neural Information Processing Systems, pages 4008–4016, 2016.

[39] Vijayaraghavan Murali, Swarat Chaudhuri, and Chris Jermaine. Neural sketch learning for conditional program generation. In ICLR, 2018.

[40] Arkadii Semenovich Nemirovsky and David Borisovich Yudin. Problem complexity and method efficiency in optimization. 1983.

[41] Emilio Parisotto, Abdel-rahman Mohamed, Rishabh Singh, Lihong Li, Dengyong Zhou, and Pushmeet Kohli. Neuro-symbolic program synthesis. arXiv preprint arXiv:1611.01855, 2016.

[42] Oleksandr Polozov and Sumit Gulwani. Flashmeta: a framework for inductive program synthesis. In Proceedings of the 2015 ACM SIGPLAN International Conference on ObjectOriented Programming, Systems, Languages, and Applications, OOPSLA 2015, part of SPLASH 2015, Pittsburgh, PA, USA, October 25-30, 2015, pages 107–126, 2015.

[43] Veselin Raychev, Pavol Bielik, Martin T. Vechev, and Andreas Krause. Learning programs from noisy data. In Proceedings of the 43rd Annual ACM SIGPLAN-SIGACT Symposium on Principles of Programming Languages, POPL 2016, St. Petersburg, FL, USA, January 20 - 22, 2016, pages 761–774, 2016.

[44] Stephane Ross and J Andrew Bagnell. Reinforcement and imitation learning via interactive no-regret learning. arXiv preprint arXiv:1406.5979, 2014.

[45] Stéphane Ross, Geoffrey Gordon, and Drew Bagnell. A reduction of imitation learning and structured prediction to no-regret online learning. In Proceedings of the fourteenth international conference on artificial intelligence and statistics, pages 627–635, 2011.

[46] Stéphane Ross, Geoffrey J. Gordon, and Drew Bagnell. A reduction of imitation learning and structured prediction to no-regret online learning. In Proceedings of the Fourteenth International Conference on Artificial Intelligence and Statistics, AISTATS 2011, Fort Lauderdale, USA, April 11-13, 2011, pages 627–635, 2011.

[47] John Schulman, Sergey Levine, Pieter Abbeel, Michael Jordan, and Philipp Moritz. Trust region policy optimization. In International Conference on Machine Learning, pages 1889–1897, 2015.

[48] John Schulman, Filip Wolski, Prafulla Dhariwal, Alec Radford, and Oleg Klimov. Proximal policy optimization algorithms. arXiv preprint arXiv:1707.06347, 2017.

[49] Xujie Si, Yuan Yang, Hanjun Dai, Mayur Naik, and Le Song. Learning a meta-solver for syntax-guided program synthesis. In 7th International Conference on Learning Representations, ICLR 2019, New Orleans, LA, USA, May 6-9, 2019, 2019.

[50] Armando Solar-Lezama, Liviu Tancau, Rastislav Bodík, Sanjit A. Seshia, and Vijay A. Saraswat. Combinatorial sketching for finite programs. In ASPLOS, pages 404–415, 2006.

[51] Wen Sun, J Andrew Bagnell, and Byron Boots. Truncated horizon policy search: Combining reinforcement learning & imitation learning. In International Conference on Learning Representations (ICLR), 2018.

[52] Wen Sun, Geoffrey J Gordon, Byron Boots, and J Bagnell. Dual policy iteration. In Advances in Neural Information Processing Systems, pages 7059–7069, 2018.

[53] Wen Sun, Arun Venkatraman, Geoffrey J Gordon, Byron Boots, and J Andrew Bagnell. Deeply aggrevated: Differentiable imitation learning for sequential prediction. In International Conference on Machine Learning (ICML), 2017.

[54] Richard S Sutton and Andrew G Barto. Reinforcement learning: An introduction. MIT press, 2018.

[55] Richard S Sutton, David A McAllester, Satinder P Singh, and Yishay Mansour. Policy gradient methods for reinforcement learning with function approximation. In Advances in neural information processing systems, pages 1057–1063, 2000.

[56] Philip S Thomas, William C Dabney, Stephen Giguere, and Sridhar Mahadevan. Projected natural actor-critic. In Advances in neural information processing systems, pages 2337–2345, 2013.

[57] Lazar Valkov, Dipak Chaudhari, Akash Srivastava, Charles Sutton, and Swarat Chaudhuri. Houdini: Lifelong learning as program synthesis. In Advances in Neural Information Processing Systems, pages 8687–8698, 2018.

[58] Abhinav Verma, Vijayaraghavan Murali, Rishabh Singh, Pushmeet Kohli, and Swarat Chaudhuri. Programmatically interpretable reinforcement learning. In International Conference on Machine Learning, pages 5052–5061, 2018.

[59] Bernhard Wymann, Eric Espié, Christophe Guionneau, Christos Dimitrakakis, Rémi Coulom, and Andrew Sumner. TORCS, The Open Racing Car Simulator. http://www.torcs.org, 2014.

[60] He Zhu, Zikang Xiong, Stephen Magill, and Suresh Jagannathan. An inductive synthesis framework for verifiable reinforcement learning. In ACM Conference on Programming Language Design and Implementation (SIGPLAN), 2019.

A.1 Preliminaries and Notations

We formally define an ambient control policy space U to be a vector space equipped with inner product  ⟨·, ·⟩ : U × U �→ R, which induces a norm  ∥u∥ =�⟨u, u⟩, and its dual norm defined as ∥v∥∗ = sup{⟨v, u⟩| ∥u∥ ≤ 1}. While multiple ways to define the inner product exist, for concreteness we can think of the example of square-integrable stationary policies with  ⟨u, v⟩ =�S u(s)v(s)ds. The addition operator + between two policies  u, v ∈ Uis defined as (u + v)(s) = u(s) + v(s) for all state  s ∈ S. Scaling λu + κvis defined similarly for scalar  λ, κ.

The cost functional of a control policy u is defined as  J(u) =� ∞0 c(s(τ), u(τ))dτ, or J(u) = �S c(s, u(s))dµu(s), where  µuis the distribution of states induced by policy u. This latter example is equivalent to the standard notion of value function in reinforcement learning.

Separate from the parametric representation issues, both programmatic policy class  Πand neural policy class F, and by extension - the joint policy class H, are considered to live in the ambient vector space U. We thus have a common and well-defined notion of distance between policies from different classes.

We make an important distinction between differentiability of J(h) in the ambient policy space (non-parametric), versus differentiability in parameterization (parametric). For example, if  Πis a class of decision-tree based policy, policies in  Πmay not be differentiable under representation. However, policies  π ∈ Πmight still be differentiable when considered as points in the ambient vector space U.

We will use the following standard notion of gradient and differentiability from functional analysis: Definition A.1 (Subgradients). The subgradient of J at h, denoted  ∂J(h), is the non-empty set


Definition A.2 (Fréchet gradient). A bounded linear operator  ∇ : H �→ His called Fréchet functional gradient of  J at h ∈ H if lim∥g∥→0J(h+g)−J(h)−⟨∇J(h),g⟩∥g∥ = 0

The notions of convexity, smoothness and Bregman divergence are analogous to finite-dimensional setting:

Definition A.3 (Strong convexity). A differentiable function  R is α−strongly convex w.r.t norm  ∥·∥if  R(y) ≥ R(x) + ⟨∇R(x), y − x⟩ + α2 ∥y − x∥2

Definition A.4 (Lipschitz continuous gradient smoothness). A differentiable function R is LR−strongly smooth w.r.t norm  ∥·∥ if ∥∇R(x) − ∇R(y)∥∗ ≤ LR ∥x − y∥

Definition A.5 (Bregman Divergence). For a strongly convex regularizer  R, DR(x, y) = R(x) −R(y) − ⟨∇R(y), x − y⟩is the Bregman divergence between x and y (not necessarily symmetric)

The following standard result for Bregman divergence will be useful:

Lemma A.1. [10] For all x, y, z we have the identity  ⟨∇R(x) − ∇R(y), x − z⟩ = DR(x, y) +DR(z, x) − DR(z, y). Since Bregman divergence is non-negative, a consequence of this identity is that  DR(z, x) − DR(z, y) ≤ ⟨∇R(x) − ∇R(y), z − x⟩

A.2 Expected Regret Bound under Noisy Policy Gradient Estimates and Projection Errors

In this section, we show regret bound for the performance of the sequence of returned programs π1, . . . , πTof the algorithm. The analysis here is agnostic to the particular implementation of algorithm 2 and algorithm 3.

Let  R be a α−strongly convex and  LR−smooth functional with respect to norm  ∥·∥ on H. The steps from algorithm 1 can be described as follows.


This procedure is an approximate functional mirror descent scheme under bandit feedback. We will develop the following result, which is a more detailed version of 4.1 in the main paper.

In the statement below, D is the diameter on  Πwith respect to defined norm  ∥·∥(i.e., D = sup ∥π − π′∥). LJis the Lipschitz constant of the functional J on  H. β, σ2are the bound on the bias and variance of the gradient estimate at each iteration, respectively.  α and Rare the strongly convex and smooth coefficients of the functional regularizer R. Finally,  ϵis the bound on the projection error with respect to the same norm  ∥·∥.

Theorem A.2 (Regret bound of returned policies). Let  π1, . . . , πTbe a sequence of programmatic policies returned by algorithm 1 and  π∗be the optimal programmatic policy. We have the expected regret bound:


Proof. At each round t, let  ∇t = E[�∇t|πt]be the conditional expectation of the gradient estimate. We will use the shorthand notation  ∇t = ∇J(πt). Denote the upper-bound on the bias of the estimate by  βt, i.e.,��∇t − ∇t��∗ ≤ βtalmost surely. Denote the noise of the gradient estimate by


The projection operator is  ϵ−approximate in the sense that ��πt − PROJECTRΠ(ft)�� =���� �PROJECTRΠ(ht) − PROJECTRΠ(ht)���� ≤ ϵwith some constant  ϵ, which reflects the statistical error of the imitation learning procedure. This projection error in general is independent of the choice of function classes  Πand F.We will use the shorthand notation  π∗t = PROJECTRΠ(ft)for the true Bregman projection of  ht onto Π.

Due to convexity of J over the space H (which includes  Π), we have for all  π ∈ Π:


We proceed to bound the RHS, starting with bounding the inner product where the actual gradient is replaced by the estimated gradient.


Equation (5) is due to the gradient update rule in F space. Equation (6) is derived from definition of Bregman divergence. Equation (7) is due to the generalized Pythagorean theorem of Bregman projection  DR(x, y) ≥ DR(x, PROJECTRΠ(x)) + DR(PROJECTRΠ(x), y). The RHS of equation (7) are decomposed into three components that will be bounded separately.

Bounding projection error. By lemma (A.1) we have



Equation (10) is due to Cauchy–Schwarz. Equation (11) is due to Lipschitz smoothness of  ∇Rand definition of  ϵ−approximate projection.

Bounding relative improvement. This follows standard argument from analysis of mirror descent



Equation (13) is from the  α−strong convexity property of regularizer R. Equation (14) is by definition of the gradient update. Combining the bounds on the three components and taking expectation, we thus have


Next, the difference between estimated gradient �∇tand actual gradient  ∇tfactors into the bound via Cauchy-Schwarz:


The results can be deduced from equations (16) and (17).

Unbiased gradient estimates. For the case when the gradient estimate is unbiased, assume the variance of the noise of gradient estimates is bounded by  σ2, we have the expected regret bound for all  pi ∈ Π


here to clarify,  LRis the smoothness coefficient of regularizer R (i.e., the gradient of R is  LR-Lipschitz,  LJis Lipschitz constant of J, D is the diameter of  Πunder norm  ∥·∥, σ2is the upper-bound on the variance of gradient estimates, and  ϵis the error from the projection procedure (i.e., imitation learning loss).


Biased gradient estimates. Assume that the bias of gradient estimate at each round is upper-bounded by  βt ≤ β. Similar to before, combining inequalities from (16) and (17), we have


bounded by  O(σ�1T + ϵ + β). Compared to the bound on (18), in the biased case, the extra regret incurred per bound is simply a constant, and does not depend on T.

A.3 Finite-Sample Analysis

In this section, we provide overall finite-sample analysis for PROPEL under some simplifying assumptions. We first consider the case where exact gradient estimate is available, before extending the result to the general case of noisy policy gradient update. Combining the two steps will give us the proof for the following statement (theorem 4.2 in the main paper)

Theorem A.3 (Finite-sample guarantee). At each iteration, we perform vanilla policy gradient estimate of  π (over H) using mtrajectories and use DAgger algorithm to collect M roll-outs. Setting


the learning rate  η =


holds with probability at least  1 − δ, with H the task horizon, A the cardinality of action space,  σ2the variance of policy gradient estimates, and k the dimension  Π’s parameterization.

Exact gradient estimate case. Assuming that the policy gradients can be calculated exactly, it is straight-forward to provide high-probability guarantee for the effect of the projection error. We start with the following result, adapted from [45] for the case of projection error bound. In this version of DAgger, we assume that we only collect a single (state, expert action) pair from each trajectory roll-out. Result is similar, with tighter bound, when multiple data points are collected along the trajectory.

Lemma A.4 (Projection error bound from imitation learning procedure). Using DAgger as the imitation learning sub-routine for our PROJECT operator in algorithm 3, let M be the number of trajectories rolled-out for learning, and H be the horizon of the task. With probability at least  1 − δ,we have


where  πis the result of PROJECT,  π∗is the true Bregman projection of h onto  Π, and  ℓmaxis the maximum value of the imitation learning loss function  DR(·, ·)

The bound in lemma A.4 is simpler than previous imitation learning results with cost information ([44, 45]. The reason is that the goal of the PROJECT operator is more modest. Since we only care about the distance between the empirical projection  πand the true projection  π∗, the loss objective in imitation learning is simplified (i.e., this is only a regret bound), and we can disregard how well policies in  Πcan imitate the expert h, as well as the performance of  J(π)relative to the true cost from the environment J(h).

A consequence of this lemma is that for the number of trajectories at each round of imitation learning M = O( log 1/δϵ2 ) + O( Hϵ ), we have DR(πt, π∗t ) ≤ ϵwith probability at least  1 − δ. Applying union bound across T rounds of learning, we obtain the following guarantee (under no gradient estimation error)

Proposition A.5 (Finite-sample Projection Error Bound). To simplify the presentation of the result, we consider  LR, D, L, αto be known constants. Using DAgger algorithm to collect  M = O( log T/δϵ2 ) +O( Hϵ )roll-outs at each iteration, we have the following regret guarantee after T rounds of our main algorithm:


with probability at least  1 − δ. Consequently, setting  η =


with probability at least  1 − δ

Note that the dependence on the time horizon of the task is sub-linear. This is different from standard imitation learning regret bounds, which are often at least linear in the task horizon. The main reason is that our comparison benchmark  π∗ does live in the space  Π, whereas for DAgger, the expert policy may not reside in the same space.

Noisy gradient estimate case. We now turn to the issue of estimating the gradient of  ∇J(π). We make the following simplifying assumption about the gradient estimation:

The  πis parameterized by vector  θ ∈ Rk (such as a neural network). The parameterization is differentiable with respect to  θ(Alternatively, we can view  Πas a differentiable subspace of F, in which case we have H = F)

At each UPDATE loop, the policy is rolled out m times to collect the data, each trajectory has horizon length H

For each visited state  s ∼ dh, the policy takes a uniformly random action a. The action space is finite with cardinality A.

The gradient  ∇hθis bounded by B

The gradient estimate is performed consistent with a generic policy gradient scheme, i.e.,


Taking uniform random exploratory actions ensures that the samples are i.i.d. We can thus apply Bernstein’s inequality to obtain the bound between estimated gradient and the true gradient. Indeed, with probability at least  1 − δ, we have that the following bound on the bias component-wise:


which leads to similar bound with respect to  ∥·∥∗(here we leverage the equivalence of norms in finite dimensional setting):


Applying union bound of this result over T rounds of learning, and combining with the result from proposition (A.5), we have the following finite-sample guarantee in the simplifying policy gradient update. This is also the more detailed statement of theorem 4.2 in the main paper.

Proposition A.6 (Finite-sample Guarantee under Noisy Gradient Updates and Projection Error). At each iteration, we perform policy gradient estimate using  m = O( (σ2+AHBβ) log T kδβ2 )trajectories

and use DAgger algorithm to collect  M = O( log T/δϵ2 ) + O( Hϵ )roll-outs. Setting the learning rate


with probability at least  1 − δ.

Consequently, we also have the following regret bound:


holds with probability at least  1 − δ, where again H is the task horizon, A is the cardinality of action space, and k is the dimension of function class  Π’s parameterization.

Proof. (For both proposition (A.6) and (A.5)). The results follow by taking the inequality from equation (19), and by solving for  ϵand  βexplicitly in terms of relevant quantities. Based on the specification of M and m, we obtain the necessary precision for each round of learning in terms of number of trajectories:


Setting the learning rate  η =�


The regret bound depends on the variance  σ2of the policy gradient estimates. It is well-known that vanilla policy gradient updates suffer from high variance. We instead use functional regularization technique, based on CORE-RL, in the practical implementation of our algorithm. The CORE-RL subroutine has been demonstrated to reduce the variance in policy gradient updates [19].

A.4 Defining a consistent approximation of  ∇HJ(π)- Proof of Proposition 4.3

We are using the notion of Fréchet derivative to define gradient of differentiable functional. Note that while Gateaux derivative can also be utilized, Fréchet derivative ensures continuity of the gradient operator that would be useful for our analysis.

Definition A.6 (Fréchet gradient). A bounded linear operator  ∇ : H �→ His called Fréchet functional gradient of  J at h ∈ H if lim∥g∥→0J(h+g)−J(h)−⟨∇J(h),g⟩∥g∥ = 0

We make the following assumption about H and F. One interpretation of this assumption is that the space of policies  Π and Fthat we consider have the property that a programmatic policy  π ∈ Π canbe well-approximated by a large space of neural policies  f ∈ F.

Assumption 1. J is Fréchet differentiable on H. J is also differentiable on the restricted subspace F. And F is dense in H (i.e., the closure F = H)

It is then clear that  ∀ f ∈ Fthe Fréchet gradient  ∇FJ(f), restricted to the subspace F is equal to the gradient of f in the ambient space H (since Fréchet gradient is unique). In general, given  π ∈ Π andf ∈ F, π + fis not necessarily in F. However, the restricted gradient on subspace F of  J(π + f)can be defined asymptotically.

Proposition A.7. Fixing a policy  π ∈ Π, define a sequence of policies  fk ∈ F, k = 1, 2, . . .that converges to  π: limk→∞ ∥fk − g∥ = 0, we then have  limk→∞ ∥∇FJ(fk) − ∇HJ(π)∥∗ = 0

Proof. Since Fréchet derivative is a continuous linear operator, we have limk→∞ ∥∇HJ(fk) − ∇HJ(π)∥∗ = 0. By the reasoning above, for  f ∈ F, the gradient ∇FJ(f)defined via restriction to the space F does not change compared to  ∇HJ(f), the gradient defined over the ambient space H. Thus we also have  limk→∞ ∥∇FJ(fk) − ∇HJ(π)∥∗ = 0. By the same argument, we also have that for any given  π ∈ Πand  f ∈ F, even if  π + f ̸∈ F, the gradient  ∇FJ(π + f)with respect to the F can be approximated similarly.

Note that we are not assuming  J(π)to be differentiable when restricting to the policy subspace  Π.

A.5 Theoretical motivation for Algorithm 2 - Proof of Proposition 4.4 and 4.5

We consider the case where  Πis not differentiable by parameterization. Note that this does not preclude  J(π) for π ∈ Πto be differentiable in the non-parametric function space. Two complications arise compared to our previous approximate mirror descent procedure. First, for each  π ∈ Π, estimating the gradient  ∇J(π)(which may not exist under certain parameterization, per section 4.3) can become much more difficult. Second, the update rule  ∇R(π) − ∇FJ(π)may not be in the dual space of F, as in the simple case where  Π ⊂ F, thus making direct gradient update in the F space inappropriate.

Assumption 2. J is convex in H.

By convexity of J in H, sub-gradients  ∂J(h)exists for all  h ∈ H. In particular,  ∂J(π)exists for all π ∈ Π. Note that  ∂J(π)reflects sub-gradient of  πwith respect to the ambient policy space H.

We will make use of the following equivalent perspective to mirror descent[10], which consists of two-step process for each iteration t

1. Solve for ht+1 = argminh∈H η⟨∂J(πt), h⟩ + DR(h, πt)

2. Solve for  πt+1 = argminπ∈Π DR(π, ht+1)

We will show how this version of the algorithm motivates our main algorithm. Consider step 1 of the main loop of PROPEL, where given a fixed  π ∈ Π, the optimization problem within H is


Due to convexity of H and the objective, problem (OBJECTIVE_1) is equivalent to:


where  τdepends on  η. Since  πis fixed, this optimization problem can be relaxed by choosing λ ∈ [0, 1], and a set of candidate policies  h = π + λf, for all  f ∈ F, such that  DR(h, π) ≤ τis satisfied (Selection of  λis possible with bounded spaces). Since this constraint set is potentially a restricted set compared to the space of policies satisfying inequality (22), the optimization problem (20) is relaxed into:


Due to convexity property of J, we have


The original problem OBJECTIVE_1 is thus upper bounded by:


Thus, a relaxed version of original optimization problem OBJECTIVE_1 can be obtained by miniziming  J(π + λf)over  f ∈ F(note that  πis fixed). This naturally motivates using functional regularization technique, such as CORE-RL algorithm [19], to update the parameters of differentiable function f via policy gradient descent update:


where the gradient of J is taken with respect to the parameters of f (neural networks). This is exactly the update step in algorithm 2 (also similar to iterative updte of CORE-RL algorithm), where the neural network policy is regularized by a prior controller  π.

Statement and Proof of Proposition 4.5

Proposition A.8 (Regret bound for the relaxed optimization objective). Assuming J(h) is L-strongly smooth over  H, i.e., ∇HJ(h) is L-Lipschitz continuous, approximating UPDATEH by UPDATEF perAlg. 2 leads to the expected regret bound:  E�1T�Tt=1 J(πt)�− J(π∗) = O�λσ�

Proof. Instead of focusing on the bias of the gradient estimate  ∇HJ(π), we will shift our focus on the alternative proximal formulation of mirror descent, under optimization and projection errors. In particular, at each iteration t, let  h∗t+1 = argminh∈H η⟨∇J(πt), h⟩ + DR(h, πt)and let the optimization error be defined as  βtwhere  ∇R(ht+1) = ∇R(h∗t+1) + βt. Note here that this is different from (but related to) the notion of bias from gradient estimate of  ∇J(π)used in theorem 4.1 and theorem A.2. The projection error from imitation learning procedure is defined similarly to theorem 4.1:  π∗t+1 = argminπ∈Π DR(π, ht+1)is the true projection, and��πt+1 − π∗t+1�� ≤ ϵ.

We start with similar bounding steps to the proof of theorem 4.1:


The bound on projection error is identical to theorem A.2:



��π∗t+1 − πt��2 + ⟨∇R(h∗t+1), π∗t+1 − πt⟩ + ⟨βt, π∗t+1 − πt⟩= −η⟨∇JH(πt), π∗t+1 − πt⟩ − α2


Note here that the gradient  ∇HJ(πt)is not the result of estimation. Combining equations (25), (26), (27), (28), we have:


Next, we want to bound  βt. Choose regularizer R to be 12 ∥·∥2(consistent with the pseudocode in algorithm 2). We have that:


which is equivalent to:


Let  f ∗t+1 = argminf∈F η⟨∇J(πt), f⟩ + 12 ∥f∥2. Taking the gradient over f, we can see that  f ∗t+1 =−η∇J(πt). Let ft+1be the minimizer of  minf∈F J(πt + λf). We then have  h∗t+1 = πt + f ∗t+1 andht+1 = π + λft+1. Thus βt = ht+1 − h∗t+1 = ft+1 − f ∗t+1.

On one hand, we have


due to optimality of  ft+1and strong smoothness property of J. On the other hand, since J is convex,


Combine with the inequality above, and subtract  J(πt)from both sides, and using the relationship f ∗t+1 = −η∇J(πt), we have that:


Since this is true  ∀ω, rearrange and choose  ωsuch that ωη − Lω22 = − λ2η, namely  ω = 1−√1−ληLLη, and complete the square, we can establish the bound that:


for B the upperbound on  ∥ft+1∥. We thus have  ∥βt∥ = O(η(λL)2). Plugging the result from equation 30 into RHS of equation 29, we have:


Since J is convex in  H, we have J(πt) − J(π) ≤ ⟨∇J(πt), πt − π⟩. Similar to theorem 4.1, setting


Note that unlike regret bound from theorem 4.1 under general bias, variance of gradient estimate and projection error,  σ2here explicitly refers to the bound on neural-network based policy gradient variance. The variance reduction of  λσ, at the expense of some bias, was also similarly noted in a recent functional regularization technique for policy gradient [19].


We generate controllers for cars in The Open Racing Car Simulator (TORCS) [59]. In its full generality TORCS provides a rich environment with input from up to 89 sensors, and optionally the 3D graphic from a chosen camera angle in the race. The controllers have to decide the values of 5 parameters during game play, which correspond to the acceleration, brake, clutch, gear and steering of the car.

Apart from the immediate challenge of driving the car on the track, controllers also have to make race-level strategy decisions, like making pit-stops for fuel. A lower level of complexity is provided in the Practice Mode setting of TORCS. In this mode all race-level strategies are removed. Currently, so far as we know, state-of-the-art DRL models are capable of racing only in Practice Mode, and this is also the environment that we use. Here we consider the input from 29 sensors, and decide values for the acceleration, steering, and braking actions.

We chose a suite of tracks that provide varying levels of difficulty for the learning algorithms. In particular, for the tracks Ruudskogen and Alpine-2, the DDPG agent is unable to reliably learn a policy that would complete a lap. We perform the experiments with twenty-five random seeds and report the median lap time over these twenty-five trials. However we note that the TORCS simulator is not deterministic even for a fixed random seed. Since we model the environment as a Markov Decision Process, this non-determinism is consistent with our problem statement.

For our Deep Reinforcement Learning agents we used standard open source implementations (with pre-tuned hyper-parameters) for the relevant domain.

All experiments were conducted on standard workstation with a 2.5 GHz Intel Core i7 CPU and a GTX 1080 Ti GPU card.

The code for the TORCS experiments can be found at: https://bitbucket.org/averma8053/propel

In Table 3 we show the lap time performance and crash ratios of PROPEL agents initialized with neural policies obtained via DDPG. As discussed in Section 5, DDPG often exhibits high variance across trials and this adversely affects the performance of the PROPEL agents when they are initialized via DDPG. In Table 4 we show generalization results for the PROPELTREE agent. As noted in Section 5, the generalization results for PROPELTREE are in between those of DDPG and PROPELPROG.

Verified Smoothness Property. For the program given in Figure 2 we proved using symbolic verification techniques, that  ∀k, �k+5i=k ∥peek(s[RPM], i + 1) − peek(s[RPM], i)∥ < 0.003 =⇒∥peek(a[Accel], k + 1) − peek(a[Accel], k)∥ < 0.63. Here the function peek(., i) takes in a history/sequence of sensor or action values and returns the value at position i, . Intuitively, the above logical implication means that if the sum of the consecutive differences of the last six RPM sensor values is less than 0.003, then the acceleration actions calculated at the last and penultimate step will not differ by more than 0.63.

Table 3: Performance results in TORCS of PROPEL agents initialized with neural policies obtained via DDPG, over 25 random seeds. Each entry is formatted as Lap-time / Crash-ratio, reporting median lap time in seconds over all the seeds (lower is better) and ratio of seeds that result in crashes (lower is better). A lap time of CR indicates the agent crashed and could not complete a lap for more than half the seeds.


Table 4: Generalization results in TORCS for PROPELTREE, where rows are training and columns are testing tracks. Each entry is formatted as PROPELPROG / DDPG, and the number reported is the median lap time in seconds over all the seeds (lower is better). CR indicates the agent crashed and could not complete a lap for more than half the seeds.


Table 5: Performance results in Classic Control problems. Higher scores are better.


B.2 Classic Control

We present results from two classic control problems, Mountain-Car (with continuous actions) and Pendulum, in Table 5. We use the OpenAI Gym implementations of these environments. More information about these environments can be found at the links: MountainCar and Pendulum.

In Mountain-Car the goal is to drive an under-powered car up the side of a mountain in as few time-steps as possible. In Pendulum, the goal is to swing a pendulum up so that it stays upright. In both the environments an episode terminates after a maximum of 200 time-steps.

In Table 5 we report the mean and standard deviation, over twenty-five random seeds, of the average scores over 100 episodes for the listed agents and environments. In Figure 6 and Figure 7 we show the improvements made over the prior by the PROPELPROG agent in MountainCar and Pendulum respectively, with each iteration of the PROPEL algorithm.


Figure 6: Score improvements in the MountainCar environment over iterations of PRO- PELPROG.

Figure 7: Score improvements in the Pendulum environment over iterations of PROPEL-PROG.

Designed for Accessibility and to further Open Science