02-14-2019
Quote:
Originally Posted by
RavinderSingh13
Wowww, Woww, Wowww. This is simply AWESOME. You ROCK Neo.
Thanks,
R. Singh
LOL... yes, I know.... ROTFL
Thanks for the compliment, Ravinder.
This solution has a long way to go yet. I'm mocking up the template now, step-by-step and coding new features as I go.
This User Gave Thanks to Neo For This Post:
6 More Discussions You Might Find Interesting
1. UNIX for Dummies Questions & Answers
I'm trying to setup a test environment to test upgrading to Solaris 8 or 9. Though when I go to install the demo version of Sun1 (forte/sun workshop as named before) it'll still ask for a 20 hexidecminal password for just a demo. huh?
Ok so I added a 20 character thing in there it passed... (2 Replies)
Discussion started by: merlin
2 Replies
2. UNIX for Dummies Questions & Answers
I have been thinking about another shell scripting project using the Arduino Diecimila board.
I was going to make a kids level slow 8 channel Logic Analyser.
I thought about using the termiinal esc code graphics characters.
This is the test code using said terminal escape codes. I have only... (2 Replies)
Discussion started by: wisecracker
2 Replies
3. What is on Your Mind?
Was working on Vue.js and stumbled upon this beautiful Vue project by Kadin Zhang
Periodicity is a dynamic periodic table built with Vue.js that animates and graphs data to aid the visualization of chemical concepts.
The code is available on GitHub (2 Replies)
Discussion started by: Neo
2 Replies
4. What is on Your Mind?
Regarding the latest version of the UserCP prototype (version 0.63) I have made a lot of major changes, including
Added a "Posts Timeline" table for the recent posts, complimenting the non-table version earlier, which has been moved off the main menu (link at the bottom of the table).
Added a... (4 Replies)
Discussion started by: Neo
4 Replies
5. What is on Your Mind?
Update!
UserCP Screeching Frog 0.7485
Created a new page for uploaded a profile picture (profile pictures are different than avatar pictures).
https://www.unix.com/usercp/#/settings/other
https://www.unix.com/members/1-albums225-picture1158.png
... (0 Replies)
Discussion started by: Neo
0 Replies
6. What is on Your Mind?
Here is my second live video screencast (in my life, using Camtasia) with voice for the new usercp:
Overview of the Vue.js UserCP @UNIX.com
Shout outs to Don Cragun, Corona688, Rudi, Wolf, Made in Germany, stomp, Ravinder, Creative Tim, PubNub and others in the video. Thanks.
If you are... (1 Reply)
Discussion started by: Neo
1 Replies
LEARN ABOUT DEBIAN
puzawa
puzawa(4rheolef) rheolef-6.1 puzawa(4rheolef)
NAME
puzawa -- Uzawa algorithm.
SYNOPSIS
template <class Matrix, class Vector, class Preconditioner, class Real>
int puzawa (const Matrix &A, Vector &x, const Vector &b, const Preconditioner &M,
int &max_iter, Real &tol, const Real& rho, odiststream *p_derr=0);
EXAMPLE
The simplest call to 'puzawa' has the folling form:
size_t max_iter = 100;
double tol = 1e-7;
int status = puzawa(A, x, b, EYE, max_iter, tol, 1.0, &derr);
DESCRIPTION
puzawa solves the linear system A*x=b using the Uzawa method. The Uzawa method is a descent method in the direction opposite to the gradi-
ent, with a constant step length 'rho'. The convergence is assured when the step length 'rho' is small enough. If matrix A is symmetric
positive definite, please uses 'pcg' that computes automatically the optimal descdnt step length at each iteration.
The return value indicates convergence within max_iter (input) iterations(0), or no convergence within max_iter iterations(1). Upon suc-
cessful return, output arguments have the following values:
x approximate solution to Ax = b
max_iter
the number of iterations performed before the tolerance was reached
tol the residual after the final iteration
IMPLEMENTATION
template < class Matrix, class Vector, class Preconditioner, class Real, class Size>
int puzawa(const Matrix &A, Vector &x, const Vector &Mb, const Preconditioner &M,
Size &max_iter, Real &tol, const Real& rho,
odiststream *p_derr, std::string label)
{
Vector b = M.solve(Mb);
Real norm2_b = dot(Mb,b);
Real norm2_r = norm2_b;
if (norm2_b == Real(0)) norm2_b = 1;
if (p_derr) (*p_derr) << "[" << label << "] #iteration residue" << std::endl;
for (Size n = 0; n <= max_iter; n++) {
Vector Mr = A*x - Mb;
Vector r = M.solve(Mr);
norm2_r = dot(Mr, r);
if (p_derr) (*p_derr) << "[" << label << "] " << n << " " << sqrt(norm2_r/norm2_b) << std::endl;
if (norm2_r <= sqr(tol)*norm2_b) {
tol = sqrt(norm2_r/norm2_b);
max_iter = n;
return 0;
}
x -= rho*r;
}
tol = sqrt(norm2_r/norm2_b);
return 1;
}
rheolef-6.1 rheolef-6.1 puzawa(4rheolef)