Incorrect Jacobian generated when Order of assignment to Outputs is changed

[Nirhar]

Minimum reproducible example:

void h(float a, float output[]){
  output[1]=a;
  output[0]=a*a;  
}
int main(){
  auto dh = clad::jacobian(h);
  float a=3;
  float output[2]={0};
  float jc[2]={0};
  dh.execute(a,output,jc);
  printf("Result = {%.2f, %.2f}",jc[0],jc[1]);
}

The above example prints: Result = {0.00, 6.00}, which is wrong

However if the order of assignment to outputs is correctly given:

void h(float a, float output[]){
  output[0]=a*a;  
  output[1]=a;
}
int main(){
  auto dh = clad::jacobian(h);
  float a=3;
  float output[2]={0};
  float jc[2]={0};
  dh.execute(a,output,jc);
  printf("Result = {%.2f, %.2f}",jc[0],jc[1]);
}

The result is correct: Result = {6.00, 1.00}

[Nirhar]

The issue seems to be in this line, which expects sequential ordering of references to the output array