# Patent

## Title
Control flow graph flattening device and method

## Abstract

Control Flow Graph flattening of a function of software code including a plurality of basic blocks having an address and at least one instruction is described. A processor uses a transition function to obtain an address in a jump table of at least one diffusion function coefficient of a subsequent block of the plurality of basic blocks. Thereafter, the processor generates and inserts branch functions, assigning a value to the label of each of the plurality of basic blocks and then obtains coefficients for each of the plurality of basic blocks.

## Background
This application claims the benefit, under 35 U.S.C. §119 of European Patent Application No. 14306695.9, filed Oct. 24, 2014 and European Patent Application No. 15305231.1, filed Feb. 16, 2015.
TECHNICAL FIELD
The present principles relate generally to computer systems and in particular to control flow graph (CFG) flattening.
BACKGROUND
This section is intended to introduce the reader to various aspects of art, which may be related to various aspects of the present principles that are described and/or claimed below. This discussion is believed to be helpful in providing the reader with background information to facilitate a better understanding of the various aspects of the present principles. Accordingly, it should be understood that these statements are to be read in this light, and not as admissions of prior art.
Control flow graph (CFG) flattening is a software obfuscation technique used to make reverse engineering of a function difficult. The initial CFG of a function, made of basic blocks (BB) and jumps between them, is transformed to a functionally equivalent flattened graph. In the flattened graph, a central basic block, the dispatcher (D), leads to all other BBs, which all call back to the dispatcher. Wang first described this technique in a PhD thesis “A Security Architecture for Survivable Systems”, Department of Computer Science, University of Virginia, and in C. Wang, J. Davidson, J. Hill & J. Knight, “Protection of software-based survivability mechanisms,” In Dependable Systems and Networks 2001, DSN 2001, International Conference on. IEEE. 
FIG. 1 illustrates an exemplary function f and its corresponding CFG comprising four BBs, A-D. FIG. 2 illustrates a flattened CFG corresponding to the CFG in FIG. 1. It can be seen that the dispatcher executes a switch function S between the execution of BBs.
The control flow during execution is dependent on a dispatcher variable X. This variable is instantiated locally to a label Xi in each basic block BBi. At runtime, each basic block BBi is assigned a new label Xi+1 thereby indicating which next basic block to be executed next. A switch block then uses the dispatcher variable X to jump indirectly, via a jump table, to the intended successive basic block.
Every basic block BBi embeds instructions to compute label Xi+1 from label Xi. Caeppaert calls this a branch function [see J. Cappaert, A General Model for Hiding Control Flow, Proceedings of the tenth annual ACM workshop on DRM, ACM 2010]. In this description the branch function is denoted Bi(Xi).
As an example, in basic block A, the statement {x=i<j? 1:2} can be expressed generally as the result of a branch function BA( ): BA(XA)=BooleanFunction(inputs)*a+b. More particularly, BA(Xi)=(i<j?)*1+(i>=j?)*2. (i<j?) is a Boolean function that returns 1 or 0, so BA( ) returns 1 or 2 depending on the inputs, i and j. A branch function can be generalized to any conditional statement having more than two term conditions.
The dispatcher uses another function, the transition function F( ). In FIG. 2, this function is responsible for converting the label Xi, to an address PAi of the next BB to execute. When flattening a CFG at source code level (e.g. for a program in C language), the transition function F( ) may be implicit. For example, the switch case in the dispatcher is converted to a sequence of instructions by the compiler. However, when flattening the CFG for a low level assembly language (e.g. X86), the transition function normally has to be written expressly.
FIG. 3 illustrates a generic sequence diagram for a flattened CFG. A basic block BBi comprises two parts: a first part 30 executed upon entry to the BB and a second part 32 that is executed before leaving the BB. The second part 32 calculates the next label Xi+1 that is used by a dispatcher 34 to compute the address of the next BB and then jumps to this address. The bold border of the first part 30 indicates that this part is part of the original, non-flattened, CFG.
Basic implementations of flattened CFGs are vulnerable to static analysis attacks. By analysing each basic block and the constants used by the branch function B( ), an attacker is able to gain knowledge of dependencies between basic blocks, and then rebuild the original CFG. For example, in block 1 in FIG. 2, the conditional statement (x=i<j? 1:2) permits to deduce a relation between basic block A and basic blocks B and C.
A more generalized technique proposed by Cappaert in the previously mentioned paper consists in expressing the label Xi+1 relatively to the label Xi: Bi(Xi+1)=Xi+ai*BooleanFunction(inputs)+bi. The meaning of (ai,bi) has slightly evolved, it no longer directly expresses a label value, but relative differences between labels. However, Cappaert observes that this representation is still vulnerable to local analysis. The differences between labels restrict the set of possibilities for Xi, and thus leak information to the attacker, making brute force attacks feasible.
More specifically, the attacker has access to the code of the transition function, and thus knows the correspondence between labels Xi and addresses PAi of the basic blocks. To deduce the label of a basic block is straightforward: the attacker just has to compare the current address of the basic block and the information in the dispatcher. So for a given BBi, the attacker can guess what the current label Xi is, and then propagate the constant values (ai,bi) backward and forward to deduce the predecessor and the successor of each basic block.
It can thus be seen that the existence of a logical inference chain between embedded constants in the branch function, the logical labels, and the address of the basic block causes a problem. Existing protections find a way of breaking this chain of logical inferences.
A first solution was proposed by Wang in the previously mentioned PhD thesis. The solution focuses on the obfuscation of the embedded constants. These are not changed, but simply expressed in term of formulae, using a global array. For this purpose Wang defines a global array g_array, containing fakes values and valid values. If, for example, g_array[10]=1 and g_array[20]=1, then the constant Ci=2 can be expressed by mathematical operation on elements of the array: Ci=(g_array[10]+g_array[20]) mod N. However, the solution is hardly more robust as the array is static and constant. By obtaining access to the array through reverse engineering, an attacker can resolve the different formulae.
To counter this attack, Wang further proposes to periodically reshuffle and redistribute the content of g_array using a redistribution function that is invoked outside of the flattened function. Ideally, this solution could be secure if the redistribution is not done externally, but internally in the flattened function during transition between basic blocks. However, Wang does not give any clues on how to achieve this. Moreover, there is a further difficulty. In a graph, multiple paths can lead to the same basic block. If for example both basic blocks A and B lead to C, then the content of the array is not deterministic in C.
A second solution was proposed by Jung Ge et al. [See Jung Ge, Soma Chaudhuri, Akhilesh Tyagi: Control Flow Based Obfuscation, in Proceedings of DRM 05, the 5th ACM Workshop on Digital Rights Management]. According to this solution called CFI-hiding scheme, based upon Aucsmith's algorithm, sets of basic blocks are shuffled, XORed together, and dynamically deciphered at execution by a monitoring process called Self-modifying scheme. At any time the only set of basic block that is in clear text is the one that is already executed. [See D. Aucsmith: Tamper Resistant Software: An Implementation, in R. J. Anderson, editor, Information Hiding, volume 1174 of Lecture Notes in Computer Science].
A third solution was provided by Apple in patent application US 2013/0232323 that describes a solution that obfuscates the control flow in a different way. All calls are replaced by indirect jumps using a jump table. The basic blocks are not encrypted, but the jump table is; only a portion of the jump table is in the clear at a given time.
A common drawback with self-modifying schemes like the one by Jung Ge et al. is that the program requires write permission on itself. Such protection techniques are sometimes not allowed on some application platforms, like the iOS appstore, for which these techniques are assimilated to polymorphic viruses. Thus a protected application would not pass the security verifiers in charge of verifying the safety of the application before publishing it. Another problem with self-modifying techniques is their vulnerability to replay attacks. The randomization induced by the protection can be overcome if the attacker has means to restart execution from a backup initial state.
A fourth solution was provided by Cappaert in the previously mentioned paper. The solution uses a one-way function, typically a hash function applied on label values Xi when entering the dispatcher. The resulting value is then matched against a series of if . . . then clauses to determine the next address.
By static analysis of the dispatcher in depth, an attacker obtains knowledge of all the possible images of label values by the hash function and of all addresses referenced. But as the hash function is not invertible, the antecedent label values cannot be obtained with certainty. By coding the label values with enough bits (for example a 128-bit random number), the input domain of the hash function cannot be enumerated by brute force. The attacker then cannot predict which antecedent Xi label values leads to a given Zi value, and, thus, the corresponding antecedent address values.
It will be appreciated that this stops logical backward inference and makes forward inference more complex. In other words, the attacker cannot statically predict the predecessor of a basic block, and the constants used in the branching function have only a meaning for logical labels, and do not express indirectly constraints to basic block addresses. Hence, differences values between labels are no longer exploitable.
Unfortunately, while this solution provides good security, it has a serious drawback in terms of performance. The drawback stems from the dispatcher. In practice, the compiler replaces the dispatcher's switch statement by a cascade of if . . . else condition statements. Given a CFG that has N basic blocks, the switch case will result in the cascade of N more conditional statements, so N more conditional basic blocks. A first consequence is that the number of total basic blocks in the CFG is multiplied by two (twice more times jump instructions). A second consequence is that each time the processor executes the dispatcher, it faces a cascade of N conditional statements, which results in the cost of potentially 2*N additional instructions. Moreover, this technique adds many comparisons and jump instructions to the execution, and also stresses the branch prediction much. For each missed branch prediction the processor started to execute instructions that then are not executed and the processor then has to ‘undo’ everything before continuing execution at the right address.
Indeed, when using a flattened CFG based upon switch case, execution time is proportional to the square of the number of basic blocks (O(n2)). The time penalty is so great for large functions that it is difficult to envisage the solution to protect a major part of the protected binary. As a result, to keep a good execution performance, the coverage of the protection has to be limited to only some sensitive portions of the code, which naturally decreases the security. The poor performance result is a direct consequence of the selected hash function; the outputs of such functions are discontinuous values that cannot be mapped directly to an indexed table.
It will thus be appreciated that there is a need for a CFG flattening solution addresses these and other problems. The present principles can provide such a solution.
SUMMARY
In a first aspect, the principles are directed to a method for control flow graph flattening of a function of software code, the function comprising a plurality of basic blocks each having an address and at least one instruction, the method, performed by a processor of a device, comprising assigning a value to at least one diffusion function coefficient of a diffusion function of each of the plurality of basic blocks; calculating changes for labels and diffusion functions indices upon transitions between the plurality of basic blocks, the changes corresponding to modifications made by a linear first branch function and a linear second branch function; creating a dispatcher basic block comprising instructions to look up an address in a jump table and to jump to the address; creating the jump table associating the label of each basic block with the address of the basic block; creating a coefficient array comprising diffusion function coefficients; replacing a jump terminal instruction by a jump to the dispatcher basic block in each of the plurality of basic blocks; creating and inserting into each of the plurality of basic blocks the first branch function and the second branch function, the first branch function taking as input a result of a diffusion function in the present basic block and outputting the label of the subsequent basic block to execute, the second branch function taking as input a diffusion function index of the present basic block and outputting a diffusion function index of the subsequent basic block based; creating and inserting a transition function into the dispatcher basic block, the transition function obtaining the address in the jump table based on at least the label of the subsequent basic block; and creating and inserting diffusion functions in the plurality of basic blocks, a diffusion function inserted in the present basic block calculating a value based on at least one diffusion function coefficient dependant on the diffusion function index and the label of the present basic block.
Various embodiments of the first aspect comprise that:

    
    
        the transition function obtains the address in the jump table based on also at least one diffusion function coefficient of the subsequent basic block;
        before creating and inserting the branch functions, assigning a value to the label of each of the plurality of basic blocks;
        obtaining coefficients for each of the plurality of basic blocks.
    
    


