[petsc-users] Create identity IS

[Florian Lindner]

Hello,

I use an index set for a row mapping but don't want to use one for the column mapping.

I try to create an idenitity index set to supply to ISLocalToGlobalMappingCreateIS and MatSetLocalToGlobalMapping:

  IS is;
  ISCreate(PETSC_COMM_WORLD, &is);
  ISSetIdentity(is);

but I get an segmentation fault at the last line.

What is the best way to create an identity index set?

Thanks,
Florian