24СʱÈÈÃŰæ¿éÅÅÐаñ    

±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1554  |  »Ø¸´: 9
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

dollapym

Ìú³æ (ÕýʽдÊÖ)

[ÇóÖú] ¼±¼±¼±£¡£¡£¡ÇóÖúÓÃBICG»òÕßrestarted GMRES ½âÏßÐÔ·½³Ì×éA x=b µÄMATLAB code ÒÑÓÐ1È˲ÎÓë

Çó´óÉñ°ïæ°¡ £¬ ºÜ¼±Èç¹ûÐèÒªµÄ»°¿ÉÒÔÔÙ¼Ó½ð±ÒµÄ£¡£¡
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

dollapym

Ìú³æ (ÕýʽдÊÖ)

ÒýÓûØÌû:
7Â¥: Originally posted by FMStation at 2016-09-07 16:56:43
https://github.com/exasc/lssp
LSSP is a Linear Solver package for SParse linear system, Ax = b.

ллÄãÌṩµÄ×ÊÁÏ£¬ËäÈ»²»ÄÜÖ±½Ó½â¾öÎÒµÄÎÊÌ⣬µ«ÊÇÈÃÎÒÁ˽âÁËÐí¶à֪ʶ
8Â¥2016-09-08 14:42:32
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 10 ¸ö»Ø´ð

FMStation

ÖÁ×ðľ³æ (ÖªÃû×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
dollapym: ½ð±Ò+10, ¡ïÓаïÖú 2016-09-07 10:31:01
dollapym: ½ð±Ò+10, ¡ïÓаïÖú 2016-09-07 10:37:42
https://www.mathworks.com/help/matlab/ref/bicg.html
CODE:
load west0479;
A = west0479;
b = full(sum(A,2));
tol = 1e-12;
maxit = 20;
[x0,fl0,rr0,it0,rv0] = bicg(A,b,tol,maxit);
semilogy(0:maxit,rv0/norm(b),'-o');
xlabel('Iteration number');
ylabel('Relative residual');

2Â¥2016-09-06 20:26:01
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

dollapym

Ìú³æ (ÕýʽдÊÖ)

ÒýÓûØÌû:
2Â¥: Originally posted by FMStation at 2016-09-06 20:26:01
https://www.mathworks.com/help/matlab/ref/bicg.html

load west0479;
A = west0479;
b = full(sum(A,2));
tol = 1e-12;
maxit = 20;
= bicg(A,b,tol,maxit);
semilogy(0:maxit,rv0/norm(b),'-o');
xla ...

²»ÊÇÖ±½ÓÓú¯Êýbicg£¬ÊÇҪдÕâ¸ö·½·¨µÄ´úÂë

·¢×ÔСľ³æIOS¿Í»§¶Ë
3Â¥2016-09-06 20:55:04
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

FMStation

ÖÁ×ðľ³æ (ÖªÃû×÷¼Ò)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

http://www.aip.de/groups/soe/local/numres/bookcpdf/c2-7.pdf
void linbcg(unsigned long n, double b[], double x[], int itol, double tol,
int itmax, int *iter, double *err)
4Â¥2016-09-06 21:46:35
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸ Î÷±±´óѧ ×Ü·Ö282 Ó¢ÓïÒ»62 Çóµ÷¼Á +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[¿¼ÑÐ] 279 ·Ö Çóµ÷¼Á +4 ˯¸öºÃ¾õ_16 2026-03-24 4/200 2026-03-27 15:05 by ×íÔÚ·çÀï
[¿¼ÑÐ] 08¿ªÍ·275Çóµ÷¼Á +4 À­Ë­²»ÖØÒª 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[¿¼ÑÐ] »¯Ñ§308·ÖÇóµ÷¼Á +8 ÄãºÃÃ÷ÌìÄãºÃ 2026-03-23 9/450 2026-03-27 14:01 by Ñî¹âÓÚÇàÔÆ
[¿¼ÑÐ] Ò»Ö¾Ô¸Ö£´ó085600£¬310·ÖÇóµ÷¼Á +5 Àîäì¿É 2026-03-26 5/250 2026-03-27 11:14 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] ¿¼Ñе÷¼Á +9 СÀ¯ÐÂ±Ê 2026-03-26 9/450 2026-03-27 11:10 by ²»³Ôô~µÄ؈
[˶²©¼ÒÔ°] ±±¾©ÁÖÒµ´óѧ˶µ¼ÕÐÉú¹ã¸æ +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
[¿¼ÑÐ] 359Çóµ÷¼Á +4 ÍõÁ˸öéª 2026-03-25 4/200 2026-03-27 08:43 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 321Çóµ÷¼Á +6 Ymlll 2026-03-24 6/300 2026-03-26 20:50 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 281Çóµ÷¼Á +3 ÑÇ¿ËÎ÷good 2026-03-26 5/250 2026-03-26 19:48 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] ²ÄÁÏ¿¼ÑÐÇóµ÷¼Á +3 Dendel 2026-03-23 6/300 2026-03-26 17:51 by fmesaito
[¿¼ÑÐ] ²ÄÁÏ¿ÆÑ§Ó빤³Ì 317Çóµ÷¼Á +4 JKSOIID 2026-03-26 4/200 2026-03-26 15:58 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 263Çóµ÷¼Á +6 yqdszhdap£­ 2026-03-22 10/500 2026-03-26 13:11 by ¹«èªåÐÒ£
[¿¼ÑÐ] ÉúÎïѧѧ˶Çóµ÷¼Á +7 СÑò˯×ÅÁË? 2026-03-23 10/500 2026-03-25 02:24 by Çå·ç·÷Ñï¡£ m
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á£¬Çóµ¼Ê¦ÊÕ +7 ÌìÌìºÃÔËÀ´Éϰ¶° 2026-03-24 7/350 2026-03-24 20:26 by peike
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄϺ½²ÄÁÏר317·ÖÇóµ÷¼Á +5 ըѽըѽըÊíÌõ 2026-03-23 5/250 2026-03-24 16:52 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 305·ÖÇóµ÷¼Á£¨Ê³Æ·¹¤³Ì£© +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖØÇì´óѧ085700×ÊÔ´Óë»·¾³£¬×Ü·Ö308Çóµ÷¼Á +7 īīĮ 2026-03-23 8/400 2026-03-23 20:36 by Creta
[¿¼ÑÐ] ²ÄÁÏѧ˶301·ÖÇóµ÷¼Á +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +4 ÄÝÄÝninicgb 2026-03-21 4/200 2026-03-21 18:39 by ѧԱ8dgXkO
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û