C/C++/C#

MPIR（）をインストールしたので使ってみた。MPIRを使うと巨大な数を扱うことができるようになる。ソースコード例#include<mpir.h>#include<iostream>int main() { mpf_set_default_...

行列演算ライブラリは巨大なのでコンパイルに時間がかかる。

coreのほかにLUも必要。#include "../Eigen/core"#include "../Eigen/LU"あとはA.inverse()とすればいい。実行結果。-0.666667 -1.33333 1-0.666667 3.66...

行列の固有値の積は行列式の値と等しくなる。これをEigenで試す。実行結果固有値 (16.7075,0)(-0.905741,0) (0.198246,0)行列式-3確かめれば\begin{align}16.7075 \times (-0....

Eigenで固有値と固有ベクトルを求めるにはEigen::EigenSolver< Eigen::MatrixXf > s(A);とした後、sに対してs.eigenvalues() とすれば固有値を、s.eigenvectors()とすれば...

\(A\)と\(A\)の転地の積のトレースはフロベニウスノルムの二乗と等しくなる。つまり\begin{align}\mathrm{tr} (A A^{T}) = \mathrm{tr} (A^{T} A) = \sum_{i=1}^{n}\...

トレースにはつぎのようなの性質がある。\begin{align}\mathrm{tr} {A1} + \mathrm{tr}{A2} &= \mathrm{tr} ( A1+A2 ) \\a\ \mathrm{tr} {A1} &= \ma...

何度もstd~と書くのはめんどくさいので汎用print関数を作る。詳細は記事がたくさんあるので割愛。vectorなど工夫しないと使えないものもあるが、数値、文字関係なく使える。template <class T> void mytprint...

CやC++では配列を動的に確保することができないのでmallocやnewを使う。メモリ開放をしないと大変なことになるので注意。実行結果要素数は？->3数を入力->2数を入力->1数を入力->42,1,4,ソースコード。#include <i...

トレースには次の性質がある。\begin{align}tr (A+B) = tr A + tr B\end{align} 今回はこれをC++で計算する。実行結果1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,1...

正方行列の対角成分の和\begin{align}tr A = \sum_{i=1}^{n} a_{ii}\end{align}をトレースという。今回はこれをC++で計算する。実行結果1,2,3,4,5,6,7,8,9,15ソースコード#in...

フロベニウスノルムは\begin{align}\| A \|_{F}=\sqrt{\sum_{i=1}^{m} \sum_{j=1}^{n} \left |a_{ij} \right |}\end{align}で計算できる。一次元配列に行列...

ベクトルを定義するときは#include "../Eigen/Dense"をインクルードしてEigen::Vector3d x(1, 2, 3);とすればベクトルが定義できるので、内積と外積をx.dot(y)x.cross(y)で計算する。...

Eigenで行列の和・差・積を試す。「+」、「-」、「*」が使えるので直観的。ソースコード#include <iostream>#include "../Eigen/core"int main() { int n = 2; Eigen::M...

Eigen(は行列用の科学技術ライブラリで、ヘッダーだけ読み込めば動く。ダウンロードしたファイルを展開するとEigenというファイルがあるのでそれを適当な場所に置き#include "../Eigen/core"としてcoreにパスを通せば...

typedefを使うと何度も宣言しなくてよく、型のように扱える。typedef void (*ECHO)(int);ソースコード#include <iostream>typedef void (*ECHO)(int);void echo(i...

内容はC言語。関数ポインタはvoid (*pecho)(int);で定義してpecho = echo;で参照設定する。関数ポインタの定義は参照元の関数の返り値と引数の型void echo(int x) { std::cout << x;}が...

static_castを使えばできる y = static_cast<double>(x);型変換した後する前と変換後で型を調べる。実行結果intdoubleソースコード#include <iostream>#include <typein...

findを使うと先頭から何番目にその文字があるかを探すことができる実行結果184467440737095516150存在しませんちなみにこれらのコードはどれも同じ動作をする。autoの動作が原因。 if (str1.find("z") ==...

引数に応じて動作を変えられる。コンストラクタの定義はこの部分 Student() { } Student(std::string name, int jpn, int math, int eng) { m_name = name; m_ma...