firmware/os/algos/calibration/diversity_checker/diversity_checker.c

/*
 * Copyright (C) 2016 The Android Open Source Project
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *      http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */

#include "calibration/diversity_checker/diversity_checker.h"

#include <stdarg.h>
#include <stdio.h>
#include <string.h>

#include "common/math/vec.h"
#include "chre/util/nanoapp/assert.h"

// Struct initialization.
void diversityCheckerInit(struct DiversityChecker* diverse_data,
                          const struct DiversityCheckerParameters* parameters) {
  CHRE_ASSERT_NOT_NULL(diverse_data);

  // Initialize parameters.
  diverse_data->threshold_tuning_param_sq =
      (parameters->threshold_tuning_param * parameters->threshold_tuning_param);
  diverse_data->max_distance_tuning_param_sq =
      (parameters->max_distance_tuning_param *
       parameters->max_distance_tuning_param);

  // Updating the threshold and max_distance using assumed local field.
  // Testing for zero and negative local_field.
  const float local_field =
      (parameters->local_field <= 0.0f) ? 1.0f : parameters->local_field;
  diversityCheckerLocalFieldUpdate(diverse_data, local_field);
  diverse_data->min_num_diverse_vectors = parameters->min_num_diverse_vectors;

  // Checking for min_num_diverse_vectors = 0.
  if (parameters->min_num_diverse_vectors < 1) {
    diverse_data->min_num_diverse_vectors = 1;
  }
  diverse_data->max_num_max_distance = parameters->max_num_max_distance;
  diverse_data->var_threshold = parameters->var_threshold;
  diverse_data->max_min_threshold = parameters->max_min_threshold;

  // Setting the rest to zero.
  diversityCheckerReset(diverse_data);

  // Debug Messages
#ifdef DIVERSE_DEBUG_ENABLE
  memset(&diverse_data->diversity_dbg, 0, sizeof(diverse_data->diversity_dbg));
#endif
}

// Reset
void diversityCheckerReset(struct DiversityChecker* diverse_data) {
  CHRE_ASSERT_NOT_NULL(diverse_data);
  // Clear data memory.
  memset(&diverse_data->diverse_data, 0, sizeof(diverse_data->diverse_data));

  // Resetting counters and data full bit.
  diverse_data->num_points = 0;
  diverse_data->num_max_dist_violations = 0;
  diverse_data->data_full = false;
}

bool diversityCheckerFindNearestPoint(struct DiversityChecker* diverse_data,
                                      float x, float y, float z) {
  // Converting three single inputs to a vector.
  const float vec[THREE_AXIS_DATA_DIM] = {x, y, z};

  // Result vector for vector difference.
  float vec_diff[THREE_AXIS_DATA_DIM];

  // normSquared result (k)
  float norm_squared_result;

  size_t i;

  // Running over all existing data points
  for (i = 0; i < diverse_data->num_points; ++i) {
    // v = v1 - v2;
    vecSub(vec_diff, &diverse_data->diverse_data[i * THREE_AXIS_DATA_DIM], vec,
           THREE_AXIS_DATA_DIM);

    // k = |v|^2
    norm_squared_result = vecNormSquared(vec_diff, THREE_AXIS_DATA_DIM);

    // if k < Threshold then leave the function.
    if (norm_squared_result < diverse_data->threshold) {
      return false;
    }

    // if k > max_distance, count and leave the function.
    if (norm_squared_result > diverse_data->max_distance) {
      diverse_data->num_max_dist_violations++;
      return false;
    }
  }
  return true;
}

void diversityCheckerUpdate(struct DiversityChecker* diverse_data, float x,
                            float y, float z) {
  CHRE_ASSERT_NOT_NULL(diverse_data);

  // If memory is full, no need to run through the data.
  if (!diverse_data->data_full) {
    // diversityCheckerDataSet() returns true, if input data is diverse against
    // the already stored.
    if (diversityCheckerFindNearestPoint(diverse_data, x, y, z)) {
      // Converting three single inputs to a vector.
      const float vec[THREE_AXIS_DATA_DIM] = {x, y, z};

      // Notice that the first data vector will be stored no matter what.
      memcpy(
          &diverse_data
               ->diverse_data[diverse_data->num_points * THREE_AXIS_DATA_DIM],
          vec, sizeof(float) * THREE_AXIS_DATA_DIM);

      // Count new data point.
      diverse_data->num_points++;

      // Setting data_full to true, if memory is full.
      if (diverse_data->num_points == NUM_DIVERSE_VECTORS) {
        diverse_data->data_full = true;
      }
    }
  }
}

bool diversityCheckerNormQuality(struct DiversityChecker* diverse_data,
                                 float x_bias, float y_bias, float z_bias) {
  CHRE_ASSERT_NOT_NULL(diverse_data);
  // If not enough diverse data points or max distance violations return false.
  if (diverse_data->num_points <= diverse_data->min_num_diverse_vectors ||
      diverse_data->num_max_dist_violations >=
          diverse_data->max_num_max_distance) {
    return false;
  }
  float vec_bias[THREE_AXIS_DATA_DIM] = {x_bias, y_bias, z_bias};
  float vec_bias_removed[THREE_AXIS_DATA_DIM];
  float norm_results;
  float acc_norm = 0.0f;
  float acc_norm_square = 0.0f;
  float max = 0.0f;
  float min = 0.0f;
  size_t i;
  for (i = 0; i < diverse_data->num_points; ++i) {
    // v = v1 - v_bias;
    vecSub(vec_bias_removed,
           &diverse_data->diverse_data[i * THREE_AXIS_DATA_DIM], vec_bias,
           THREE_AXIS_DATA_DIM);

    // norm = ||v||
    norm_results = vecNorm(vec_bias_removed, THREE_AXIS_DATA_DIM);

    // Accumulate for mean and VAR.
    acc_norm += norm_results;
    acc_norm_square += norm_results * norm_results;

    if (i == 0) {
      min = norm_results;
      max = norm_results;
    }
    // Finding min
    if (norm_results < min) {
      min = norm_results;
    }

    // Finding max.
    if (norm_results > max) {
      max = norm_results;
    }
    // can leave the function if max-min is violated
    // no need to continue.
    if ((max - min) > diverse_data->max_min_threshold) {
      return false;
    }
  }
  float inv = 1.0f / diverse_data->num_points;
  float var = (acc_norm_square - (acc_norm * acc_norm) * inv) * inv;

  // Debug Message.
#ifdef DIVERSE_DEBUG_ENABLE
  diverse_data->diversity_dbg.diversity_count++;
  diverse_data->diversity_dbg.var_log = var;
  diverse_data->diversity_dbg.mean_log = acc_norm * inv;
  diverse_data->diversity_dbg.max_log = max;
  diverse_data->diversity_dbg.min_log = min;
  memcpy(&diverse_data->diversity_dbg.diverse_data_log,
         &diverse_data->diverse_data,
         sizeof(diverse_data->diversity_dbg.diverse_data_log));
#endif
  return (var < diverse_data->var_threshold);
}

void diversityCheckerLocalFieldUpdate(struct DiversityChecker* diverse_data,
                                      float local_field) {
  if (local_field > 0) {
    // Updating threshold based on the local field information.
    diverse_data->threshold =
        diverse_data->threshold_tuning_param_sq * (local_field * local_field);

    // Updating max distance based on the local field information.
    diverse_data->max_distance = diverse_data->max_distance_tuning_param_sq *
                                 (local_field * local_field);
  }
}