Skip to content

Commit

Permalink
fix some file's encoding
Browse files Browse the repository at this point in the history
  • Loading branch information
LovelyA72 committed Nov 10, 2021
1 parent a5dba28 commit e64d017
Show file tree
Hide file tree
Showing 3 changed files with 221 additions and 221 deletions.
28 changes: 14 additions & 14 deletions SpaceWorld/src/matlabfunctions.cpp
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
// Matlabから移植した関数の寄せ集め
// Matlabから移植した関数の寄せ集め
#include <math.h>
#include "world.h"

Expand Down Expand Up @@ -42,7 +42,7 @@ void interp1(double *t, double *y, int iLen, double *t1, int oLen, double *y1)
s = (double *)malloc(sizeof(double) * oLen);
k = (int *)malloc(sizeof(int) * oLen);

// 初期設定
// 初期設定
for(i = 0;i < iLen-1;i++) h[i] = t[i+1]-t[i];
for(i = 0;i < oLen;i++) {p[i] = i; k[i] = 0;}

Expand All @@ -67,7 +67,7 @@ void interp1(double *t, double *y, int iLen, double *t1, int oLen, double *y1)
long decimateForF0(double *x, int xLen, double *y, int r)
{
// int r = 11;
int nfact = 9; // 多分これは固定でOK
int nfact = 9; // 多分これは固定でOK
double *tmp1, *tmp2;
tmp1 = (double *)malloc(sizeof(double) * (xLen + nfact*2));
tmp2 = (double *)malloc(sizeof(double) * (xLen + nfact*2));
Expand Down Expand Up @@ -99,7 +99,7 @@ void filterForDecimate(double *x, int xLen, double *y, int r)
{
double w[3], wt;
w[0] = w[1] = w[2] = 0.0;
double a[3], b[2]; // フィルタ係数 (r依存)
double a[3], b[2]; // フィルタ係数 (r依存)

switch(r)
{
Expand Down Expand Up @@ -163,7 +163,7 @@ void filterForDecimate(double *x, int xLen, double *y, int r)
}
}

// matlabに順ずる丸め
// matlabに順ずる丸め
double round(double x)
{
if(x > 0)
Expand All @@ -172,14 +172,14 @@ double round(double x)
return (int)(x-0.5);
}

// 差分
// 差分
void diff(const double *x, int x_length, double *y) {
for (int i = 0; i < x_length - 1; ++i) y[i] = x[i + 1] - x[i];
}


// サンプリング間隔が等間隔に限定し高速に動作するinterp1.
// 基本的には同じだが,配列の要素数を明示的に指定する必要がある.
// サンプリング間隔が等間隔に限定し高速に動作するinterp1.
// 基本的には同じだが,配列の要素数を明示的に指定する必要がある.
void interp1Q(double x, double shift, const double *y, int x_length,
const double *xi, int xi_length, double *yi) {
double *xi_fraction = new double[xi_length];
Expand All @@ -203,7 +203,7 @@ void interp1Q(double x, double shift, const double *y, int x_length,
delete[] delta_y;
}

// xorshift法と中心極限定理との組み合わせ
// xorshift法と中心極限定理との組み合わせ
float randn(void)
{
static unsigned int x = 123456789;
Expand All @@ -228,8 +228,8 @@ float randn(void)
return (float)tmp / 268435456.0f - 6.0f;
}

// fftfilt関数の移植
// yは,fftl分の長さを確保すること.
// fftfilt関数の移植
// yは,fftl分の長さを確保すること.
/*
void fftfilt(double *x, int xlen, double *h, int hlen, int fftl, double *y)
{
Expand Down Expand Up @@ -272,7 +272,7 @@ void fftfilt(double *x, int xlen, double *h, int hlen, int fftl, double *y)
fftw_destroy_plan(inverseFFT);
}
*/
// 2次元配列 (n*n)の逆行列を計算.メモリは確保しておくこと
// 2次元配列 (n*n)の逆行列を計算.メモリは確保しておくこと
void inv(double **r, int n, double **invr)
{
int i,j,k;
Expand All @@ -287,7 +287,7 @@ void inv(double **r, int n, double **invr)
}
for(i = 0;i < n;i++) invr[i][i] = 1.0;

// 配列の初期化
// 配列の初期化
//
for(i = 0;i < n;i++)
{
Expand All @@ -302,7 +302,7 @@ void inv(double **r, int n, double **invr)
}
}

// これで半分完了
// これで半分完了
for(i = n-1;i >= 0;i--)
{
for(j = 0;j < i;j++)
Expand Down
Loading

0 comments on commit e64d017

Please sign in to comment.