カテゴリー: テント村

ひさびさにJuliaを起こしてみたが，使い方をすっかり忘れてしまった！functionの書き方どころか，コメントの付け方さえ忘れている．二重階乗というプロジェクトがあったので，コピーしてここから始めることにしたのだが… 環境も変わってしまっているような感じだ．using Primes が通らないので，再インストールしている．Julia は間に GitHub が入っているのでいまいち何がどうなっているのかよく分からない．というか，その前に Juliaの拡張が入っていない．というか，２，３日前に Windows 11 をインストールしているのだが，それが裏目に出なければよいのだが…　エクスプローラの詳細表示で行間が広くなっている．こういう変化はうれしくない．

GeoGebraのIF文がどうやっても通らない．こんなつまらないところで躓くとは… 何かシンプルなサンプルはないだろうか？

▲q1とp1の交点T2(C_1)は辺CA上に乗っていない．わずかながらずれている．作図ミスと思われる．q1とp1はいずれも円だ．q1(l2, V), p1(T1, V)．Vは円周上の点でV(c1, j)だが，これが間違っているのではないか？jは(L, C)という直線でLは直線i上の点だ．i(A, B)でこれは間違いようがない．どうもLは点Aを選択したつもりで新たな点として追加されたものと思われる．⇒これで正しく動作するようになった．Point(…)を総点検する必要がある．Lは削除しておこう．

周長２等分線図の包絡線図を出せるようになった．できれば，この表示を切り替えるできるようにしたいのだが… オブジェクトの表示・非表示を切り替えるボタンを設置できるようになっている．チェックボックスにはオブジェクト・リストを選択できるようになっているので同時に切り替えが可能だ．ただし，オブジェクトの属性をここで切り替えることはできないので，直線を２種類作らなくてはならない．⇒それはそれでよいのではないか？しかし，このスイッチは表示・非表示を切り替えるだけなので，AとA’を切り替えることはできない…

リストを作るときにはオブジェクトを矩形で囲んで選択する．これは編集できるので，なんとかなりそうだ．ボタンを設置してスクリプトを書くという手も考えられる．オブジェクトの表示・非表示と属性（色と線の太さ）を変えたいのだが… 切断線は３本のはずなので，個別に書き換えてもそれほどの手間にはならないのではないか？

GeoGebraには軌跡の方程式を求める機能が備わっている！「数式処理 (CAS) エンジンでグレブナー基底を用いた計算を行い，軌跡の方程式を求めている」という．なんてこった！⇒ただし，求められるのは点の軌跡だけのようだ．つまり，直線の方程式を求めることはできない．まぁ，点の軌跡を得られればそれから割り出すことは可能だろう．

駆動点はなんらかの動線上の点でなくてはならない．P点をある円周上の点に設定し直すことはできるのではないか？その上で，小三角形上の点の軌跡の方程式を求めれば，その点とPを結ぶ直線の式は導き出せるはずだ．ただし，GeoGebraでは文字式は使えないので，方程式の係数はすべて数値になっているはずだ．それでも，方程式の形式がわかるだけでも有用性はあるだろう．⇒どうだろう？役に立つ場合もあるかもしれないが，あまり役には立たない気もする．

得られるのは軌跡の（駆動点とは独立の）方程式であり，駆動点との関係は捨象されてしまっている．このような式が与えられたとしても，それをどう利用すればよいのか？にわかには思いつかない．以下の記事ではGeoGebraを高校数学の授業に活用することを念頭に動作確認しているが，この機能にはまだまだかなりの限界があることが示されている．

https://www.kurims.kyoto-u.ac.jp/~kyodo/kokyuroku/contents/pdf/1951-15.pdf

久留島喜内の道具箱は結構作り込んであるつもりなのだが，まだまだだ．このツールは本来べき乗剰余の周期列を見るためのものだが，対象数の逆数をB進数表記で見ることもできるようになっている．べき乗と逆数には強い関連があるためだ．しかし，ときには対象数の逆というより，除数を分母としてみたい場合もある．そのような見地からすれば，現在の動作はγ＝１の場合を見ているに過ぎない．

また，現状ではB進数表記にもかなり制約があるように思われる．対象数が大きくなると簡単に表示不能になってしまう．もっと桁数の大きい数値を扱えるようにする必要がある．1/αの表示形式にもやや問題がある．間仕切りの｜は却って邪魔になっているような気がする※．仕様も使い方も忘れてしまったところがあり，stripeなどの意味もわからない．せめてヒントくらい出るようになっていないと… いや，ほとんど忘れてしまっていると言っても言い過ぎではない．マトリックスが２種出るようになっているが，なぜ２つも出しているのかも忘れてしまった．

