vgvassilev
commented
2 years ago Hi @FROL256,
This works for me. Can you give a standalone example which crashes? I use Ubuntu clang version 12.0.1-++20210630032957+fed41342a82f-1~exp1~20210630133716.134
cat array.cpp
#include "clad/Differentiator/Differentiator.h"
void h(float a[3], float b, float output[]) {
output[0] = a[0] * a[0] * a[0];
output[1] = a[0] * a[0] * a[0] + b * b * b;
output[2] = 2 * (a[0] + b);
}
int main() {
auto dh = clad::gradient(h);
dh.dump();
return 0;
}
```bash
vvassilev@vv-nuc ~/workspace/builds/clad $ /usr/lib/llvm-12/bin/clang -x c++ -std=c++11 -Xclang -add-plugin -Xclang clad -Xclang -plugin-arg-clad -Xclang -fdump-derived-fn -Xclang -load -Xclang /home/vvassilev/workspace/builds/clad/lib/clad.so array.cpp -I/home/vvassilev/workspace/sources/clad/test/Hessian/../../include -lstdc++ -lm
void h_grad(float a[3], float b, float output[], clad::array_ref _d_a, clad::array_ref _d_b, clad::array_ref _d_output) {
float _t0;
float _t1;
float _t2;
float _t3;
float _t4;
float _t5;
float _t6;
float _t7;
float _t8;
float _t9;
float _t10;
float _t11;
float _t12;
_t2 = a[0];
_t1 = a[0];
_t3 = _t2 * _t1;
_t0 = a[0];
output[0] = _t3 * _t0;
_t6 = a[0];
_t5 = a[0];
_t7 = _t6 * _t5;
_t4 = a[0];
_t10 = b;
_t9 = b;
_t11 = _t10 * _t9;
_t8 = b;
output[1] = _t7 * _t4 + _t11 * _t8;
_t12 = (a[0] + b);
output[2] = 2 * _t12;
{
float _r_d2 = _d_output[2];
float _r12 = _r_d2 * _t12;
float _r13 = 2 * _r_d2;
_d_a[0] += _r13;
* _d_b += _r13;
_d_output[2] -= _r_d2;
_d_output[2];
}
{
float _r_d1 = _d_output[1];
float _r4 = _r_d1 * _t4;
float _r5 = _r4 * _t5;
_d_a[0] += _r5;
float _r6 = _t6 * _r4;
_d_a[0] += _r6;
float _r7 = _t7 * _r_d1;
_d_a[0] += _r7;
float _r8 = _r_d1 * _t8;
float _r9 = _r8 * _t9;
* _d_b += _r9;
float _r10 = _t10 * _r8;
* _d_b += _r10;
float _r11 = _t11 * _r_d1;
* _d_b += _r11;
_d_output[1] -= _r_d1;
_d_output[1];
}
{
float _r_d0 = _d_output[0];
float _r0 = _r_d0 * _t0;
float _r1 = _r0 * _t1;
_d_a[0] += _r1;
float _r2 = _t2 * _r0;
_d_a[0] += _r2;
float _r3 = _t3 * _r_d0;
_d_a[0] += _r3;
_d_output[0] -= _r_d0;
_d_output[0];
}
}
vvassilev@vv-nuc ~/workspace/builds/clad $ ./a.out
The code is:
void h_grad(float a[3], float b, float output[], clad::array_ref _d_a, clad::array_ref _d_b, clad::array_ref _d_output) {
float _t0;
float _t1;
float _t2;
float _t3;
float _t4;
float _t5;
float _t6;
float _t7;
float _t8;
float _t9;
float _t10;
float _t11;
float _t12;
_t2 = a[0];
_t1 = a[0];
_t3 = _t2 * _t1;
_t0 = a[0];
output[0] = _t3 * _t0;
_t6 = a[0];
_t5 = a[0];
_t7 = _t6 * _t5;
_t4 = a[0];
_t10 = b;
_t9 = b;
_t11 = _t10 * _t9;
_t8 = b;
output[1] = _t7 * _t4 + _t11 * _t8;
_t12 = (a[0] + b);
output[2] = 2 * _t12;
{
float _r_d2 = _d_output[2];
float _r12 = _r_d2 * _t12;
float _r13 = 2 * _r_d2;
_d_a[0] += _r13;
* _d_b += _r13;
_d_output[2] -= _r_d2;
_d_output[2];
}
{
float _r_d1 = _d_output[1];
float _r4 = _r_d1 * _t4;
float _r5 = _r4 * _t5;
_d_a[0] += _r5;
float _r6 = _t6 * _r4;
_d_a[0] += _r6;
float _r7 = _t7 * _r_d1;
_d_a[0] += _r7;
float _r8 = _r_d1 * _t8;
float _r9 = _r8 * _t9;
* _d_b += _r9;
float _r10 = _t10 * _r8;
* _d_b += _r10;
float _r11 = _t11 * _r_d1;
* _d_b += _r11;
_d_output[1] -= _r_d1;
_d_output[1];
}
{
float _r_d0 = _d_output[0];
float _r0 = _r_d0 * _t0;
float _r1 = _r0 * _t1;
_d_a[0] += _r1;
float _r2 = _t2 * _r0;
_d_a[0] += _r2;
float _r3 = _t3 * _r_d0;
_d_a[0] += _r3;
_d_output[0] -= _r_d0;
_d_output[0];
}
}
```
Hello! I took an example of jacobian from main page (h). If i make the first parameter array, clad crashes :(
Input:
Error:
PS: reproduces both with current version in master and release v0.9. Thank you very much! :)