トップページ⇒math
1001コメント517KB

〓 Mathematica 伍 〓

レス数が950を超えています。1000を超えると書き込みができなくなります。
0001１３２人目の素数さん2011/11/11(金) 09:02:28.82
　　　　　　　　_....._{{ 〃
　　　　 　, - ' ,..､､.ヾ{{フ'⌒`ヽ､
　　　　／　 ,:', -‐‐` ´ '´⌒ヽ　ヾ:､
.　　　,'　　 ,'´　,ィ ,ｨ ,' , 　 ｀ヽ',　 ',-<
　 　 ,' 　　.i　 /|. /.| { i,　 i,　　}.　 }_,,）)
　　　!　|　 ! .,'-.{ !　!|; |`､.}ﾞ!.!　|.　 !　ヽ.
　　　',　',　|Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l !　| 　 ﾞ',
　　　 ヽ､', l:!Kノ}. 　 　 f:_.)iﾞi: ﾘ !　l　ル　
　　　　 | l!iヽﾞ- '　, 　　.!__:ﾉ ﾞ ,ﾘ l ﾘ'´ 　 　　
.　　　　 ',|!!､ 　 　r‐┐ 　 ｀ ノ'. /,ｲ　　　　／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣ 　
　　　　　 'i!ﾞ､ヽ､　ﾞｰ'　　_, ィ,:',:''´　　　 ＜　 Mathematica に関する話題はここに書いてね！
　　　　　　ﾞ:､ｨ､jヾｰ::: 'iﾍ　.ノ',ﾘ. 　 　 　　 ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
　　　,､- '´　ヽ､ﾞ､　　 {　｀＞"､ 　
　　/＼＼　　　　',　　 }　　 //｀ヽ

過去スレ
〓Mathematica〓
http://science2.2ch.net/test/read.cgi/math/985023298/
〓 Mathematica 2 〓
http://science3.2ch.net/test/read.cgi/math/1078534285/
〓 Mathematica 3 〓
http://science6.2ch.net/test/read.cgi/math/1121413040/
〓 Mathematica 四 〓
http://kamome.2ch.net/test/read.cgi/math/1197270001/ (鯖飛びでログ消滅)
http://mimizun.com/log/2ch/math/1197270001/
〓 Mathematica 5 〓
http://kamome.2ch.net/test/read.cgi/math/1285859504/
0894１３２人目の素数さん2014/08/30(土) 14:11:32.46
>>893
RegionMeasure[ImplicitRegion[(c1^2+c2^2)^(-1/2)>=(c1^2/a1^2+c2^2/a2^2)^(-1/2)+1,{{a1,0,1},{a2,0,1},{c1,-1,1},{c2,-1,1}}]]
Integrate[1,Element[{a1,a2,c1,c2},ImplicitRegion[(c1^2+c2^2)^(-1/2)>=(c1^2/a1^2+c2^2/a2^2)^(-1/2)+1,{{a1,0,1},{a2,0,1},{c1,-1,1},{c2,-1,1}}]]]
orz
0895１３２人目の素数さん2014/09/06(土) 07:49:39.29
215（x^4+y^2+54）
とか
689（x^3+y^7+xy+xy^4+125）
の式から、係数「215」、「689」を抽出するには、どうしたら良いでしょうか？
0896１３２人目の素数さん2014/09/06(土) 18:16:46.06
Coefficient で係数を得てから、非零の係数から最大公約数でも計算したら
0897１３２人目の素数さん2014/09/07(日) 20:52:42.26
回答有難う御座います。
やってみます。
0898１３２人目の素数さん2014/09/07(日) 21:34:21.14
立体を切断するにはどうすればいいんだ?

knife.mは最近のバージョンでは使えないんだろ
PlotRangeだと斜めとか切れないし
0899１３２人目の素数さん2014/09/08(月) 00:34:23.87
knife.m現在はファイルすら見つからないね
著者の人は相変わらずMathematica使ってるみたいだけど
興味が折り紙に向かってる模様
http://www.jssst.or.jp/award/detail/id=94

せっかく作ったパッケージなのになんで公開やめちゃったんだろう
新しいバージョンで動かなくなってメンテやめちゃったのかな
0900１３２人目の素数さん2014/09/14(日) 15:11:45.59
30万円もするようなもの、そう簡単に買えるか。
0901１３２人目の素数さん2014/09/14(日) 22:00:53.73
学生だったら２、３万で買えるし、home editionが５万くらいでしょ？
0902１３２人目の素数さん2014/09/16(火) 18:26:30.32
https://codeiq.jp/ace/code_teller/q1086
「設問チラ見」で問題が見れる

この問題をMathematicaで解いて使用言語「Wolfram言語」として送ったった。

届いた評価フィードバックがこちら。

【解答内容】
最適解でした。

【コメント】
詳細な説明ありがとうございます。
解の全列挙も良いアピールですね。Wolfram言語の需要がどの程度あるかは予測できかねますが、良いオファーがあることをお祈り申し上げます。

【参考評価値（1〜5点）】
(1)提出フォーマット 5.0
(2)入出力 5.0
(3)可読性 5.0
(4)汎用性・拡張性 5.0
(5)判断力・思考力 5.0

総合評価 5.0

10行ぐらいのコードだったんだが、コードの内容には一切触れず。
なんかよくわからんけど合ってるよ、おめでとう、ぐらいの感じだな。
Mathematicaの実行環境がないということで実行すらしてないだろう。
Wolfram Programming Cloudの存在教えてやりゃーよかったかな。
0903１３２人目の素数さん2014/09/17(水) 18:14:38.85
>>902

挑戦済み0人になってるけど。
wolfram langageだとカウントもされんのか。
0904１３２人目の素数さん2014/09/17(水) 22:23:23.03
>>898
>立体を切断するにはどうすればいいんだ?
>
>knife.mは最近のバージョンでは使えないんだろ
>PlotRangeだと斜めとか切れないし


Ver.10で立体の共通部分とか
できるようになったからやってみそ
0905１３２人目の素数さん2014/09/17(水) 22:37:28.26
>>903

＞wolfram langageだとカウントもされんのか。
ワロタｗ
なんか受付期間を変えて何度もやってるみたい
https://codeiq.jp/ace/code_teller/
https://codeiq.jp/ace/code_teller/q986

この問題は言語不問だからWolfram言語でもOK
もちろん言語指定の問題もある

https://codeiq.jp/ace/code_teller/q1076

この問題をC#で解いたらオファーが3個くらいきた
0906１３２人目の素数さん2014/09/17(水) 22:42:07.67
>>904

立体の共通部分はまだダメだ
平面だけでした
0907１３２人目の素数さん2014/09/25(木) 21:03:59.71
>>900-901
今は一般人でも5000円もあれば使える時代に何を言ってるのやらｗ
ただしハードが遅いけどね
0908１３２人目の素数さん2014/09/27(土) 16:08:52.78
>>905だけどWolfram言語で解いた問題にもオファーきた
しかも2社

プログラミング言語総選挙にはC#で応募しといた
https://codeiq.jp/wp/vote/
0909１３２人目の素数さん2014/10/04(土) 08:03:17.17
いきなりですが、
素粒子論の専門家でmathematicaを使用している人には、どのような先生がいるのでしょうか？
0910１３２人目の素数さん2014/10/04(土) 13:23:50.72
増田 哲也 先生
0911１３２人目の素数さん2014/10/04(土) 13:34:58.54
＞910
有難う御座います。
更に、ご存じないでしょうか？
0912age2014/10/04(土) 14:25:46.19
物理板できけ
0913１３２人目の素数さん2014/10/04(土) 14:31:12.17
増田先生の専門は数理物理
0914１３２人目の素数さん2014/10/22(水) 11:24:33.26
パラメータ付けされた曲線の長さをマセマティカで計算する方法教えていただけませんか、、、
0915１３２人目の素数さん2014/10/22(水) 11:35:29.20
>>914
定積分の計算じゃないの
0916１３２人目の素数さん2014/10/22(水) 13:41:46.73
マセマティカで論文を生成する方法教えて
0917１３２人目の素数さん2014/10/23(木) 06:23:29.80
>>914
ArcLength
0918１３２人目の素数さん2014/10/27(月) 21:22:49.19
閉曲線が自分と交わってる点の求め方ってMathematicaでなにかありますか？
閉曲線のパラメータ表示はわかってるとします。
グラフみればなんとなくどこかわかるけど、グラフが複雑になるとわけわからなくて困ってます。
0919１３２人目の素数さん2014/11/09(日) 08:42:42.46
こんにちは、
添付は、電子・ミューオン散乱の計算です。
答えは間違っています。どこを修正すれば、良いでしょうか？

http://www.geocities.jp/dirac_equation/index1.html
0920１３２人目の素数さん2014/11/09(日) 22:43:38.91
>>919
y
0921１３２人目の素数さん2014/11/10(月) 19:56:11.98
ご回答有難う御座います。

>>920
ほんとですね。確かに、ｙの部分が、間違っておりました。

更に、電子陽電子の散乱を計算しました。
答えが全く合いません。
原因は、マンデルスタム変数の部分が、コンプトン散乱で使用していた
ものと異なるからだと思っております。
コンプトン散乱の場合は、電子はｍ、光子はｍ＝０
電子陽電子の散乱の場合は、両者ともｍとして計算するためです。

添付にコンプトン散乱の場合のマンデルスタム変数の計算をつけております。
s+t+u=2m^2
が導出出来ません。どこが間違っているのでしょうか？

http://www.geocities.jp/dirac_equation/index1.html
0922１３２人目の素数さん2014/11/10(月) 20:34:06.13
>>921
Mathematica文法で {a,b,c,d}^2 はどう計算されるか知ってる？
0923１３２人目の素数さん2014/11/10(月) 22:35:00.45
お返事有難う御座います。
>>922
{a,b,c,d}.{a,b,c,d}にすれば、良いわけですね。
全く、考えていなかったです。
電子・ミューオン散乱の計算は、相対論的極限で,m=M＝０のときは
答えが合っていますが、そうでない時は、導出できません。
電子、陽電子散乱も、陽電子の質量が考慮されておりませんので
正しい答えが導出出来ません。（多分）
要するに、下記の変換に、ミューオンのＭ（または陽電子のｍ）を反映させて
今度は、s+t+u=4m^2等の関係になるようにすることが必要です。
でも、そのＭ（またはｍ）を、（同じ場所のはずですが）、どこに追加すれば良いかは謎です。
直感的には、対称性からして
p0→(s+m^2+M^2)/(2Sqrt[s]),
p3→(s-m^2+M^2)/(2Sqrt[s]),
t→m^2+M^2-s-u
ではないかと、思っております。
0924１３２人目の素数さん2014/11/10(月) 22:45:43.70
>>923
＞ではないかと、思っております。
計算したら、全くはずれてました。(泣)
0925１３２人目の素数さん2014/11/10(月) 23:16:18.52
>>921
>>740
0926１３２人目の素数さん2014/11/11(火) 21:25:34.51
お返事有難う御座います。
>>740
承知致しました。
下記に
１．コンプトン散乱のマンデルスタム変数を満たす確認計算
２．電子・陽電子散乱のマンデルスタム変数を満たす確認計算
３．電子・陽電子散乱の計算
を入れております。

直感的は正しいかもしれません。（一部訂正）
p0→(s+m^2+M^2)/(2Sqrt[s]),
p3→(s-m^2-M^2)/(2Sqrt[s]),
t→２＊（m^2+M^2）-s-u

なぜか、３の計算が合いません。なぜでしょうか？

http://www.geocities.jp/dirac_equation/index1.html
0927１３２人目の素数さん2014/11/11(火) 21:40:48.10
>>926
>p0→(s+m^2+M^2)/(2Sqrt[s]),
>p3→(s-m^2-M^2)/(2Sqrt[s]),
この粒子の質量はどうなる？
0928１３２人目の素数さん2014/11/11(火) 23:42:50.79
>>926
>j0 -> p3, j1 -> -p3*Sqrt[1 - z^2], j2 -> 0, j3 -> -p3*z
この粒子の質量は？
0929１３２人目の素数さん2014/11/12(水) 00:48:11.77
>>927
お返事有難う御座います。

すいません。基本的なことを間違っていました。
この計算は、電子・陽電子ではなく、電子と電子の散乱でした。

＞この粒子の質量はどうなる？
全体的には４元運動量として、保存されます。質量だけに注目すると、ローレンツ変換により、速度に応じて、変化すると思います。
数値では解らないです。

mathematicaの20141112の「電子・電子散乱のマンデルスタム変数の確認　その１」の計算は間違っているでしょうか？いいところまで、満たしております。
１つだけ、“０”になるべきところ、“２ｍ＾２”になります。８０点です。


「電子・電子散乱のマンデルスタム変数の確認　その２」は、ｐとｋ（またqとｋ）で、質量を分担させましたが、結果は、全く、マンデルスタム変数の条件を満たさないです。
０点です。

どこから、考えれば良いでしょうか？
0930１３２人目の素数さん2014/11/12(水) 00:52:39.68
>>929
>どこから、考えれば良いでしょうか？
特殊相対論を勉強するところから。
0931１３２人目の素数さん2014/11/12(水) 19:42:48.56
mandelstam variables と mathematicaでググれば？
自力も大切ですが、すでに聞きまくりですから(^^;)
0932１３２人目の素数さん2014/11/12(水) 22:10:32.91
>>930
お返事有難う御座います。
P0^2-p1^2-p2^2-p^2=m^2
になるということですか。

>p0→(s+m^2+M^2)/(2Sqrt[s]),
>p3→(s-m^2-M^2)/(2Sqrt[s]),
この場合、質量がm^2+M^2になって、間違い。
すると
ｍ→ｍ＋Ｍ
に変更すればよろしいでしょうか？
0933１３２人目の素数さん2014/11/12(水) 22:49:48.89
>>932
>>857
0934１３２人目の素数さん2014/11/13(木) 23:11:44.57
お返事有難う御座います。
>>931
日本語文献は、ほぼ完璧に探しましたが、英文は探しませんでした。
英文の検索しました、まだ完全に探しておりませんが、ほしい式は見つからないです。


>>933
そうですか。
下記ＨＰの本の式（８２．２）は正しいでしょうか？
特に２番目の式が、怪しいのでは？と思っています。

http://www.geocities.jp/dirac_equation/index1.html
0935１３２人目の素数さん2014/11/13(木) 23:44:41.29
>>934
間違ってるとしたら、どこがどう間違っていて、正しくはどうだと思うのか書いてみ
0936１３２人目の素数さん2014/11/15(土) 14:00:23.02
こんにちは、
>>931
探しましたが、ないです。
下記の最後のプログラムが近いですが、今の計算とは少し違うので、使えないです。

http://www.phy.pmf.unizg.hr/~kkumer/articles/feynman_for_beginners.pdf#search='mandelstam+variables+++mathematica'

もし、参考になるような式があれば、教えて下さい。

>>935

http://en.wikipedia.org/wiki/Mandelstam_variables

を見ると、間違っていないことがわかりました。

ｐの電子とｋの電子が衝突して、ｑとjの状態になって、飛び去る。
ｐのp0^2-p1^2-p2^2-p3^3を計算すると“ｍ＾２”になり
ｋのそれは、”0”
qのそれは、”m＾２”
ｊのそれは、“０”
になります。
ｋが０になるのは、光子質量が０であるからであり、ｋが“ｍ＾２”になるように計算すれば
良いのでしょうか？

いろいろと試しましたが、わからないです。かなり、コンプトン散乱の場合と、変更があるのではないでしょうか？

http://www.geocities.jp/dirac_equation/index1.html
0937１３２人目の素数さん2014/11/15(土) 20:04:20.80
こんにちは、
http://www.geocities.jp/dirac_equation/index1.html
の
20141115の計算　その１
で正しい否かを教えて下さい。
0938１３２人目の素数さん2014/11/15(土) 22:34:19.17
毎度のことながら、正しい計算結果が出ません。
http://www.geocities.jp/dirac_equation/index1.html

20041115 電子・電子散乱
どこに、問題があるでしょうか？
0939１３２人目の素数さん2014/11/16(日) 08:29:42.99
重心系と実験室系を考えていませんでした。
答えが合わないはずです。
0940１３２人目の素数さん2014/11/16(日) 10:57:58.70
計算が合わない要因
１、電子・電子散乱のマンデルスタム変数を使った変換式が間違っている
　　（形からして、重心系の変換式のはず。82.8を代入したら一致するはず）
２、計算結果は正しいが、s,t,uの３種類の変数が混在しているため、見かけ上
異なるように見える
３、その他の箇所で、計算間違いをしている
等々、どこが悪いのか？その原因が解らないので、対策が立てられない状況です。
0941１３２人目の素数さん2014/11/16(日) 23:27:19.42
>>938
>どこに、問題があるでしょうか？
sl[j]の定義式より後は全部問題あり。「何一つ正しい式を書いてない」レベル。
0942１３２人目の素数さん2014/11/17(月) 21:00:48.13
ご指摘有難う御座います。
承知しました。計算をやり直してみます。
ｐ1,p2,p3,p4とｐ、ｑ、ｋ、ｊが、対応してないかもしれません。
0943１３２人目の素数さん2014/11/24(月) 10:03:39.52
こんにちは、
Compton散乱を計算し直しております。
下記ＨＰのCompton散乱は、ＨＰのHeitlerの本を元に計算してます。
プログラムComptonは、式（３９）に式（４）を入れて、式（４０）を導いてます。
プログラムCompton　R2は、行列を一部変更したのですが、式（４０）に該当する式に、ｍが残ってしまい、グラフ化しません。グラフ化することは可能でしょうか？
変数、k0,k,u,γに対して、式（４）とγ=k0/uの関係式しかないので、計算は無理でしょうか？

http://www.geocities.jp/dirac_equation/index1.html
0944１３２人目の素数さん2014/11/24(月) 19:05:30.29
>>943
ガンマ行列の性質も知らん人にはQEDの計算は無理
0945１３２人目の素数さん2014/11/25(火) 20:30:33.71
＞944
すいません。説明不足です。
プログラムCompton　R2のガンマ行列もどきは、ガンマ行列になっていません。
そもそも、プログラムCompton　R2は、間違っていますが、このような場合、最後に“ｍ”が残ると、変数γだけの式にはならないでしょうか？
変数、k0,k,u,γに対して、式（４）とγ=k0/u（γ=k0/ｍ）の関係式しかないので、計算は無理でしょうか？
0946１３２人目の素数さん2014/11/26(水) 16:53:01.39
>>945
最後の式をよく読めばどうしてmが残ってるか分かる
0947１３２人目の素数さん2014/11/26(水) 22:28:18.93
お返事有難う御座います。」
＞最後の式をよく読めばどうしてmが残ってるか分かる
ご指摘有難う御座います。確認します。（現状、よくわからないです）
＞特殊相対論を勉強するところから。
今、酔っています。現状で、仲間（賢い人一杯いると思います）に「特殊相対論」って一言、言ったら、「病院で見てもらえ！」って（目で）言われます（ました）。
でも、ワイスコップの人の心を埋めるのは、「科学」「宗教」「芸術」であるいう言葉を信じます。（確信しています。）実は、数歩、進んでいる。
あくまでも、物理「科学」大好き、（あほですが）よろしくお願いいたします。
0948１３２人目の素数さん2014/11/27(木) 12:22:31.02
＞最後の式をよく読めばどうしてmが残ってるか分かる
すいません。最後の式とは、ハイトラーの本の式（４０）のことでしょうか？
それとも（４３）でしょうか？
0949１３２人目の素数さん2014/11/27(木) 16:37:23.94
>>948
貴方が自分で出力して
>最後に“ｍ”が残る
と言ってる式
0950１３２人目の素数さん2014/11/27(木) 22:21:11.18
＞＞949
見て考えたのですが、わかりません。
ｍを消すことは可能でしょうか？
0951１３２人目の素数さん2014/11/27(木) 23:07:47.44
>>950
FullSimplify[(m^2)^(3/2)/m^3] はどうなるかをを考えてみ
0952１３２人目の素数さん2014/11/28(金) 07:17:52.06
＞＞951
有難う御座います。考えます。
0953１３２人目の素数さん2014/11/28(金) 21:41:47.33
お返事有難う御座います。
＞FullSimplify[(m^2)^(3/2)/m^3] はどうなるかをを考えてみ
ボキャブラリーというか、ExpnadAllで展開させることしか、思い浮かばないです。
まず、このFullSimplify[(m^2)^(3/2)/m^3]を、どうしたら、“１”になるのか？わからないです。
±があるので、正確には、これ以上どうしょうも、ないのでしょうか？
（ルートｍ＾２の解は、＋ｍと-ｍの２つが考えられる。）
このコンプトン散乱の場合、±の-は、一応、無視したいです。
その場合でも、どうしたら、良いのか？わからないです。
0954１３２人目の素数さん2014/11/28(金) 21:44:09.26
不細工ですが
√ｍ＾２→ｍ　に、力ずくで、置き換えれば良いかもしれませんね。
0955１３２人目の素数さん2014/11/28(金) 23:29:22.97
>>953
高校数学やり直し
0956１３２人目の素数さん2014/11/29(土) 10:50:53.66
ワロタ
0957１３２人目の素数さん2014/11/29(土) 10:51:45.15
FullSimplify[(m^2)^(3/2)/m^3, m > 0]
0958１３２人目の素数さん2014/11/29(土) 12:19:06.14
こんにちは、
ご回答有難う御座います。
今度は、図に出来ません。
原因は、虚数だからでしょうか？
（しかし、式全体に“Ｉ”を掛けても、図にできません。）
目盛りが、一致してないからでしょうか？

http://www.geocities.jp/dirac_equation/index1.html
0959１３２人目の素数さん2014/11/29(土) 13:22:18.28
>>958
この計算の答えは、存在しない（間違っている）
ということでしょうか
0960１３２人目の素数さん2014/11/30(日) 20:52:35.08
>>958
>>955
「複素数」って知ってる？
0961１３２人目の素数さん2014/11/30(日) 21:48:59.39
そうか。考え付かなかった。複素共役から、それなりのグラフ化が可能かも、しれませんね。
0962１３２人目の素数さん2014/11/30(日) 23:18:21.15
>>959
Mathematicaが出力したものが「この計算の答え」じゃないなら何なんだ
>>945によれば、そもそも入力が物理的にはデタラメだということは理解してるんだろ？
0963１３２人目の素数さん2014/12/01(月) 00:12:38.31
ノート見てないから適当言うけど、m って質量じゃないの？
0964１３２人目の素数さん2014/12/23(火) 15:31:03.39
こんにちは
乱数で表示された数のうち，1つだけ数えたいのですが・・・
1回1回実行するのが面倒なのでなんとかして，1度に複数回こなせるようなものを考えたいんですけどできません。
どうすればいいんでしょうか・・・？
0965１３２人目の素数さん2014/12/23(火) 19:33:03.97
>>964
>乱数で表示された数のうち，1つだけ数えたい
意味不明。それを1回実行するMathematiaコードを書いてみて。
0966１３２人目の素数さん2014/12/23(火) 22:21:46.85
このスレはコンプトンをやってる人に個人的に教えてあげるスレですか？
0967１３２人目の素数さん2014/12/24(水) 10:00:12.32
>>965
>意味不明
そりゃそうですよね・・・
すいませんでした。
一回実行するコードはこんな感じです。

dice010402 = Table[Random[Integer, {1, 4}],{22}]
dice010602 = Table[Random[Integer, {1, 6}], {3}]
count1[x_] := Which[x == 1, 1, x != 1, 0]
Sum[count1[dice010402[[k]]], {k, 1, 22}]
Sum[count1[dice010602[[k]]], {k, 1, 3}]

これを実行すると一回分の結果しか出ませんが，それを一度に複数回こなせるようなものに改良したいんです。
どうすればいいのでしょうか？
0968１３２人目の素数さん2014/12/24(水) 17:42:18.10
>>967
こういったかんじ？
Tally@RandomInteger[{1, 4}, 22]
Tally@RandomInteger[{1, 6}, 3]
0969１３２人目の素数さん2014/12/25(木) 13:58:44.35ID:U2/O3mDM
>>968
そんなかんじです
それを改良したいんですけど…
0970１３２人目の素数さん2014/12/25(木) 14:59:44.26ID:y7AcYpT7
>>969
改良とは、例題での1の度数だけでなく、総ての目の度数を
しらべることではないの？

968ではそうしてるけど･･･(出力フォーマットは別として)
0971１３２人目の素数さん2014/12/30(火) 15:20:41.68ID:yc9gAwbH
Jane Styleから書き込みえできません。
0972１３２人目の素数さん2014/12/30(火) 15:23:19.58ID:yc9gAwbH
>>971
書き込み出来ました。

こんにちは、
以下を教えて下さい。
１．ガンマ行列のトレース計算をしました。
http://www.geocities.jp/dirac_equation/index1.html

この計算で正しいのでしょうか？
ps={p0,-I*p1,-I*p2,-I*p3};
qs={q0,-I*q1,-I*q2,-I*q3};
ks={k0,-I*k1,-I*k2,-I*k3};
js={j0,-I*j1,-I*j2,-I*j3};
のように、空間成分に-Iが付くのでしょうか？

２．電子電子散乱の計算で、
http://www.geocities.jp/dirac_equation/landau.pdf
の式（82.6）の上の２つの式が、式（82.6）に一致することを計算しました。
http://www.geocities.jp/dirac_equation/index1.html

２つの式の内、
y01 = 2/t^2*((ps.ks)*(ps.ks) + (ps.js)*(ps.js) + 2*m^2*(m^2 - (ps.qs)))
y02 = 2/(t*u)*((ps.ks) - 2*m^2) *(ps.ks)
ｙ02は、式（82.6）に一致しますが、ｙ01は、一致しません。
どこが間違っているのでしょうか？
0973１３２人目の素数さん2014/12/30(火) 21:38:48.59ID:J3Yd0htO
特殊相対論の勉強を頑に拒むのはなぜ？
0974１３２人目の素数さん2015/02/11(水) 16:44:01.60ID:Yq7k8Ekq
てｓ
0975１３２人目の素数さん2015/02/11(水) 16:46:02.28ID:Yq7k8Ekq
てｓ
0976１３２人目の素数さん2015/02/11(水) 16:58:23.97ID:Yq7k8Ekq
てｓ
0977１３２人目の素数さん2015/02/17(火) 10:37:44.40ID:jGAH4sK+
1つのノートブック内に3つのセルを作成しました．

セルA：　計算パラメータの定義を集めたもの．
セルB：　計算パラメータの変更に依存しない部分の計算を集めたもの．
セルC：　計算パラメータに依存する計算を集めたもの．

この状況で，各セルを適宜再実行しながら計算するコードをセルDに書きたいと考えています．
そこで「あるセルから他のセルを呼び出して実行」ができるとよいのですが，可能でしょうか？
0978１３２人目の素数さん2015/02/17(火) 17:35:27.68ID:EdiRDBQD
>>977
Block とか Module とかじゃいかんの？
0979１３２人目の素数さん2015/02/18(水) 18:18:02.68ID:Af217Zem
>>978
うまく出来そうです．
ありがとうございます．
0980１３２人目の素数さん2015/02/23(月) 02:56:52.76ID:Jtdu54/Z
Needs["Benchmarking`"];
BenchmarkReport[]

1.33だった
0981１３２人目の素数さん2015/02/23(月) 18:50:29.06ID:x7brCvuk
>>980
System: Microsoft Windows (64-bit)
Date: February 23, 2015
Mathematica Version: 10.0.2
Benchmark Result: 1.95
0982１３２人目の素数さん2015/02/23(月) 19:05:54.69ID:Jtdu54/Z
どんだけいいマシン使ってんだよｗ
比較用に列挙されてるやつの一番いいやつでも1.89じゃん

マシンのスペック書いて！
最近は定格で4GHz超えてるマシンも出てきたみたいやね
0983１３２人目の素数さん2015/02/23(月) 22:26:55.60ID:x7brCvuk
>>982
CPU: Intel(R) Core(TM) i7-3970X CPU @ 3.50GHz
RAM: 64.0 GB
0984１３２人目の素数さん2015/02/23(月) 22:31:26.42ID:x7brCvuk
CPU: Intel(R) Xeon(R) CPU E5-2690 v2 @ 3.00GHz
RAM: 256 GB

System: Microsoft Windows (64-bit)
Date: February 23, 2015
Mathematica Version: 10.0.1
Benchmark Result: 1.41



明らかに得手不得手が出てるな
0985１３２人目の素数さん2015/02/23(月) 23:56:40.43ID:2NtpOD86
>>983
自分のは
Intel Core i7
Processor Speed: 2.7 GHz
RAM 16GB
だけどなんで64GBもメモリを積めるんですか？
0986１３２人目の素数さん2015/02/24(火) 00:22:09.26ID:1T2TwbB0
>>985
ﾏｻﾞﾎﾞが違うと思われ
俺が使ってるやつ↓
http://www.asus.com/jp/Motherboards/P9X79/
0987１３２人目の素数さん2015/02/24(火) 02:02:33.70ID:XH+kD/gH
メモリ256GBは初めて聞いた
CPUクロック数の上昇もそろそろ限界に達してるから
世の中のHPCの動向がそうであるように
MathematicaもそろそろGPUへの対応をした方がいいね
CUDAコード出力するだけじゃなくて組み込みレベルで。
0988１３２人目の素数さん2015/02/24(火) 20:16:28.04ID:9BvsOY5B
>>984
Mathematicaのデフォは4coreがmaxでなかったっけ？、
だとしたら10coreに上げてもclockを下げた分不利…

core extension導入でcore数を増やせるらしいけど、
ベンチマーク上は微増らしい、つか明示的に並列ｐ
書かんと速度上がらんぽい(コア管理がユーサー側

これじゃMathematicaを使う意味ないじゃんと想う
0989１３２人目の素数さん2015/02/25(水) 00:49:54.91ID:4QJ9DHPX
一部GPU対応しているはずだけどね
ちなみに近頃のHPC
ttp://www.hpc.co.jp/HPC7000-XI8192R8P_features.html
0990１３２人目の素数さん2015/02/25(水) 09:02:28.35ID:GacQMgDm
三年百六日。
0991１３２人目の素数さん2015/02/25(水) 20:51:50.18ID:GA/9W6kV
DateDifference[{2011, 11, 11}, DateValue[{"Year", "Month", "Day"}]]

->1202 days(*平均1ﾚｽ/day以下*)
0992１３２人目の素数さん2015/02/26(木) 00:18:27.64ID:UjfUV2xI
今はMathematica10だからスレッドもpart10になるぐらいのスピードにしないとな
1バージョン1スレッド
0993１３２人目の素数さん2015/02/26(木) 02:32:45.39ID:UjfUV2xI
このスレにはいないと思うが、これできる人いたらお願い。
http://wc2014.2ch.net/test/read.cgi/math/1424876583/
レス数が950を超えています。1000を超えると書き込みができなくなります。