In a second aspect, the principles are directed to a device for control flow graph flattening of a function of software code, the function comprising a plurality of basic blocks each having an address and at least one instruction, the device comprising a hardware processor configured to assign a value to at least one diffusion function coefficient of a diffusion function of each of the plurality of basic blocks; calculate changes for labels and diffusion functions indices upon transitions between the plurality of basic blocks, the changes corresponding to modifications made by a linear first branch function and a linear second branch function; create a dispatcher basic block comprising instructions to look up an address in a jump table and to jump to the address; create the jump table associating a label of each basic block with the address of the basic block; create a coefficient array comprising diffusion function coefficients; replace a jump terminal instruction by a jump to the dispatcher basic block in each of the plurality of basic blocks; create and insert into each of the plurality of basic blocks the first branch function and the second branch function, the first branch function taking as input a result of a diffusion function in the present basic block and outputting the label of the subsequent basic block to execute, the second branch function taking as input a diffusion function index of the present basic block and outputting a diffusion function index of the subsequent basic block based; create and insert a transition function into the dispatcher basic block, the transition function obtaining the address in the jump table based on at least the label of the subsequent basic block; and create and insert diffusion functions in the plurality of basic blocks, a diffusion function inserted in the present basic block calculating a value based on at least one diffusion function coefficient dependant on the diffusion function index and the label of the present basic block; and an interface configured to output the function after control flow graph flattening.
Various embodiments of the second aspect comprise that the processor is further configured to:

    
    
        make the transition function obtain the address in the jump table based on also at least one diffusion function coefficient of the subsequent basic block;
        assign a value to the label (X) of each of the plurality of basic blocks, before creating and inserting the branch functions;
        obtain coefficients for each of the plurality of basic blocks.
    
    





