くれなゐの雑記

身の回りの人や自分が困っていたことを記事にしています

OpenFOAMの終了時間予測ツールを作ってみた

動機(ポエム) よく先生にこのような質問をされる 「その計算はいつ終わるんですか」 と 僕はこたえる 「今まで通りならだいたい6時間くらいで終わります」 と じゃあメッシュを細かくしたらどうか 計算条件が変わっても同じでも同じとは言えるのか そういう…

高専から大学編入をしてやったこと変わったこと

藤童子(@fwarashi)だったりkurenaifだったりと呼ばれています. 高専の方は5年生だったらギリギリ私のことを知っているかもしれません 神戸高専の出身で,神戸大学に編入してそのまま神戸大学の大学院にいます. ちょっと技術的な話のストックがマニアックな…

RECRUIT HOLDINGS SUMMER INTERNSHIP 2017 の ENGINEER コースに参加してきました

学業ガチ勢おじさんですが、気が変わったのでインターンに行ってきました。 結果として、とても楽しいインターンでよかったです。 動機 来年の就活が不安で仕方がないというツイートをポストしたら、知り合いから就活のしかたについて「インターンに行く。 …

OpenFOAMのcase fileをgitで管理する

OpenFOAMのcase fileをgitで管理する お久しぶりです。くれなゐです gitをソースコードの管理以外で使ったことがない方っていませんか? gitは非常に便利なツールで、様々な履歴を残すことができます。 今日はOpenFOAMのケースファイルを例に、ソースコード…

ipython notebookの環境整備

ipython notebookの準備 anacondaを入れる終わり。 ipythonの各種プラグイン設定 以下の記事を見ながらやったらpep8, autopep8, vimが入りました。qiita.comqiita.comqiita.com

SnackDown Elimination 2017 PLUSMUL

問題 https://www.codechef.com/SNCKEL17/problems/PLUSMUL 問題要約(相方担当) N個の数列が与えられる それぞれの数の間に"*" か "+"を入れた式を作る(通りの式ができる) 個の式の総和 MOD を出力せよ. 概要 なんとなくDPで解けそうだったので適当にやった…

いろんな言語でSI単位系の計算をやってみよう(Cpp, python, rust, Julia)(WIP)

