Skip to content

Instantly share code, notes, and snippets.

Avatar

mk-mode komasaru

View GitHub Profile
@komasaru
komasaru / calc.cpp
Created Jun 20, 2022
C++ source code to calculate a simple regression curve(5d).
View calc.cpp
#include "calc.hpp"
/*
* @brief 単回帰曲線(5次)の計算
*
* @param[ref] a (double)
* @param[ref] b (double)
* @param[ref] c (double)
* @param[ref] d (double)
* @param[ref] e (double)
@komasaru
komasaru / regression_curve_5d.rb
Last active Jun 20, 2022
Ruby script to calculate a simple regression curve.(5d)
View regression_curve_5d.rb
#! /usr/local/bin/ruby
#*********************************************
# Ruby script to calculate a simple regression curve.
# : y = a + b * x + c * x^2 + d * x^3 + e * x^4 + f * x^5
# : 連立方程式を ガウスの消去法(ピボット選択)で解く方法
#*********************************************
#
class Array
def reg_curve_5d(y)
# 以下の場合は例外スロー
@komasaru
komasaru / gauss_elimination_pivot.cpp
Last active Apr 15, 2022
C++ source code to solve simultaneous equations with Gauss elimination method(pivot).
View gauss_elimination_pivot.cpp
/***********************************************************
連立方程式の解法 ( ガウスの消去法(ピボット選択) )
DATE AUTHOR VERSION
2021.09.16 mk-mode.com 1.00 新規作成
Copyright(C) 2021 mk-mode.com All Rights Reserved.
***********************************************************/
#include <cmath> // for fabs
#include <cstdlib> // for EXIT_XXXX
@komasaru
komasaru / gauss_elimination_pivot.f95
Created Apr 14, 2022
Fortran95 source code to solve simultaneous equations with Gauss elimination method(pivot).
View gauss_elimination_pivot.f95
!************************************************************
! Simultaneous equations solving by Gauss-Elimination(Pivot) method
!
! DATE AUTHOR VERSION
! 2022.04.14 mk-mode.com 1.00 新規作成
!
! Copyright(C) 2022 mk-mode.com All Rights Reserved.
!************************************************************
!
module const
@komasaru
komasaru / gauss_elimination_pivot.rb
Created Apr 13, 2022
Ruby script to solve simultaneous equations with Gauss elimination method(pivot).
View gauss_elimination_pivot.rb
#! /usr/local/bin/ruby
#*********************************************
# 連立方程式の解法 ( ガウスの消去法(ピボット選択) )
#*********************************************
#
class GaussElimination
def initialize
# 係数
@a = [
[1.0, 2.0, 7.0, 6.0, 6.0],
@komasaru
komasaru / date_loop.sh
Created Jul 2, 2021
Bash script to loop by date.
View date_loop.sh
#!/bin/bash
DATE_S="20201226"
DATE_E="20210105"
# Example-1
for (( dt=$DATE_S ; $dt < $DATE_E ; dt=`date -d "$dt 1 day" '+%Y%m%d'` )) ; do
echo $dt
done
echo "---"
@komasaru
komasaru / make_eop.cpp
Last active Jun 8, 2021
C++ source code to make an EOP(Earth Orientation Parameters) file.
View make_eop.cpp
/***********************************************************
IERS の Buttelin A テキストデータから
EOP(Polar Motion etc.) テキストファイルを生成
* 予め [こちら](ftp://ftp.iers.org/products/eop/rapid/) か
らダウンロードしておいたものを使用する。
(IAU 2000A 章動理論によるデータ
finals2000A.all", "finals2000A.daily")
* 1日のデータに速報値(区分"I")と確定値がある場合は、確定値
を優先。
* 2ファイルで重複する日付のデータは "finals2000A.daily" を
@komasaru
komasaru / blh2enu.cpp
Last active Sep 16, 2021
C++ source code to convert a coordinate from WGS84(BLH) to ENU.
View blh2enu.cpp
/***********************************************************
BLH -> ENU 変換
: WGS84 の緯度(Beta)/経度(Lambda)/楕円体高(Height)を
ENU (East/North/Up; 地平) 座標に変換する。
* 途中、 ECEF(Earth Centered Earth Fixed; 地球中心・地
球固定直交座標系)座標への変換を経由。
DATE AUTHOR VERSION
2021.05.06 mk-mode.com 1.00 新規作成
@komasaru
komasaru / ecef2blh.cpp
Last active May 6, 2021
C++ source code to convert a coordinate from ECEF to WGS84(BLH).
View ecef2blh.cpp
/***********************************************************
ECEF -> BLH 変換
: ECEF(Earth Centered Earth Fixed; 地球中心・地球固定直交座標系)座標を
WGS84 の緯度(Latitude)/経度(Longitude)/楕円体高(Height)に変換する。
DATE AUTHOR VERSION
2021.05.02 mk-mode.com 1.00 新規作成
Copyright(C) 2021 mk-mode.com All Rights Reserved.
@komasaru
komasaru / blh2ecef.cpp
Last active May 6, 2021
C++ source code to convert a coordinate from WGS84(BLH) to ECEF.
View blh2ecef.cpp
/***********************************************************
BLH -> ECEF 変換
: WGS84 の緯度(Beta)/経度(Lambda)/楕円体高(Height)を
ECEF(Earth Centered Earth Fixed; 地球中心・地球固定直交座標系)座標に
変換する。
DATE AUTHOR VERSION
2021.04.30 mk-mode.com 1.00 新規作成
Copyright(C) 2021 mk-mode.com All Rights Reserved.