BRIEF DESCRIPTION OF DRAWINGS
Preferred features of the present principles will now be described, by way of non-limiting example, with reference to the accompanying drawings, in which
FIG. 1 illustrates an exemplary function and its corresponding CFG;
FIG. 2 illustrates a flattened CFG corresponding to the CFG in FIG. 1;
FIG. 3 illustrates a generic sequence diagram for a flattened CFG;
FIG. 4 illustrates a computer program with a flattened CFG;
FIG. 5 illustrates a method for CFG flattening according to the present principles;
FIG. 6 illustrates an exemplary unflattened CFG; and
FIG. 7 illustrates an exemplary CFG flattening device.



DESCRIPTION OF EMBODIMENTS
In its current implementation, Wang's solution does not resist well against static analysis and leak of information. The constant coefficients used in every branch function of each BB should not be predictable by static analysis, so they should not be assigned constant values, even obfuscated from outside call sites. One aspect of the present principles is to overcome limitations of Wang's solution by providing a stronger obfuscation of coefficients.
The present principles rely on the use of a live variable X and a diffusion function D, instantiated randomly at compilation time for each BB of the code to be protected; where necessary for the sake of clarity, for BBi the variable is denoted Xi and the diffusion function Di. The variable and the diffusion function can be made uncorrelated and unpredictable by static analysis and are:

    
    
        a current logical label Xi in the range {1 . . . n}, where n is the number of BBs of the code to protect; and
        a current diffusion function D( )ε{DP1( ) . . . DPm( )} that modifies the values of the current logical label Xi. In practice, the diffusion function is identified using a diffusion function index Pi that is updated after use in the basic blocks.
    
    