C++ C++なのでコンパイル時になんやかんやできるやつを探してたらBoostであった Boost.Units を使えば良い.(http://www.boost.org/doc/libs/1_64_0/doc/html/boost_units.html) sample 使用例として,簡単な問題を解いてみる. h=634m地点から初速v0=2.4m/s…

Juliaで競技プログラミング(もうしない, ABC061)

Motivation 研究でちょっと触っているのでABCで使ってみた 解説は色々みたらわかるので割愛 A問題 http://abc061.contest.atcoder.jp/tasks/abc061_a A,B,C = parse.(split(readline(STDIN))) if A <= C && C <= B println("Yes") else println("No") end ht…

rustをインストールから行列で連立方程式を解くところまで

Motivation なぜrust? 最近C++がつらい rustが流行っているらしい https://imoz.jp/note/rust-functions.html 環境 ubuntu 16.04 64bit rustのインストール このコマンドで簡単にインストールする curl https://sh.rustup.rs -sSf | sh 以下のようなメッセー…

OpenFOAM-2.4.0をmacのgccでmakeする

Motivation Macはclangをつかってコンパイルをしている。 どうやらコンパイラによって微妙に挙動が違うらしくMacでgccが使いたいと先生に言われて色々やった。 OpenFOAM-2.4.0向け 概要 流れ的には gccを入れる openmpiを入れる flexを入れる OpenFOAMを入れ…

decomposeParでmanual切りをする

注意: バージョンによってファイル形式が異なる可能性があります.(OpenFOAM-2.4.0) Motivation 需要は少ないですが,manualでdecomposeParをする必要がありました. その手法について資料がないので解説します ケースファイル tutorials/incompressible/icoF…

OpenFOAMで計算結果以外をコピーするスクリプト

constantとおかsystemとかを残してコピーします 数字のファイルの判別方法が雑なのでミスってるかも… 一応これによって生じた不利益等は責任は負いません. #!/bin/bash if [ $# -ne 2 ]; then echo "usage: $0 source dist [first step]" exit 1 fi mkdir $…

pvpythonでアニメーション作成した

追記 2017-03-26 追記 --- pvbatchの並列実行 可視化対象 今回は適当にキャビティ流れを対象とします. tutorials/incompressible/icoFoam/cavity 実行 blockMesh icoFoam 可視化をするソースコード とりあえず実行 OpenFOAMの可視化はpvpythonを使って行い…

gdbがちょっと便利になるcgdb

OpenFOAMをデバッグしようとしてClionやらEclipseやら試してたんだけどどうもうまく動かない できればIDEでやりたいけどgdbを便利にしたものは無いものかと調べた結果がこれ インストール方法 ubuntuならapt-getで入った. macならbrewで入るらしい 使い方 …

epslatexなgnuplotで保存するたびにpdfを更新してそれを見ながら編集する(gnuplot, omake)

デモ ちょっとタイムラグと画質がきになりますがこんな感じです 今回はtexのフォントを使いたかったのでplatexを1回挟んでます makeplt_tex demo .plt(gnuplotのファイル) -> .eps+.tex -> .pdfファイル 以下のスクリプトを使います makeplt_texと名付けて…

OpenFOAMでgdbを使ってデバッグっぽいことをしてみる

OpenFOAMでgdbを使いたい Intro Motivation 先日 OpenCAE勉強会でgdbが便利だという話だという話をしました。 詳細について教えてほしいということなのでブログの記事を書くことにしました。 この記事では, 概要 デバッグオプションをいれたOpenFOAMのmake …

文書作成ソフト使い方まとめ

Kobe University Advent Calendar 2016 - Adventar の 16日目を担当させていただきました。さて, そろそろ卒論シーズンですねもうかきはじめている人もいるかと思いますが、皆様何を使っていますか? 文書作成ソフトですとWordが簡単で使いやすいみたいなのを…

ubuntuにsourcecodeProを入れる

タイトル通り こちらの記事をコピペするだけで行けました ubuntu 16.04askubuntu.com #!/bin/bash mkdir /tmp/adodefont cd /tmp/adodefont wget https://github.com/adobe-fonts/source-code-pro/archive/2.010R-ro/1.030R-it.zip unzip 1.030R-it.zip mkdi…

CODE FESTIVAL 2016 qual A. D - マス目と整数 / Grid and Integers

解説に乗ってるけどちょっとわかりにくかったので自分の言葉に置き換えて整理します。 あとイラレの練習( http://code-festival-2016-quala.contest.atcoder.jp/data/other/code-festival-2016-quala/editorial.pdf ) 問題 code-festival-2016-quala.contest…

CODE FESTIVAL 2016 qual A. C - 次のアルファベット / Next Letter

問題 code-festival-2016-quala.contest.atcoder.jp 解き方 この操作をして, 辞書順を小さくするためには, 'b'より大きいものを一周させて'a'にする以外ない.辞書順最小にするためには, 手前から見ていって, 'a'にできるものは'a'にして, それ以外は無視する…

Codeforces #712 Div.2 D. Memory and Scores

問題 codeforces.com 要約 が与えられる. : Memoryさんの初期点数 : 相手の初期点数 : 1ターンにの点数を得ることができる : 2tターン行うMemoryさんがかつことができるのは何通りか(MOD 1e9+7)以下のようなdpを組む その時間(time)で, Memoryさんのスコア-…

Codeforces #712 Div.2 C. Memory and De-Evolution

問題 codeforces.com 要約 正三角形2つが与えられる. それぞれの辺の長さは, [tex:x, y (y やること 非縮退三角形とは, 以下の条件を満たす三角形(つまり普通の三角形) 三角形の辺の長さを とすると, 少し変形して, 辺に着目すると, 以下のような式になる辺…

Codeforces #712 Div.2 B. Memory and Trident

問題 codeforces.com 要約 R,L,U,D (Right, Left, Up, Down)の4つからなる文字列が与えられる. それぞれの方向に1ずつ移動する Memoryさんは, 文字を書き換える能力がある. 文字を書き換えて, 最終的に原点に戻ってくるようにしたい. なん文字書き換えればよ…

Codeforces #712 Div.2 A. Memory and Crow

問題 codeforces.com 要約 数列 が与えられる. を満たすような を求めよ やること 後ろから考えると速い. が思いの外大きかったのでTLEに気をつけようを求めるんだけど, 総和の部分は以下のソースコードのように反転させながら足していかないとTLEする(2敗) …

c++でpython風のinput()を作った(戻り値型推論的なやつ)

Motivation int i = input(); string s = input(); みたいなのをしたい気分になった SourceCode 以下のやつをコピペすれば動きます struct input_returnner { template<typename T>operator T() const { T t; cin >> t; return t;} }; input_returnner input() { return </typename>…

AGC002 D. Stamp Rally

問題 agc002.contest.atcoder.jp やること 二分木探索で, 最大のスコアを探索する midまでの辺のfromとtoのノードを, union_findでuniteする union_findでは同じグループの数をuniteするときに一緒に計算しておく を含むグループの数と, を比較し, high, low…

東プレ REALFORCE 108UBK 変荷重キーボード/静電容量無接点/108キー/USB SJ08B0 レビュー

買ったやつ 東プレ REALFORCE 108UBK 変荷重キーボード/静電容量無接点/108キー/USB SJ08B0出版社/メーカー: 東プレ発売日: 2008/02/01メディア: Personal Computers購入: 11人 クリック: 1,063回この商品を含むブログ (22件) を見る 使った感想 非常に日本…

LOGICOOL ワイヤレストラックボール M570t レビュー

久々のレビュー記事 マウスをトラックボールにしてみました。他でもよく書かれていることとあまり書かれてないことを書きたいと思います。 買ったやつ LOGICOOL ワイヤレストラックボール M570t出版社/メーカー: ロジクール発売日: 2013/08/02メディア: Pers…

Codeforces #367 Div.2 C. Hard problem

problem codeforces.com 問題要約 コスト と 文字列 が与えられる。 をするために, コスト がかかるうまいことをして, 文字列を辞書順に昇順()となるようにしたい.昇順にできるとき, 最小の総コストを求めて, どう頑張ってもできない時は-1を出力する 解き方…

Codeforces #366 Div.2 C. Thor

Question Problem - C - Codeforces 問題要約 個のクエリが渡される クエリは以下の書式で渡される t x ここで, - の時: アプリケーションが通知を生成する - の時: アプリケーションが生成した通知をすべて読む (過去に読んだやつももう一度読む) - の時: …