※本家AMSツイッタの常連投稿者の pink @fourtyonerocks が4/13 = 0.(307692)という書き方をしている．これは読み易いような気がする．

前日のログで，「動作のおかしいところがある」としているが，やはりおかしいと思う．これらの点は⊿XYZの辺の延長と⊿ABCの辺の延長の交点だが，飛び出す点と飛び出さない点がある．飛び出さない点は上図にコーナーを回って隣の辺に移動している．この動作の方が正しいと思う．いや，それでは定義に反するものになるのではないか？たとえば，C1とC2はC点から等距離でなくてはならないはずだ．むしろそうなっていないA2の方が間違っている．A2はIntersect(L1, t1)となっているが，t1というのは⊿ABCそのものだ．

L1は線分ではなく，jとbの交点だ．いや，L1という点とL1という直線があるのではないか？L1というのはC点の別名のようだ．しかし，この修正はかなり難しい．１点を修正するとその副作用で複数の点が未定義になってしまう．A1, A2というラベルの付いている点は表示用で，オリジナルの点は直線上の点になっている．A1, A2のラベルの付いている点を非表示にして，オリジナルの点を表示するようにした．多分これで動くはずだ．⇒すべてのオブジェクトを表示設定して修正した．

三角形の周長２等分線図（包絡線）を GeoGebra に投稿した
https://www.geogebra.org/m/xa7nsrpg

動作のおかしいところがある．

Aの左にグリーンの◯が飛び出している．いや，これはノーマルなのかもしれない．辺ACは短いので周長/4の点が延長に出てしまうのだろう．まぁ，これは現在の仕様では仕方ないのかも知れない．２個の点を線分で結んでおけばよいのではないか？

周長２等分線図（全域版）をXに投稿した．これを見たのかどうかは分からないが，おしろい(Face Powder)とOmolesho Quadry Abiodunがフォローしてくれた．これでフォロワーが４人になった．

https://twitter.com/Cabbage1946/status/1787243696062341419

結構いい線まで行っているのだが，てこずっている．要領は大体分かったので，作り込んでゆけばよいはずだ．目標は，P点を全領域で動かせるようになることだ．つまり，P点がデカルト平面上のどこにあっても，与えられた三角形の周長等分線が引けるというのがゴールだ．P点が三角形の内部にあるときは，場合によっては等分線が２本になる場合もある．このような場合にも正しく動作するようにしなくてはならない．図面に書き込んでいる円の個数は１８個，多分これ以上増えることはないと思われるが，直線ないし線分の数はその何倍もあり，それをまた上回る数の点がある．おそらく，最終的にはオブジェクトの総数は数百個に達するのではないかという気がする…

この図面にはおそらく６個の象限があると推定されるが，現在のところそのうちの2.5象限しかカバーできていない．まだ未整備のところが多々あるためと考えられる．それらを潰してゆくしかないだろう．blueの領域がまったく出ていないのはなぜだろう．すでにF点は決めてあるのだが，N点が定まっていない．まず，ここから片付けることにしよう．F点というのは頂角の２等分線上の点でこの点とPの中点を中心とする円がN点の位置を決定することになっている．N点はL2上の点であり，L2というのは，３本の基線のうちの一つだ．

大体全域で動くようになった．ただし，まだ問題点がいくつかある．内部の三角形と小円が交差するとき，２点で爻わる場合，どちらを取るかという問題だ．どちらかを消さなくてはならないが，それを図形から判断できるだろうか？スクリプトを使って切り替えることができれば話は簡単なのだが… 実際，現在の実装では目視で明らかに２等分になっていない場合がある．⇒１本線を消したら正しく動作するようになった．⇒等分線が３本出てくる場合がある．この動作は正しい．

ナムの三角形周長２等分線作図問題の一般解をGeoGebraで書こうとしているのだが，なかなか難しい．できるだけ画面をすっきりさせたいのだが，余分なゴミが出てきてしまう．これはポイントが落ちてしまっているためだろう．たとえば，下図には３本千切れた尻尾出ている．

これの固着できない点はどこかの線分の延長上にあったのではないだろうか？⇒単純に消しておくだけでよいのではないか？どうも，そもそも仕様通りの動作になっていない．実装が悪いのか，設計が間違っているのか？もう一度作り直しするしかない．