It is assumed that X0 and P0 are retrieved securely or passed as parameters in any suitable prior art manner; how this is done is beyond the scope of the present principles. X and P will be used without indices except when this may cause confusion, e.g. when referring to both a current value and a future value.
During CFG flattening to protect a binary, a flattening tool inserts the following in the intermediate representation to obtain the flattened CFG, illustrated in FIG. 4 for one generic BB and a dispatcher:

    
    
        a jump table;
        a dispatcher basic block configured to determine the next basic block to execute;
        a transition function F(X): {1 . . . n}→N, called by the dispatcher, that translates, using the jump table, the logical label X to an address PA of the BB to which execution will jump is going to jump;
        a diffusion function DPi( ), to be executed when entering BBi. DPi( ) outputs a local label Z used to compute a new logical label X. DPi( ) does not need to be one-way. The only preferred restriction on the choice of DPi is that it ensures some diffusion property in {1 . . . n} and does not decrease entropy.
        a first branch function B1i( ) inlined in each BBi that computes the next logical label X (in {1 . . . N}) depending on the local label Z, one or more local variables defined in the basic block and at least one conditional statement Yi( ); and
        a second branch function B2i( ) inlined in each BBi that computes the next diffusion function DPi+1, depending on the current diffusion function DPi and the conditional statement Yi( ) or, put another way, computes the next diffusion function index Pi+1 based on the current function index Pi.
    
    


The first branch functions and the second branch functions are advantageously different between at least some, preferably all of the basic blocks.
As mentioned, FIG. 4 illustrates a computer program with a flattened CFG, the program comprising a dispatcher 48 and a basic block comprising an entry part 40, an instruction part 42 comprising the original instructions of the BB, a first branch function part 44 and a second branch function part 46. FIG. 4 also illustrates a jump table 45 and a diffusion parameter table 47. It will be understood that the program comprises a plurality of basic blocks, but that only one is illustrated for reasons of clarity. In addition, the modulus function has been left out of some formulae so as not to burden the description; the skilled person will know when the modulus is needed.
During execution, a jump is made from the dispatcher 48 to the entry part 40 of the BB and the local label Z is calculated using the current diffusion function DPi operating on the current logical label X. Then the instruction part 42 is executed and the conditional statement Yi( ) is evaluated, advantageously using one or more Boolean functions. Execution moves on to the first branch function part 44 preferably comprising at least one constant, in the example two constants ai, bi. The first branch function B1i is then evaluated taken the local label Z as input and preferably using the at least one constant and the evaluated conditional statement Yi( ); in the example: X←B1i(Z)=Z*ai*Yi( )+bi*Z mod (2n). Then the second branch function part 46 is executed in which the next diffusion function Di+1 (i.e., the next diffusion function index Pi+1) is calculated using the second branch function B2i( ) preferably using at least one constant taking as input the current diffusion function Di; in the example: Pi←B2i(Pi)=Pi+ciYi( )+di, where ci and di are constants. Although not illustrated, the second branch function part 46 finishes by a jump to the dispatcher 48.
In an advantageous embodiment (illustrated in FIG. 4), the transition function F( ) further depends on the at least a diffusion parameter of the diffusion function Di such as the current value of the diffusion function index P. The jump table is then obfuscated with the diffusion function index Pi using any suitable obfuscation method such as, but not limited to, encryption of table cells with the diffusion function index Pi or use of the diffusion function index Pi as a seek offset in the table.
A role of the diffusion function D( ) is to ensure that any constants in the first branch function B1i( ) are relative to the current instantiation of DPi( ). As a result, the constants can have a different meaning in each BBi and a global static analysis of all constants within the program do not reveal any information.
There are several appropriate diffusion functions that can be used, such as for example any of the following functions or any combination of them:
1. an affine function in Z or Z/nZ;
2. a permutation, a rotation or a S-BOX; and
3. a linear feedback shift register (LFSR).
The cost in performance and memory footprint is different for each function:
1. a table of affine functions parameter for all functions;
2. a list of permutation tables, rotation tables or S-BOXes; and
3. a list of polynomials in GF(2n).
An affine function is preferred as it represents a good trade-off between performance and memory footprint.
The second branch function B2i( ) can be any suitable function from Z to Z, or Z to Z/nZ (e.g. modulo n). As there are no security constraints on it, a simple function can be used, such as the one used in the example: Pi←B2i(Pi)=Pi+ciYi( )+di.
The diffusion function DPi( ) and the branch function B1i( ) should preferably be resistant to a comparison attack. The following functions are examples of such functions:

