transform.c

   1 /*
   2  * Copyright (C) 2014-2015 Intel Corporation.
   3  */
   4
   5 #include <stdlib.h>
   6 #include <math.h>
   7 #include <utils/Log.h>
   8 #include <cutils/properties.h>
   9 #include <hardware/sensors.h>
  10 #include "calibration.h"
  11 #include "common.h"
  12 #include "description.h"
  13 #include "transform.h"
  14 #include "utils.h"
  15 #include "filtering.h"
  16 #include "enumeration.h"
  17
  18 #define GYRO_MIN_SAMPLES 5 /* Drop first few gyro samples after enable */
  19
  20
  21 /*----------------------------------------------------------------------------*/
  22
  23
  24 /* Macros related to Intel Sensor Hub */
  25
  26 #define GRAVITY         9.80665
  27
  28 /* 720 LSG = 1G */
  29 #define LSG             1024.0
  30 #define NUMOFACCDATA    8.0
  31
  32 /* Conversion of acceleration data to SI units (m/s^2) */
  33 #define CONVERT_A       (GRAVITY_EARTH / LSG / NUMOFACCDATA)
  34 #define CONVERT_A_X(x)  ((float(x) / 1000) * (GRAVITY * -1.0))
  35 #define CONVERT_A_Y(x)  ((float(x) / 1000) * (GRAVITY * 1.0))
  36 #define CONVERT_A_Z(x)  ((float(x) / 1000) * (GRAVITY * 1.0))
  37
  38 /* Conversion of magnetic data to uT units */
  39 #define CONVERT_M       (1.0 / 6.6)
  40 #define CONVERT_M_X     (-CONVERT_M)
  41 #define CONVERT_M_Y     (-CONVERT_M)
  42 #define CONVERT_M_Z     (CONVERT_M)
  43
  44 #define CONVERT_GAUSS_TO_MICROTESLA(x)  ((x) * 100)
  45
  46 /* Conversion of orientation data to degree units */
  47 #define CONVERT_O       (1.0 / 64)
  48 #define CONVERT_O_A     (CONVERT_O)
  49 #define CONVERT_O_P     (CONVERT_O)
  50 #define CONVERT_O_R     (-CONVERT_O)
  51
  52 /* Conversion of gyro data to SI units (radian/sec) */
  53 #define CONVERT_GYRO    (2000.0 / 32767 * M_PI / 180)
  54 #define CONVERT_GYRO_X  (-CONVERT_GYRO)
  55 #define CONVERT_GYRO_Y  (-CONVERT_GYRO)
  56 #define CONVERT_GYRO_Z  (CONVERT_GYRO)
  57
  58 #define BIT(x) (1 << (x))
  59
  60 inline unsigned int set_bit_range (int start, int end)
  61 {
  62         int i;
  63         unsigned int value = 0;
  64
  65         for (i = start; i < end; ++i)
  66                 value |= BIT(i);
  67
  68         return value;
  69 }
  70
  71 inline float convert_from_vtf_format (int size, int exponent, unsigned int value)
  72 {
  73         int divider = 1;
  74         int i;
  75         float sample;
  76         float mul = 1.0;
  77
  78         value = value & set_bit_range(0, size * 8);
  79
  80         if (value & BIT(size*8-1)) {
  81                 value =  ((1LL << (size * 8)) - value);
  82                 mul = -1.0;
  83         }
  84
  85         sample = value * 1.0;
  86
  87         if (exponent < 0) {
  88                 exponent = abs(exponent);
  89                 for (i = 0; i < exponent; ++i)
  90                         divider = divider * 10;
  91
  92                 return mul * sample/divider;
  93         }
  94
  95         return mul * sample * pow(10.0, exponent);
  96 }
  97
  98 /* Platform sensor orientation */
  99 #define DEF_ORIENT_ACCEL_X      -1
 100 #define DEF_ORIENT_ACCEL_Y      -1
 101 #define DEF_ORIENT_ACCEL_Z      -1
 102
 103 #define DEF_ORIENT_GYRO_X       1
 104 #define DEF_ORIENT_GYRO_Y       1
 105 #define DEF_ORIENT_GYRO_Z       1
 106
 107 /* G to m/s^2 */
 108 #define CONVERT_FROM_VTF16(s,d,x)      convert_from_vtf_format(s,d,x)
 109 #define CONVERT_A_G_VTF16E14_X(s,d,x)  (DEF_ORIENT_ACCEL_X * convert_from_vtf_format(s,d,x) * GRAVITY)
 110 #define CONVERT_A_G_VTF16E14_Y(s,d,x)  (DEF_ORIENT_ACCEL_Y * convert_from_vtf_format(s,d,x) * GRAVITY)
 111 #define CONVERT_A_G_VTF16E14_Z(s,d,x)  (DEF_ORIENT_ACCEL_Z * convert_from_vtf_format(s,d,x) * GRAVITY)
 112
 113 /* Degree/sec to radian/sec */
 114 #define CONVERT_G_D_VTF16E14_X(s,d,x)  (DEF_ORIENT_GYRO_X * convert_from_vtf_format(s,d,x) * M_PI / 180)
 115 #define CONVERT_G_D_VTF16E14_Y(s,d,x)  (DEF_ORIENT_GYRO_Y * convert_from_vtf_format(s,d,x) * M_PI / 180)
 116 #define CONVERT_G_D_VTF16E14_Z(s,d,x)  (DEF_ORIENT_GYRO_Z * convert_from_vtf_format(s,d,x) * M_PI / 180)
 117
 118 /* Milli gauss to micro tesla */
 119 #define CONVERT_M_MG_VTF16E14_X(s,d,x) (convert_from_vtf_format(s,d,x) / 10)
 120 #define CONVERT_M_MG_VTF16E14_Y(s,d,x) (convert_from_vtf_format(s,d,x) / 10)
 121 #define CONVERT_M_MG_VTF16E14_Z(s,d,x) (convert_from_vtf_format(s,d,x) / 10)
 122
 123
 124 static int64_t sample_as_int64 (unsigned char* sample, datum_info_t* type)
 125 {
 126         uint64_t u64;
 127         int i;
 128         int zeroed_bits = type->storagebits - type->realbits;
 129         uint64_t sign_mask;
 130         uint64_t value_mask;
 131
 132         u64 = 0;
 133
 134         if (type->endianness == 'b')
 135                 for (i=0; i<type->storagebits/8; i++)
 136                         u64 = (u64 << 8) | sample[i];
 137         else
 138                 for (i=type->storagebits/8 - 1; i>=0; i--)
 139                         u64 = (u64 << 8) | sample[i];
 140
 141         u64 = (u64 >> type->shift) & (~0ULL >> zeroed_bits);
 142
 143         if (type->sign == 'u')
 144                 return (int64_t) u64; /* We don't handle unsigned 64 bits int */
 145
 146         /* Signed integer */
 147
 148         switch (type->realbits) {
 149                 case 0 ... 1:
 150                         return 0;
 151
 152                 case 8:
 153                         return (int64_t) (int8_t) u64;
 154
 155                 case 16:
 156                         return (int64_t) (int16_t) u64;
 157
 158                 case 32:
 159                         return (int64_t) (int32_t) u64;
 160
 161                 case 64:
 162                         return (int64_t) u64;
 163
 164                 default:
 165                         sign_mask = 1 << (type->realbits-1);
 166                         value_mask = sign_mask - 1;
 167
 168                         if (u64 & sign_mask)
 169                                 return - ((~u64 & value_mask) + 1);     /* Negative value: return 2-complement */
 170                         else
 171                                 return (int64_t) u64;                   /* Positive value */
 172         }
 173 }
 174
 175
 176 static void reorder_fields (float* data, unsigned char map[MAX_CHANNELS])
 177 {
 178         int i;
 179         float temp[MAX_CHANNELS];
 180
 181         for (i=0; i<MAX_CHANNELS; i++)
 182                 temp[i] = data[map[i]];
 183
 184         for (i=0; i<MAX_CHANNELS; i++)
 185                 data[i] = temp[i];
 186 }
 187
 188
 189 static void clamp_gyro_readings_to_zero (int s, sensors_event_t* data)
 190 {
 191         float x, y, z;
 192         float near_zero;
 193
 194         x = data->data[0];
 195         y = data->data[1];
 196         z = data->data[2];
 197
 198         /* If we're calibrated, don't filter out as much */
 199         if (sensor[s].cal_level > 0)
 200                 near_zero = 0.02; /* rad/s */
 201         else
 202                 near_zero = 0.1;
 203
 204         /* If motion on all axes is small enough */
 205         if (fabs(x) < near_zero && fabs(y) < near_zero && fabs(z) < near_zero) {
 206
 207                 /*
 208                  * Report that we're not moving at all... but not exactly zero as composite sensors (orientation, rotation vector) don't
 209                  * seem to react very well to it.
 210                  */
 211
 212                 data->data[0] *= 0.000001;
 213                 data->data[1] *= 0.000001;
 214                 data->data[2] *= 0.000001;
 215         }
 216 }
 217
 218
 219 static void process_event_gyro_uncal (int s, int i, sensors_event_t* data)
 220 {
 221         gyro_cal_t* gyro_data;
 222
 223         if (sensor[s].type == SENSOR_TYPE_GYROSCOPE) {
 224                 gyro_data = (gyro_cal_t*) sensor[s].cal_data;
 225
 226                 memcpy(&sensor[i].sample, data, sizeof(sensors_event_t));
 227
 228                 sensor[i].sample.type = SENSOR_TYPE_GYROSCOPE_UNCALIBRATED;
 229                 sensor[i].sample.sensor = s;
 230
 231                 sensor[i].sample.data[0] = data->data[0] + gyro_data->bias_x;
 232                 sensor[i].sample.data[1] = data->data[1] + gyro_data->bias_y;
 233                 sensor[i].sample.data[2] = data->data[2] + gyro_data->bias_z;
 234
 235                 sensor[i].sample.uncalibrated_gyro.bias[0] = gyro_data->bias_x;
 236                 sensor[i].sample.uncalibrated_gyro.bias[1] = gyro_data->bias_y;
 237                 sensor[i].sample.uncalibrated_gyro.bias[2] = gyro_data->bias_z;
 238
 239                 sensor[i].report_pending = 1;
 240         }
 241 }
 242
 243 static void process_event_magn_uncal (int s, int i, sensors_event_t* data)
 244 {
 245         compass_cal_t* magn_data;
 246
 247         if (sensor[s].type == SENSOR_TYPE_MAGNETIC_FIELD) {
 248                 magn_data = (compass_cal_t*) sensor[s].cal_data;
 249
 250                 memcpy(&sensor[i].sample, data, sizeof(sensors_event_t));
 251
 252                 sensor[i].sample.type = SENSOR_TYPE_MAGNETIC_FIELD_UNCALIBRATED;
 253                 sensor[i].sample.sensor = s;
 254
 255                 sensor[i].sample.data[0] = data->data[0] + magn_data->offset[0][0];
 256                 sensor[i].sample.data[1] = data->data[1] + magn_data->offset[1][0];
 257                 sensor[i].sample.data[2] = data->data[2] + magn_data->offset[2][0];
 258
 259                 sensor[i].sample.uncalibrated_magnetic.bias[0] = magn_data->offset[0][0];
 260                 sensor[i].sample.uncalibrated_magnetic.bias[1] = magn_data->offset[1][0];
 261                 sensor[i].sample.uncalibrated_magnetic.bias[2] = magn_data->offset[2][0];
 262
 263                 sensor[i].report_pending = 1;
 264         }
 265 }
 266
 267 static void process_event (int s, sensors_event_t* data)
 268 {
 269         /*
 270          * This gets the real event (post process - calibration, filtering & co.) and makes it into a virtual one.
 271          * The specific processing function for each sensor will populate the necessary fields and set up the report pending flag.
 272          */
 273
 274          int i;
 275
 276          /* Go through out virtual sensors and check if we can use this event */
 277          for (i = 0; i < sensor_count; i++)
 278                 switch (sensor[i].type) {
 279                         case SENSOR_TYPE_GYROSCOPE_UNCALIBRATED:
 280                                 process_event_gyro_uncal(s, i, data);
 281                                 break;
 282                         case SENSOR_TYPE_MAGNETIC_FIELD_UNCALIBRATED:
 283                                 process_event_magn_uncal(s, i, data);
 284                                 break;
 285                         default:
 286                                 break;
 287                 }
 288 }
 289
 290
 291 static int finalize_sample_default (int s, sensors_event_t* data)
 292 {
 293         /* Swap fields if we have a custom channel ordering on this sensor */
 294         if (sensor[s].quirks & QUIRK_FIELD_ORDERING)
 295                 reorder_fields(data->data, sensor[s].order);
 296
 297         sensor[s].event_count++;
 298
 299         switch (sensor[s].type) {
 300                 case SENSOR_TYPE_ACCELEROMETER:
 301                         /* Always consider the accelerometer accurate */
 302                         data->acceleration.status = SENSOR_STATUS_ACCURACY_HIGH;
 303                         if (sensor[s].quirks & QUIRK_BIASED)
 304                                 calibrate_accel(s, data);
 305                         denoise(s, data);
 306                         break;
 307
 308                 case SENSOR_TYPE_MAGNETIC_FIELD:
 309                         calibrate_compass (s, data);
 310                         denoise(s, data);
 311                         break;
 312
 313                 case SENSOR_TYPE_GYROSCOPE:
 314
 315                         /* Report medium accuracy by default ; higher accuracy levels will be reported once, and if, we achieve  calibration. */
 316                         data->gyro.status = SENSOR_STATUS_ACCURACY_MEDIUM;
 317
 318                         /*
 319                          * We're only trying to calibrate data from continuously firing gyroscope drivers, as motion based ones use
 320                          * movement thresholds that may lead us to incorrectly estimate bias.
 321                          */
 322                         if (sensor[s].selected_trigger !=
 323                                 sensor[s].motion_trigger_name)
 324                                         calibrate_gyro(s, data);
 325
 326                         /*
 327                          * For noisy sensors drop a few samples to make sure we have at least GYRO_MIN_SAMPLES events in the
 328                          * filtering queue. This improves mean and std dev.
 329                          */
 330                         if (sensor[s].filter_type) {
 331                                 if (sensor[s].selected_trigger !=
 332                                     sensor[s].motion_trigger_name &&
 333                                     sensor[s].event_count < GYRO_MIN_SAMPLES)
 334                                                 return 0;
 335
 336                                 denoise(s, data);
 337                         }
 338
 339                         /* Clamp near zero moves to (0,0,0) if appropriate */
 340                         clamp_gyro_readings_to_zero(s, data);
 341                         break;
 342
 343                 case SENSOR_TYPE_LIGHT:
 344                 case SENSOR_TYPE_AMBIENT_TEMPERATURE:
 345                 case SENSOR_TYPE_TEMPERATURE:
 346                 case SENSOR_TYPE_INTERNAL_ILLUMINANCE:
 347                 case SENSOR_TYPE_INTERNAL_INTENSITY:
 348                         /* Only keep two decimals for these readings */
 349                         data->data[0] = 0.01 * ((int) (data->data[0] * 100));
 350
 351                         /* ... fall through ... */
 352
 353                 case SENSOR_TYPE_PROXIMITY:
 354                         /* These are on change sensors ; drop the sample if it has the same value as the previously reported one. */
 355                         if (data->data[0] == sensor[s].prev_val.data)
 356                                 return 0;
 357
 358                         sensor[s].prev_val.data = data->data[0];
 359                         break;
 360                 case SENSOR_TYPE_STEP_COUNTER:
 361                         if (data->u64.step_counter == sensor[s].prev_val.data64)
 362                                 return 0;
 363                         sensor[s].prev_val.data64 = data->u64.data[0];
 364                         break;
 365         }
 366
 367         /* If there are active virtual sensors depending on this one - process the event */
 368         if (sensor[s].ref_count)
 369                 process_event(s, data);
 370
 371
 372         return 1; /* Return sample to Android */
 373 }
 374
 375
 376 static float transform_sample_default (int s, int c, unsigned char* sample_data)
 377 {
 378         datum_info_t* sample_type = &sensor[s].channel[c].type_info;
 379         int64_t s64 = sample_as_int64(sample_data, sample_type);
 380         float scale = sensor[s].scale ? sensor[s].scale : sensor[s].channel[c].scale;
 381
 382         /* In case correction has been requested using properties, apply it */
 383         scale *= sensor[s].channel[c].opt_scale;
 384
 385         /* Apply default scaling rules */
 386         return (sensor[s].offset + s64) * scale;
 387 }
 388
 389
 390 static int finalize_sample_ISH (int s, sensors_event_t* data)
 391 {
 392         float pitch, roll, yaw;
 393
 394         /* Swap fields if we have a custom channel ordering on this sensor */
 395         if (sensor[s].quirks & QUIRK_FIELD_ORDERING)
 396                 reorder_fields(data->data, sensor[s].order);
 397
 398         if (sensor[s].type == SENSOR_TYPE_ORIENTATION) {
 399
 400                 pitch = data->data[0];
 401                 roll = data->data[1];
 402                 yaw = data->data[2];
 403
 404                 data->data[0] = 360.0 - yaw;
 405                 data->data[1] = -pitch;
 406                 data->data[2] = -roll;
 407         }
 408
 409         /* Add this event to our global records, for filtering purposes */
 410         record_sample(s, data);
 411
 412         return 1; /* Return sample to Android */
 413 }
 414
 415
 416 static float transform_sample_ISH (int s, int c, unsigned char* sample_data)
 417 {
 418         datum_info_t* sample_type = &sensor[s].channel[c].type_info;
 419         int val         = (int) sample_as_int64(sample_data, sample_type);
 420         float correction;
 421         int data_bytes  = (sample_type->realbits)/8;
 422         int exponent    = sensor[s].offset;
 423
 424         /* In case correction has been requested using properties, apply it */
 425         correction = sensor[s].channel[c].opt_scale;
 426
 427         switch (sensor_desc[s].type) {
 428                 case SENSOR_TYPE_ACCELEROMETER:
 429                         switch (c) {
 430                                 case 0:
 431                                         return correction * CONVERT_A_G_VTF16E14_X(data_bytes, exponent, val);
 432
 433                                 case 1:
 434                                         return correction * CONVERT_A_G_VTF16E14_Y(data_bytes, exponent, val);
 435
 436                                 case 2:
 437                                         return  correction * CONVERT_A_G_VTF16E14_Z(data_bytes, exponent, val);
 438                         }
 439                         break;
 440
 441                 case SENSOR_TYPE_GYROSCOPE:
 442                         switch (c) {
 443                                 case 0:
 444                                         return correction * CONVERT_G_D_VTF16E14_X(data_bytes, exponent, val);
 445
 446                                 case 1:
 447                                         return correction * CONVERT_G_D_VTF16E14_Y(data_bytes, exponent, val);
 448
 449                                 case 2:
 450                                         return  correction * CONVERT_G_D_VTF16E14_Z(data_bytes, exponent, val);
 451                         }
 452                         break;
 453
 454                 case SENSOR_TYPE_MAGNETIC_FIELD:
 455                         switch (c) {
 456                                 case 0:
 457                                         return correction * CONVERT_M_MG_VTF16E14_X(data_bytes, exponent, val);
 458
 459                                 case 1:
 460                                         return correction * CONVERT_M_MG_VTF16E14_Y(data_bytes, exponent, val);
 461
 462                                 case 2:
 463                                         return correction * CONVERT_M_MG_VTF16E14_Z(data_bytes, exponent, val);
 464                         }
 465                         break;
 466
 467                 case SENSOR_TYPE_LIGHT:
 468                         return (float) val;
 469
 470                 case SENSOR_TYPE_ORIENTATION:
 471                         return correction * convert_from_vtf_format(data_bytes, exponent, val);
 472
 473                 case SENSOR_TYPE_ROTATION_VECTOR:
 474                         return correction * convert_from_vtf_format(data_bytes, exponent, val);
 475         }
 476
 477         return 0;
 478 }
 479
 480
 481 void select_transform (int s)
 482 {
 483         char prop_name[PROP_NAME_MAX];
 484         char prop_val[PROP_VALUE_MAX];
 485         int i                   = sensor[s].catalog_index;
 486         const char *prefix      = sensor_catalog[i].tag;
 487
 488         sprintf(prop_name, PROP_BASE, prefix, "transform");
 489
 490         if (property_get(prop_name, prop_val, ""))
 491                 if (!strcmp(prop_val, "ISH")) {
 492                         ALOGI(  "Using Intel Sensor Hub semantics on %s\n", sensor[s].friendly_name);
 493
 494                         sensor[s].ops.transform = transform_sample_ISH;
 495                         sensor[s].ops.finalize = finalize_sample_ISH;
 496                         return;
 497                 }
 498
 499         sensor[s].ops.transform = transform_sample_default;
 500         sensor[s].ops.finalize = finalize_sample_default;
 501 }
 502
 503
 504 float acquire_immediate_float_value (int s, int c)
 505 {
 506         char sysfs_path[PATH_MAX];
 507         float val;
 508         int ret;
 509         int dev_num = sensor[s].dev_num;
 510         int i = sensor[s].catalog_index;
 511         const char* raw_path = sensor_catalog[i].channel[c].raw_path;
 512         const char* input_path = sensor_catalog[i].channel[c].input_path;
 513         float scale = sensor[s].scale ? sensor[s].scale : sensor[s].channel[c].scale;
 514         float offset = sensor[s].offset;
 515         int sensor_type = sensor_catalog[i].type;
 516         float correction;
 517
 518         /* In case correction has been requested using properties, apply it */
 519         correction = sensor[s].channel[c].opt_scale;
 520
 521         /* Acquire a sample value for sensor s / channel c through sysfs */
 522
 523         if (sensor[s].channel[c].input_path_present) {
 524                 sprintf(sysfs_path, BASE_PATH "%s", dev_num, input_path);
 525                 ret = sysfs_read_float(sysfs_path, &val);
 526
 527                 if (!ret)
 528                         return val * correction;
 529         }
 530
 531         if (!sensor[s].channel[c].raw_path_present)
 532                 return 0;
 533
 534         sprintf(sysfs_path, BASE_PATH "%s", dev_num, raw_path);
 535         ret = sysfs_read_float(sysfs_path, &val);
 536
 537         if (ret == -1)
 538                 return 0;
 539
 540         /*
 541          * There is no transform ops defined yet for raw sysfs values.
 542          * Use this function to perform transformation as well.
 543          */
 544         if (sensor_type == SENSOR_TYPE_MAGNETIC_FIELD)
 545                 return CONVERT_GAUSS_TO_MICROTESLA ((val + offset) * scale) * correction;
 546
 547         return (val + offset) * scale * correction;
 548 }
 549
 550 uint64_t acquire_immediate_uint64_value (int s, int c)
 551 {
 552         char sysfs_path[PATH_MAX];
 553         uint64_t val;
 554         int ret;
 555         int dev_num = sensor[s].dev_num;
 556         int i = sensor[s].catalog_index;
 557         const char* raw_path = sensor_catalog[i].channel[c].raw_path;
 558         const char* input_path = sensor_catalog[i].channel[c].input_path;
 559         float scale = sensor[s].scale ? sensor[s].scale : sensor[s].channel[c].scale;
 560         float offset = sensor[s].offset;
 561         int sensor_type = sensor_catalog[i].type;
 562         float correction;
 563
 564         /* In case correction has been requested using properties, apply it */
 565         correction = sensor[s].channel[c].opt_scale;
 566
 567         /* Acquire a sample value for sensor s / channel c through sysfs */
 568
 569         if (sensor[s].channel[c].input_path_present) {
 570                 sprintf(sysfs_path, BASE_PATH "%s", dev_num, input_path);
 571                 ret = sysfs_read_uint64(sysfs_path, &val);
 572
 573                 if (!ret)
 574                         return val * correction;
 575         };
 576
 577         if (!sensor[s].channel[c].raw_path_present)
 578                 return 0;
 579
 580         sprintf(sysfs_path, BASE_PATH "%s", dev_num, raw_path);
 581         ret = sysfs_read_uint64(sysfs_path, &val);
 582
 583         if (ret == -1)
 584                 return 0;
 585
 586         return (val + offset) * scale * correction;
 587 }