printf('------------ \n');
vector p =  point(1, 'P', 0);
addpoint(0, 0);
addpoint(0, p);
addprim(0, 'polyline', 0, 1);

vector ort_axis = 0;
float min_dist = 0;

for (int i = 0; i < 6; ++i)
{
    vector tmp = 0;
    tmp[i / 2] = (i % 2 * 2) - 1;
    
    float tmp_dist = distance(p, tmp);
    if (min_dist < tmp_dist)
    {
        min_dist = tmp_dist;
        ort_axis = floor(tmp * -1);
    }
}

printf('%d \n', ort_axis );

vector v = normalize ( chv("vec") );
vector dirs [ ] = array ( {1, 0, 0}, {0, 1, 0 }, {0, 0, 1} );

float maxdv = -1;
vector nearestdir = 0;

for ( int i = 0; i < len ( dirs ); ++i )
{
    float dv = dot ( v, dirs[ i ] );
    float adv = abs ( dv );
    if ( adv > maxdv )
    {
        maxdv = adv;
        nearestdir = dv > 0 ? dirs [ i ] : -dirs [ i ];
    }
}

int pt = addpoint ( 0, nearestdir );
setpointattrib ( 0, "N", pt, nearestdir );

vector v = chv("vector");
float longest = max( abs(v) );
vector axis = {0,0,0};
for ( int i=0; i<3; i++){
    if ( abs(v[i]) == longest ) {
        axis[i] = sign( v[i] );
        break;
    }
}