On Jul 19, 2012, at 4:28 AM, Nicolas TARDIEU wrote: > > Thanks for your answer, Matt. > The problem is that I am using ML as a preconditioner and it does not support > 64 bit indices (information reported in the PETSc configure phase). > I am afraid I am stuck with this wrong nonzeros number.
You cannot use ML for such large problems, if ML it is only built for 32 bit integers. It is that simple. Barry > > Nicolas > <Mail Attachment.gif> > Nicolas TARDIEU > Ing. Chercheur > EDF - R&D Dpt AMA > > nicolas.tardieu at edf.fr > T?l. : 01 47 65 39 05 > <Mail Attachment.gif> Un geste simple pour l'environnement, > n'imprimez ce message que si vous en avez l'utilit?. > > > > > knepley at gmail.com > Envoy? par : petsc-users-bounces at mcs.anl.gov > 17/07/2012 16:59 > Veuillez r?pondre ? > petsc-users at mcs.anl.gov > > A > petsc-users at mcs.anl.gov > cc > Objet > Re: [petsc-users] Number of nonzeros > 2^31-1 > > > > > > On Tue, Jul 17, 2012 at 6:20 AM, Nicolas TARDIEU <nicolas.tardieu at edf.fr> > wrote: > > Dear PETSc users, > > I am solving an unstructured finite element-based problem with 500 millions > unkwowns with PETSc. > The number of nonzeros is greater than 2^31-1. > Here is the KSPView I get : > ---------------------------------------------------------------------------------------------------------------- > > Matrix Object: 700 MPI processes > type: mpiaij > rows=499125000, cols=499125000 > total: nonzeros=-2147483648, allocated nonzeros=-2147483648 > total number of mallocs used during MatSetValues calls =0 > using I-node (on process 0) routines: found 242923 nodes, limit used is > 5 > ---------------------------------------------------------------------------------------------------------------- > > As you can see, The number of nonzeros is <0. > I would like to check that this is due to the number of nonzeros being > greater than 2^31-1. > > Here is a short description of the size of differents types : > ---------------------------------------------------------------------------------------------------------------- > > Compiled with full precision matrices (default) > sizeof(short) 2 sizeof(int) 4 sizeof(long) 8 sizeof(void*) 8 > sizeof(PetscScalar) 8 > ---------------------------------------------------------------------------------------------------------------- > > > Is there a workaround? Should I worry about the result of my simulation? > > You could configure using --with-64-bit-indices. > > Matt > > Thanks in advance, > Nicolas > <Mail Attachment.gif> > Nicolas TARDIEU > Ing. Chercheur > EDF - R&D Dpt AMA > > nicolas.tardieu at edf.fr > T?l. : 01 47 65 39 05 > <Mail Attachment.gif> Un geste simple pour l'environnement, > n'imprimez ce message que si vous en avez l'utilit?. > > > Ce message et toutes les pi?ces jointes (ci-apr?s le 'Message') sont ?tablis > ? l'intention exclusive des destinataires et les informations qui y figurent > sont strictement confidentielles. Toute utilisation de ce Message non > conforme ? sa destination, toute diffusion ou toute publication totale ou > partielle, est interdite sauf autorisation expresse. > > Si vous n'?tes pas le destinataire de ce Message, il vous est interdit de le > copier, de le faire suivre, de le divulguer ou d'en utiliser tout ou partie. > Si vous avez re?u ce Message par erreur, merci de le supprimer de votre > syst?me, ainsi que toutes ses copies, et de n'en garder aucune trace sur > quelque support que ce soit. Nous vous remercions ?galement d'en avertir > imm?diatement l'exp?diteur par retour du message. > > Il est impossible de garantir que les communications par messagerie > ?lectronique arrivent en temps utile, sont s?curis?es ou d?nu?es de toute > erreur ou virus. > ____________________________________________________ > > This message and any attachments (the 'Message') are intended solely for the > addressees. The information contained in this Message is confidential. Any > use of information contained in this Message not in accord with its purpose, > any dissemination or disclosure, either whole or partial, is prohibited > except formal approval. > > If you are not the addressee, you may not copy, forward, disclose or use any > part of it. If you have received this message in error, please delete it and > all copies from your system and notify the sender immediately by return > message. > > E-mail communication cannot be guaranteed to be timely secure, error or > virus-free. > > > > > -- > What most experimenters take for granted before they begin their experiments > is infinitely more interesting than any results to which their experiments > lead. > -- Norbert Wiener > > > > Ce message et toutes les pi?ces jointes (ci-apr?s le 'Message') sont ?tablis > ? l'intention exclusive des destinataires et les informations qui y figurent > sont strictement confidentielles. Toute utilisation de ce Message non > conforme ? sa destination, toute diffusion ou toute publication totale ou > partielle, est interdite sauf autorisation expresse. > > Si vous n'?tes pas le destinataire de ce Message, il vous est interdit de le > copier, de le faire suivre, de le divulguer ou d'en utiliser tout ou partie. > Si vous avez re?u ce Message par erreur, merci de le supprimer de votre > syst?me, ainsi que toutes ses copies, et de n'en garder aucune trace sur > quelque support que ce soit. Nous vous remercions ?galement d'en avertir > imm?diatement l'exp?diteur par retour du message. > > Il est impossible de garantir que les communications par messagerie > ?lectronique arrivent en temps utile, sont s?curis?es ou d?nu?es de toute > erreur ou virus. > ____________________________________________________ > > This message and any attachments (the 'Message') are intended solely for the > addressees. The information contained in this Message is confidential. Any > use of information contained in this Message not in accord with its purpose, > any dissemination or disclosure, either whole or partial, is prohibited > except formal approval. > > If you are not the addressee, you may not copy, forward, disclose or use any > part of it. If you have received this message in error, please delete it and > all copies from your system and notify the sender immediately by return > message. > > E-mail communication cannot be guaranteed to be timely secure, error or > virus-free. >