DPi(X):Z→Z:DPi(X)=2*coeffi*X+1,

where the coefficient, coeffi, is obtained from the diffusion parameter table 47; and

B1i( ):Z→Z/nZ:B1i(Z)=ai*Z·Yi( )+bi*Z mod 2n 

Assume that BBA (label Xa) has two successors BBb and BBc, with labels Xb and Xc. B1a(Xa) can return Xb or Xc depending on whether Yi( ) returns 1 or 0. In this example, Zi=DPi(Xa), Xb=(ai*Zi*1+bi·Zi) mod 2n=>Xb=(ai*Zi+bi*Zi) mod 2n and Xc=(ai*Zi*0+bi*Zi) mod 2n=>Xc=bi*Zi mod 2n. Subtracting the equations for Xb and Xc gives Xc−Xb=ai*Zi mod 2n. Since ai*Zi is not constant but depends on the variables DPi and Xi, a comparison attach does not succeed.
The images Zi of the diffusion function serves as input to the first branch function B1i. In other words, the codomain of DPi( ) is the input domain Din for B1( ). The codomain Dout of B1( ) regroups random label values. The choice of a diffusion function DPi should then satisfy a reachability property for B1( ): for any values (x,y) from (Din,Dout), there should be a coefficient for B1( ) such that y is the image of x by B1( ). There are many such functions. Using two constants, it is always possible to satisfy the requirement. For example, if Xi+1=B1i(Zi)=(ai*Zi·Yi( )+bi*Zi) mod 2n, there must be at least one solution (ai,bi) that satisfies the equation, whatever the values of Zi and Xi+1 are since Zi is odd.
Another way to write B1( ) is to use the Boolean function invYi( ) which returns the inverse result of Yi( ), e.g. invYi( )=0 if Yi( )=1. This makes it possible to write:




  
    
      
        
          X
          ←
          
            B
            ⁢
            
              
            
            ⁢
            
              1
              i
            
            ⁢
            
              (
              Z
              )
            
          
        
        =
        
        ⁢
        
          
            (
            
              
                
                  a
                  i
                  ′
                
                *
                Z
                *
                
                  
                    Y
                    i
                  
                  ⁡
                  
                    (
                    
                      
                    
                    )
                  
                
              
              +
              
                
                  b
                  i
                  ′
                
                *
                Z
                *
                
                  
                    invY
                    i
                  
                  ⁡
                  
                    (
                    
                      
                    
                    )
                  
                
              
            
            )
          
          ⁢
          mod
          ⁢
          
            
          
          ⁢
          
            2
            n
          
        
      
    
  
  
    
      
        
          =
          
          ⁢
          
            
              
                a
                i
                ′
              
              *
              Z
              ⁢
              
                
              
              ⁢
              mod
              ⁢
              
                
              
              ⁢
              
                2
                n
              
              ⁢
              
                
              
              ⁢
              if
              ⁢
              
                
              
              ⁢
              
                
                  Y
                  i
                
                ⁡
                
                  (
                  
                    
                  
                  )
                
              
            
            =
            1
          
        
        ,
      
    
  
  
    
      
        =
        
        ⁢
        
          
            
              b
              i
              ′
            
            *
            Z
            ⁢
            
              
            
            ⁢
            mod
            ⁢
            
              
            
            ⁢
            
              2
              n
            
            ⁢
            
              
            
            ⁢
            if
            ⁢
            
              
            
            ⁢
            
              
                Y
                i
              
              ⁡
              
                (
                
                  
                
                )
              
            
          
          =
          0.
        
      
    
  