真ん中の小さい三角形は包絡線図でかなり重要な図形になってくるが，⊿ABCの形状によっては描画できないことがある．線分との交点ではなく，延長上の交点を取ってみる．

また，CalcTrajectoryでエラーが出ている．⇒処理本体を関数化したため，変数が見えなくなっている．変数をトップスコープに移した．

▲PNGファイルの保存パスが存在しないというエラーになる．

png(“.\\PNG\\” * titleName())

現在の作業フォルダが蜜天GUIになっているためではないか？その下のmittenを開くようにした方がよい．⇒このフォルダは事前に作っておく必要がある．あるいは，むしろデスクトップに生成するようにした方が分かり易いかもしれないが…⇒これは後から考えることにして，とりあえずフォルダを作っておこう．

一応一通り動くようになった．⇒InfinitePoolGameをモジュール化したらまた動かなくなってしまった．isinteractiveをisinteractive()に変えたためだ．isinteractiveは変数ではないので，これまでそれが通っていたことの方がおかしいのだが，デバッガで走らせることはinteractiveには入らないのだろう．REPLからinclude(“src\\InfinitePoolGame.jl”)では想定通りの動作，つまり，isinteractiveはtrueになる．定数を切って切り分けるしかなさそうだ．技術的にかなり難しい問題が出てきた．

ようやく元のように動くようになった．特性関数を動的に切り替えるというのが目下の課題だ．julia には eval という関数があるのでなんとかなるのではないかと思っているのだが… まず，関数ポインタのようなものがあるかどうかを探さなくてはならない．⇒プロットを表示するのは以外と簡単にできた．

include(“InfinitePoolGame.jl”)

でInfinitePoolGame内の関数にアクセスできる．とりあえずは，これで十分だ．本当はimporないしusingでロードしたいのだが，まだパッケージになっていないためロードできない．includeは実行時に何回でも実行できるので，最悪の場合（他に打つ手がない場合）は，InfinitePoolGame.jlを直接書き換えて再ロードするということが考えられる．多分，これは不可能ではないと思われる．juliaは正規表現をサポートしているので，ソースファイルの書き換えはできるはずだ．一度ここでバックアップを取っておこう．

スタンドアローン型のJuliaアプリを作ろうとしていたのだが，断念した．ともかく，先に進むことにしよう．画面上にはまだボタン１つしか出ていないが，少なくとも，①関数入力用のテキストボックス，②グラフの範囲を設定するための数値ボックスｘ4，③実行ボタンが必要だ．ただし，それに入る前にこれまで開発したPlotsを使ったグラフ表示を組み込んで動作確認しておいた方がよい．

なぜだろう？Plotsが読み込めなくなっている．Juliah本体でも読めない．⇒pkgで読み直した．Rootsも入っていない．⇒追加した．最終版と思われるものを動かしているのだが，エラーになった．グラフは一応出ているが，軌跡は描画されていない．CalcTrajectoryの中でエラーが起きている．なぜだろう？メニューに中国語が出ている．

いままではGotoDefinitionが先頭にあったような気がするのだが… おかしい．元に戻っている．Julia には言語設定というのがあったろうか？CalcProjectoryの中でplot!を実行して例外が起きている．

ErrorException(“Cannot convert Float64 to series data for plotting”)

これまではこれで動いていたはずなのだが，Plotsをリロードしたことでバージョンが変わってしまったのだろうか？少しバージョンを戻したら動いた．多分，この辺りは仕掛りのままになっていたのだろう．蜜蜂天国 2024-04-05の前は2024-04-01で４日も空いている．ログを見てみよう．4月5日には「find_zeroが動作していない」というのがあり，修正して動作するようになったが，それ以上の記載はない．その一つ前は2024/03/29で「た軌跡の描線が出なくなっている」から始まる．

2024/04/05という日付のログは２本あり，２本目のログでは「蜜蜂天国のまとめに入っているところだが，足踏み状態だ．いろいろなパターンを試してみたいのだが，そのためにはもうちょっと準備が必要だ．」として，「まとめ」を打ち切り，「JuliaでGUIアプリ」という方向に走り出している．2024-04-01版ではCalcTrajectoryではplot!を実行する代わりにpush(plt,…)している．確かにこの違いは大きい．⇒論理を戻して動作するようになった．ただし，ウィンドウが出てこない．PNGとGIFは保存されている．どこが違うのだろう？⇒出た！show=:trueのコロンが落ちていた．うまく行ったようだ．アニメを作らないオプションでも動作するようになった．