2016-08-21 01:36:02 +08:00
|
|
|
/****************************************************************************
|
|
|
|
* libc/misc/lib_xorshift.c
|
|
|
|
*
|
|
|
|
* Copyright (C) 2016 Gregory Nutt. All rights reserved.
|
|
|
|
* Author: David S. Alessio
|
|
|
|
*
|
|
|
|
* Redistribution and use in source and binary forms, with or without
|
|
|
|
* modification, are permitted provided that the following conditions
|
|
|
|
* are met:
|
|
|
|
*
|
|
|
|
* 1. Redistributions of source code must retain the above copyright
|
|
|
|
* notice, this list of conditions and the following disclaimer.
|
|
|
|
* 2. Redistributions in binary form must reproduce the above copyright
|
|
|
|
* notice, this list of conditions and the following disclaimer in
|
|
|
|
* the documentation and/or other materials provided with the
|
|
|
|
* distribution.
|
|
|
|
* 3. Neither the name NuttX nor the names of its contributors may be
|
|
|
|
* used to endorse or promote products derived from this software
|
|
|
|
* without specific prior written permission.
|
|
|
|
*
|
|
|
|
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
|
|
|
|
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
|
|
|
|
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
|
|
|
|
* FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
|
|
|
|
* COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
|
|
|
|
* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
|
|
|
|
* BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS
|
|
|
|
* OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED
|
|
|
|
* AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
|
|
|
|
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
|
|
|
|
* ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
|
|
|
|
* POSSIBILITY OF SUCH DAMAGE.
|
|
|
|
*
|
|
|
|
****************************************************************************/
|
|
|
|
|
|
|
|
/* This random number generator is simple, fast and portable.
|
|
|
|
* Ref: https://en.wikipedia.org/wiki/Xorshift
|
|
|
|
*/
|
|
|
|
|
|
|
|
/****************************************************************************
|
|
|
|
* Included Files
|
|
|
|
****************************************************************************/
|
|
|
|
|
|
|
|
#include <nuttx/config.h>
|
2016-08-21 07:27:41 +08:00
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
#include <sys/types.h>
|
2016-08-21 01:36:02 +08:00
|
|
|
#include <stdint.h>
|
2016-08-21 21:47:53 +08:00
|
|
|
#include <assert.h>
|
2016-08-21 07:27:41 +08:00
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
#include <nuttx/lib/xorshift128.h>
|
2016-08-21 07:27:41 +08:00
|
|
|
|
2016-08-21 01:36:02 +08:00
|
|
|
/****************************************************************************
|
|
|
|
* Public Functions
|
|
|
|
****************************************************************************/
|
|
|
|
|
|
|
|
/****************************************************************************
|
2016-08-21 21:47:53 +08:00
|
|
|
* Name: xorshift128
|
2016-08-21 01:36:02 +08:00
|
|
|
*
|
|
|
|
* Description:
|
2016-08-21 21:47:53 +08:00
|
|
|
* Generate one 32-bit pseudo-random number.
|
2016-08-21 01:36:02 +08:00
|
|
|
*
|
2016-08-21 21:47:53 +08:00
|
|
|
* NOTE: Because the PRNG state is passed as a parameter, this function is
|
|
|
|
* fully re-entrant and may be called from an interrupt handler.
|
2016-08-21 01:36:02 +08:00
|
|
|
*
|
2016-08-21 21:47:53 +08:00
|
|
|
* The downside to this is that users of the PRNG might not get as much
|
|
|
|
* entropy as if it were a common state structure.
|
2016-08-21 01:36:02 +08:00
|
|
|
*
|
|
|
|
* Input Parameters:
|
2016-08-21 21:47:53 +08:00
|
|
|
* state - The current XorShift128 state.
|
2016-08-21 01:36:02 +08:00
|
|
|
*
|
|
|
|
* Returned Value:
|
|
|
|
* The generated pseudo-random number
|
|
|
|
*
|
|
|
|
****************************************************************************/
|
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
uint32_t xorshift128(FAR struct xorshift128_state_s *state)
|
2016-08-21 01:36:02 +08:00
|
|
|
{
|
|
|
|
uint32_t t;
|
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
DEBUGASSERT(state != NULL);
|
2016-08-21 07:27:41 +08:00
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
t = state->x;
|
2016-08-21 01:36:02 +08:00
|
|
|
t ^= t << 11;
|
|
|
|
t ^= t >> 8;
|
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
state->x = state->y;
|
|
|
|
state->y = state->z;
|
|
|
|
state->z = state->w;
|
2016-08-21 01:36:02 +08:00
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
state->w ^= state->w >> 19;
|
|
|
|
state->w ^= t;
|
2016-08-21 07:27:41 +08:00
|
|
|
|
2016-08-21 21:47:53 +08:00
|
|
|
return state->w;
|
2016-08-21 01:36:02 +08:00
|
|
|
}
|