Thus, the equation is satisfied if there is an “a” such that: a*Z mod 2n=Xj ∀Z,Xj. This equation has a solution if Z and 2n are coprime, in which case Z is a generator of Z/2nZ. An easy way to insure this property is to define the diffusion function so that it returns an odd number. If Z is odd, e.g Z=2*b+1, then Z and 2n are coprime. 2 is the unique divisor of 2n, and 2 does not divide 2*b+1.

It is thus possible to choose:

DPi:Z→Z:DPi(X)=2*coeffPi*X+1.

A more generic solution is: DPi: Z→Z: Di(X)=G(X) & 0xFFFFFFFC+0x1, where G is a function in Z→Z. F can for example be a polynomial function, a permutation in {1 . . . N}→{1 . . . N}, or a rotation. It is noted that such a generic function DPi is not necessary injective. This does not impact the solution of the present disclosure as it does not require unicity of DPi values.
The diffusion functions can also be set to be almost any kind of function through the use of pointers. In this case, the coefficients are pointers to functions, an example being DPi(X)=2*(*coeffPi)(X)+1, where (*coeffPi)(X) indicates, with reference to the diffusion parameter table 47, that the function at address coeffi takes X as input. The expression “diffusion parameter” is intended to cover both a value obtained from at least one of the coefficients of the diffusion function and an address of a diffusion function.
In an embodiment, for each BBi, DPi(X)=2*coeffPi*X+1. This can be implemented by declaring the global array 47 g_DiffusionCoeff[ ] that comprises the coefficients coeffi for iε{0 . . . M−1}. The embodiment uses an integer live variable P such as Pi ε{0,M−1}. Pi acts as an index in the global array 47 g_DiffusionCoeff[ ] to retrieve the current coefficient coeffPi. Thus the following function is inlined in basic block DPi: Zi=2*g_DiffusionCoeff[Pi]*X+1.
The skilled person will appreciate that FIG. 4 is merely an example and that various modifications may be made. For example, the constants can be defined anywhere within the BB as long as this is done before they are used, there can be a different number of constants used in the formulae, the formulae used to obtain the derived values may be different, the first branch function B1i( ) and the second branch function B2i( ) may be different, and so on.
FIG. 5 illustrates a method for CFG flattening according to the present principles. The method is advantageously performed by a CFG flattening device and takes as input a computer program with N basic blocks, where N is an integer greater than 1. The method can flatten the CFG for all the functions of the computer program or for a subset of these functions. An example for the unflattened CFG comprising basic blocks A-E, illustrated in FIG. 6, will be given for some steps.
In step 502, the label for each basic block is given or assigned a distinct value, preferably in {1 . . . N}. This can be done by shuffling the basic blocks using the Fisher-Yates shuffle algorithm or the like.

    
    
        BBA: XA=4
        BBB: XB=1
        BBC: XC=2
        BBD: XD=3
        BBE: XE=5
    
    


In step 504, the diffusion function index Pi for each basic block is given or assigned a value in {0 . . . M−1}. As there is no uniqueness constraint on Pi, it is possible to pick random values.
In the example, the diffusion function indices will be given abstract values P0-P2 (i.e. M=3) that in reality would correspond to integers.

    
    
        BBA: PA=P1
        BBB: PB=P2
        BBC: PC=P0
        BBD: PD=P2
        BBE: PE=P1
    
    


In step 506, the actions upon transitions between the basic blocks are calculated, i.e. the changes for the labels X and the diffusion function indices Pi. In other words, this corresponds to the modifications made by the first branch function B1i( ) and the second branch function B2i( )
In the example, the calculations for the diffusion function indices is performed modulo 3, since M=3 in the example.

BBA→BBB:XB=XA−3,PB=PA+1 mod 3

BBA→BBC:XC=XA−2,PC=PA+2 mod 3

BBC→BBD:XD=XC+1,PD=PC+2 mod 3

BBC→BBE:XE=XC+3,PE=PC+1 mod 3

