On Nov 20, 2007 10:23 AM, Viktor Cerovski <[EMAIL PROTECTED]> wrote:
> I think there is a problem with GS1.

Oops!  My mistake.

What you really want is something like (unwrap the following):

98 77 58 26 43 (] , [EMAIL PROTECTED]) 11 17 80 56 43 (] , [EMAIL PROTECTED]) 
30 18
77  6 49 (] , [EMAIL PROTECTED]) 49 32 47 75 17 (] , [EMAIL PROTECTED]) ,:norm 
76 50
56 90  6

The first norm operation can be accomplished with
  [EMAIL PROTECTED]:]`_1:`]}~ mat

The recurrence relationship can be accomplished with
   (] , [EMAIL PROTECTED])/

But you need to ensure that the right argument has a rank of 2, and the
you will start with a rank 1 right argument.  While there should be more
elegant ways of turning the initial vector into a "row vector", this will do:
   ,:^:(2 - [EMAIL PROTECTED])

Thus, I should have posted something like
   GS3=: [:   ((] , [EMAIL PROTECTED]) ,:^:(2 - [EMAIL PROTECTED]))/   [EMAIL 
PROTECTED]:@]`_1:`]}~

My GS1 and GS2 were totally wrong, because they provided vector
left arguments to orth, and not matrix arguments.  Partially, I did
not fully understand the domain, but also I did not read your code
carefully enough.

Thanks,

-- 
Raul
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm

Reply via email to