Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]


Groups > comp.soft-sys.math.mathematica > #3311

Re: Replacement rules for large matrixes

From Gabriel Landi <gtlandi@gmail.com>
Newsgroups comp.soft-sys.math.mathematica
Subject Re: Replacement rules for large matrixes
Date 2011-06-26 10:29 +0000
Organization Steven M. Christensen and Associates, Inc and MathTensor, Inc.
Message-ID <iu71ma$1bd$1@smc.vnet.net> (permalink)
References <201106250927.FAA20147@smc.vnet.net>

Show all headers | View raw


Dear Daniel,

Thank you for the reply. My matrix is constructed from a system of
polynomial equations using CoefficientArrays.
Indeed, replacement does not work. What I was doing was (sadly)

SparseArray[Normal@matrix /.rules].

Gabriel

On Sat, Jun 25, 2011 at 6:27 AM, dr DanW <dmaxwarren@gmail.com> wrote:

> Is your matrix a nested list or a SparseArray[]?  From your description, I
> am guessing it is a nested list.  Try working with your matrix as a
> SparseArray[] instead.
>
> Next problem: Replace[] does not work on SparseArray[]'s.  However, you can
> use BlockRules[] (code below) to do the replacement.
>
> Code for BlockRules[] by Andrew Moylen at Wolfram, from an earlier posting
> on this group that I can't find right now:
>
> SetAttributes[BlockRules, HoldRest];
>
> BlockRules[rules_, expr_] := Block @@ Append[Apply[Set, Hold[rules], {2}],
>    Unevaluated[expr]]
>
> Enjoy,
> Daniel
>
>

Back to comp.soft-sys.math.mathematica | Previous | Next | Find similar | Unroll thread


Thread

Re: Replacement rules for large matrixes Gabriel Landi <gtlandi@gmail.com> - 2011-06-26 10:29 +0000

csiph-web