In step 508, the dispatcher basic block is created. As already mentioned, the dispatcher comprises instructions to look up an address in the jump table and to jump to this address of the dispatcher basic block.
In step 510, the global jump table is created. The global jump table associates, for each basic block, the label X with the address of the basic block BBi.
In step 512, the array 47 of diffusion functions coefficients is generated. The array 47 of diffusion functions coefficients comprises coefficients used by the diffusion functions Di( ) and preferably also dummy values.
In step 514, the jump terminal instruction to another basic block is replaced by a jump to the dispatcher in each basic block of the original CFG.
In step 516, integer variables X and D and allocated in the stack of the target function of the original CFG and they are respectively initialized to X0 and D0. To do this, four intermediate language instructions are added in the first basic block of the function. Two integer allocation instructions are inserted to declare variables X and D into the stack and two affectation instructions are inserted to initialize these variables respectively to X0 and D0.
In step 518, the first branch functions B1i( ) and the second branch functions B2i( ) are created and inserted along with necessary constants into the basic blocks before the jump terminal instructions.
In step 520, the transition function F( ) is created and inserted into the dispatcher.
In step 522, the instructions for the diffusion functions Di( ) are created and inserted along with necessary constants into the basic blocks.
The constants for the branch functions B1i( ) and B2i( ) are inserted with the branch functions. In the example, these constants are {ai, bi} and {ci, di}. For each basic block BBi, the CFG flattening device knows the current values of Xi and Zi as well as the values Xi+1 and Zi+1 of its successors. Determining the coefficients can thus be done by resolving the equations Xi+1=B1i(Zi), DPi+1=B2i(DPi). The resolutions of the equations are trivial as these functions are linear and invertible.
As an example, imagine that basic block BB3 with label Xi=3 jumps to BB5: Xi=5 for k>2, and to BB6: Xi=6 otherwise. The conditional statement can be expressed by a Boolean function Yi(k)=1 if k>2, 0 else. Then it is possible to write:




  
    
      
        
          X
          
            i
            +
            1
          
        
        =
        
        ⁢
        
          
            B
            3
          
          ⁡
          
            (
            
              X
              i
            
            )
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            X
            i
          
          +
          
            
              
                Y
                i
              
              ⁡
              
                (
                k
                )
              
            
            ×
            
              (
              
                5
                -
                3
              
              )
            
          
          +
          
            
              (
              
                1
                -
                
                  
                    Y
                    i
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
              )
            
            ×
            
              (
              
                6
                -
                3
              
              )
            
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            X
            i
          
          +
          
            
              
                Y
                i
              
              ⁡
              
                (
                k
                )
              
            
            ×
            2
          
          +
          
            
              (
              
                1
                -
                
                  
                    Y
                    i
                  
                  ⁡
                  
                    (
                    k
                    )
                  
                
              
              )
            
            ×
            3
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            X
            i
          
          +
          
            
              
                Y
                i
              
              ⁡
              
                (
                k
                )
              
            
            ×
            
              (
              
                -
                1
              
              )
            
          
          +
          3
        
      
    
  




which means that the coefficients of branch function BB3( ) are ai=−1, bi=3. The coefficients for Di can be obtained in the same manner.

The different Xi can be obtained by enumerating the possible solutions and then validating one of them. Since a modulus function is used, there may be more than one solution, but it is certain that there is at least one function that may be found this way.
In the same manner, the CFG flattening device solves DPi+1=B2i(DPi) to obtain coefficients ci,di. The coefficients ci,di are inserted as local constants inside the basic block BBi and referenced by the second branch function B2i( ).
The branch functions B1i( ) and B2i( ) can be adapted for a basic block having 1 to K successors. For each successor BBj of BBi, there is a conditional statement that expresses the condition to branch from BBi to BBj, given some input. Let Yij( ) denote the Boolean function in BBi, iε{1 . . . N}, jε{1 . . . K}:




  
    
      
        
          
            Y
            ij
          
          ⁡
          
            (
            input
            )
          
        
        =
        
        ⁢
        
          
            1
            ⁢
            
              
            
            ⁢
            if
            ⁢
            
              
            
            ⁢
            
              (
              input
              )
            
          
          ==
          
            value
            j
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          0
          ⁢
          
            
          
          ⁢
          
            else
            .
          
        
      
    
  




Then Xi+1=B1i(Xi)=Xi+Σj=1K(Yij(input)*a′ij).

FIG. 7 illustrates an exemplary CFG flattening device 700 in which the present method is implemented. The device 700 comprises at least one hardware processing unit (“processor”) 710, memory 720, a user interface 730 for interacting with a user, and a communications interface 740 for communication with further devices over a connection such as the Internet. The skilled person will appreciate that the illustrated device is very simplified for reasons of clarity and that real devices in addition would comprise features such as power supplies and persistent storage.
The processor 710 is configured to obtain a computer program with an at least partly unflattened CFG, flatten the CFG of the computer program using the present principles (advantageously the method illustrated in FIG. 5), and output a computer program with an at least partly flattened CFG. The flattened computer program is advantageously stored on a non-transitory storage medium (not shown).
It will be appreciated that the flattening of the CFG according to the present principles can be made during compilation at the intermediate language level, but that it can also be made at source or binary level.
Moreover, the present principles can provide a generic and automatic solution that can work for any kind of CFG, whatever it looks like (for instance, with loop circuits, multiple basic block predecessors, multiple basic block successors, etc. . . . ).
Each feature disclosed in the description and (where appropriate) the claims and drawings may be provided independently or in any appropriate combination. Features described as being implemented in hardware may also be implemented in software, and vice versa. Reference numerals appearing in the claims are by way of illustration only and shall have no limiting effect on the scope of the claims.

## Claims

1. A method for control flow graph flattening of a function of software code, the function comprising a plurality of basic blocks each having an address and at least one instruction, the method, performed by a processor of a device, comprising:
assigning a value to at least one diffusion function coefficient of a diffusion function of each of the plurality of basic blocks;
calculating changes for labels and diffusion functions indices upon transitions between the plurality of basic blocks, the changes corresponding to modifications made by a linear first branch function and a linear second branch function;
creating a dispatcher basic block comprising instructions to look up an address in a jump table and to jump to the address;
creating the jump table associating the label of each basic block with the address of the basic block;
creating a coefficient array comprising diffusion function coefficients;
replacing a jump terminal instruction by a jump to the dispatcher basic block in each of the plurality of basic blocks;
creating and inserting into each of the plurality of basic blocks the first branch function and the second branch function, the first branch function taking as input a result of a diffusion function in the present basic block and outputting the label of the subsequent basic block to execute, the second branch function taking as input a diffusion function index of the present basic block and outputting a diffusion function index of the subsequent basic block based;
creating and inserting a transition function into the dispatcher basic block, the transition function obtaining the address in the jump table based on at least the label of the subsequent basic block; and
creating and inserting diffusion functions in the plurality of basic blocks, a diffusion function inserted in the present basic block calculating a value based on at least one diffusion function coefficient dependant on the diffusion function index and the label of the present basic block.

2. The method of claim 1, wherein the transition function obtains the address in the jump table based on also at least one diffusion function coefficient of the subsequent basic block.

3. The method of claim 1, further comprising, before creating and inserting the branch functions, assigning a value to the label of each of the plurality of basic blocks.

4. The method of claim 1, further comprising obtaining coefficients for each of the plurality of basic blocks.

5. A device for control flow graph flattening of a function of software code, the function comprising a plurality of basic blocks each having an address and at least one instruction, the device comprising a hardware processor configured to:
assign a value to at least one diffusion function coefficient of a diffusion function of each of the plurality of basic blocks;
calculate changes for labels and diffusion functions indices upon transitions between the plurality of basic blocks, the changes corresponding to modifications made by a linear first branch function and a linear second branch function;
create a dispatcher basic block comprising instructions to look up an address in a jump table and to jump to the address;
create the jump table associating a label of each basic block with the address of the basic block;
create a coefficient array comprising diffusion function coefficients;
replace a jump terminal instruction by a jump to the dispatcher basic block in each of the plurality of basic blocks;
create and insert into each of the plurality of basic blocks the first branch function and the second branch function, the first branch function taking as input a result of a diffusion function in the present basic block and outputting the label of the subsequent basic block to execute, the second branch function taking as input a diffusion function index of the present basic block and outputting a diffusion function index of the subsequent basic block based;
create and insert a transition function into the dispatcher basic block, the transition function obtaining the address in the jump table based on at least the label of the subsequent basic block; and
create and insert diffusion functions in the plurality of basic blocks, a diffusion function inserted in the present basic block calculating a value based on at least one diffusion function coefficient dependant on the diffusion function index and the label of the present basic block; and
an interface configured to output the function after control flow graph flattening.

6. The device of claim 5, wherein the processor is further configured to make the transition function obtain the address in the jump table based on also at least one diffusion function coefficient of the subsequent basic block.

7. The device of claim 5, wherein the processor is further configured to assign a value to the label of each of the plurality of basic blocks, before creating and inserting the branch functions.

8. The device of claim 5, wherein the processor is further configured to obtain coefficients for each of the plurality of basic blocks.